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ABSTRACT 


Vodrazka,  Walter  Charles.   Ph.D.,  Purdue  University, 
June  196  8.   SUBCLASS IFI CATION  OF  THE  STATE  HIGHWAY  SYSTEM  OF 
INDIANA  BASED  ON  SYNTHESIS  OF  INTERCITY  TRAVEL.   Major 
Professor:   Harold  L.  Michael. 


This  research  report  is  concerned  with  the  development 
of  a  procedure  to  select  the  important  highways  of  the  fu- 
ture in  Indiana.   The  principal  aims  of  the  study  v;ere  to 
subclassify  the  State  Highv/ay  System  of  Indiana  and  to 
designate  those  highways  which  should  be  considered  for  re- 
construction to  freeway  standards. 

The  method  selected  to  achieve  these  aims  was  a  state- 
wide study  of  intercity  travel  desire  designed  to  synthesize 
travel  patterns  on  the  State  Highway  System  of  Indiana  as 
they  would  exist  were  the  Interstate  system  of  highways 
complete  and  fully  operational.   A  completed,  fully  opera- 
tional Interstate  system  was  thus  a  major  assumption  of 
this  study. 

The  Intercity  Travel  Desire  Factor  model  chosen  for 
this  study  was  the  product  of  the  square  roots  of  the  popu- 
lations of  two  interacting  cities  divided  by  the  square  of 
the  minimum  path  distance  between  the  cities.   The  resulting 
number  for  each  city  pair  was  assigned  to  each  highway  link 
making  up  the  minimum  path  between  the  cities.   A  cumulative 


total  was  maintained  for  each  Indiana  highway  link  in  the 
network  as  each  city  pair  was  considered. 

The  final  factor  associated  with  each  highway  link  was 
considered  to  be  a  measure  of  both  the  relative  importance 
of  the  link  and  of  the  anticipated  traffic  volume  on  the 
link.   More  than  300,000  city  interactions  were  calculated 
in  determining  the  factors  that  eventually  were  used  in 
achieving  the  study  objectives. 

The  largo  size  of  the  coded  highway  network  necessitated 
the  development  of  a  tree  type  decomposition  algorithm  for 
minimum  paths  m  large  networks.   In  addition  to  enabling 
a  solution  of  large  networks,  use  of  the  algorithm  resulted 
in  a  savings  of  computer  time  as  well. 

The  adequacy  of  the  Intercity  Travel  Desire  Factor  to 
synthesize  travel  was  determined  by  a  centroid  analysis  which 
consisted  of  a  regression  of  the  sum  of  the  AADT ' s  for  all 
links  entering  a  city  on  various  measures  of  the  sum  of  the 

factors  for  these  same  links.   The  model  was  accepted  as 

2 
adequate  on  the  basis  of  an  R   of  0.872  for  data  from  390 

Indiana  centroids. 

A  link  analysis  v.'as  then  performed  to  develop  a  rela- 
tionship between  the  link  volume  and  the  link  factor.   This 
was  determined  by  a  regression  of  the  minimum  AADT  for  a 
link  on  the  factor  associated  with  a  link.   The  relationship 

developed,  which  expressed  link  volum:e  as  a  function  of  the 

2 
logarithm  of  the  link  factor,  had  an  R   of  0.919. 


XI 


This  relationship  coupled  with  other  considerations  such 
as  system  connectivity  was  used  to  subclassify  the  State 
Highway  System  of  Indiana  into  four  subsystems:   the 
Principal,  Primary,  Secondary,  and  Collector  State  Highv/ay 
Systems . 

The  Principal  Highway  System  consisted  of  the  presently 
designated  Interstate  Highway  System  and  all  highways  which 
should  be  considered  for  reconstruction  to  freeway  standards 
by  1982. 


HIGHWAY  CLASSIFICATION 


Definition 


Functional  highway  classification  is  defined  as  the 
grouping  of  roads  and  streets  into  classes,  or  systems,  ac- 
cording to  the  character  of  service  they  v;ill  be  expected 
* 

to  provide  (  3)  . 

Many  other  definitions  of  highway  classification  may 
be  found  ranging  from  the  very  simple  such  as  highway 
classification  is  the  process  of  grouping  roads  and  streets 
according  to  the  requirements  of  traffic  (12)  to  more  com- 
prehensive statements  of  policy  such  as  the  following: 

The  functional  classification  of  public  roads 
is  the  procedure  by  '.vhich  portions  of  the  road 
and  street  network  of  a  State  are  grouped  ac- 
cording to  their  predominant  function  or  purpose 
in  the  overall  system,  and  the  responsibility 
for  each  grouping  is  determined  and  assigned  to 
the  level  of  government  having  the  predominant 
interest  in  providing  for  their  function.  (17) 

Some  of  the  desirable  characteristics  of  highway 

classification  procedures  and  results  are  contained  in  the 

definitions  employed  by  other  agencies:   the  highways  in  each 

system  or  grouping  must  have  similar  functional  usage  (4,  7, 

15,  16)  and  render  comparable  service  (15);  the  number  of 


Numbers  in  parentheses  refer  to  entries  in  th*-  List  of 
References . 


systems  into  which  roads  are  grouped  should  be  minimal  (16) 
but  the  size  of  each  system  should  be  manageable  (14);  and 
all  systems  are  to  be  systematically  interconnected  through- 
out the  areas  they  serve  (7).   The  definitions  of  two 
agencies  state  only  that  classification  is  the  process  by 
which  responsibility  for  each  mile  of  road  and  street  is 
assigned  to  one  unit  of  government  or  other  (1,  8) . 

This  last  definition  is  probably  the  best  one,  if  a  set 
of  criteria  on  which  the  assignment  of  responsibility  can  be 
made  is  available,  because  the  advantages  and  benefits  de- 
rived from  proper  highway  classf ication  accrue  to  the  govern- 
mental units  in  their  efforts  to  provide  adequate  highway 
facilities  for  the  public. 

Purpose  and  Necessity 

Highway  classification  has  been  described  as  the  basis, 
framework,  foundation,  and  first  requisite;  as  a  device  and 
tool;  and  as  being  indispensible .   It  makes  possible,  pro- 
vides for,  aids,  enables,  stimulates,  promotes,  im.proves, 
facilitates,  decreases,  increases,  and  has  advantages  and 
benefits . 

The  fundamental  purpose  of  proper  highway  classifica- 
tion is  that  responsibility  for  a  particular  system  be 
assigned  to  that  jurisdiction  most  interested  in  the  service 
the  system  provides. 

Each  jurisdiction  has  two  responsibilities:   to  best 
serve  the  economy  of  the  state  and  to  best  serve  the  public 


in  their  need  for  travel  to  places  of  employment,  to  haul 
livestock  and  farm  produce,  to  deliver  goods  and  mail,  to 
reach  tourist  attractions,  markets,  shopping  centers,  as 
well  as  business,  industrial,  and  social  destinations. 

There  are  many  reasons  why  classification  is  necessary. 
Not  all  highways  are  alike  either  in  respect  to  the  service 
provided  or  in  their  desirable  construction  and  maintenance 
standards.   Available  funds  are  limited.   For  example,  the 
average  expenditure  per  mile  of  rural  highway  in  South 
Dakota  in  1959  was  $700  (13).   Obviously  little  could  be 
accomplished  if  funds  were  allocated  equally  per  mile  of 
road.   Classification  provides  a  means  of  determining  the 
important  routes  upon  which  funds  should  be  concentrated  and 
the  extent  of  the  improvements  thereon. 

Classification  is  also  necessary  for  the  following 
reasons : 

1.  To  formulate  a  long  range  program  of  development 
which  is  orderly,  efficient,  and  economical  (10,  12,  16, 
18)  . 

2.  To  provide  for  a  systematic,  interlocked  arrange- 
ment of  finance,  management,  construction,  and  maintenance 

(9)  . 

3.  To  develop  sound  management  programs  (18). 

4.  To  provide  a  firm,  equitable  fiscal  structure  (12, 
18) .   The  main  problem  here  is  that  all  levels  of  government 
want  a  larger  share  of  the  available  funds  (13). 


5.  To  evaluate  the  relative  importance  of  highways  as 
a  basis  for  highway  improvement  programs  (19). 

6.  To  provide  for  modern  traffic  requirements  which 
demand  a  functional  highway  system  (25). 

Principles  and  Objectives 
Three  basic  principles  of  highway  classification  are 
(15)  : 

1.  A  road  or  street  section  can  be  placed  in  only  one 
system. 

2.  The  authority  for  a  system  can  be  vested  in  one 
government  agency  such  as  the  state  or  in  government  agen- 
cies of  the  same  level  such  as  the  several  counties  of  a 
state. 

3.  Each  government  agency  is  completely  responsible 
for  the  management,  safety,  im.provement,  operation,  and 
maintenance  of  all  roads  in  its  jurisdiction. 

The  systems  of  roads  resulting  from  the  application  of 
a  classification  technique  must  meet  several  objectives. 
Among  these  are: 

1.   Service .   The  needs  of  traffic  must  be  served  by 
each  route  through  systems  integrated  for  mutual  support 
to  provide  service  between  all  areas  of  the  state  (10,  18). 
Most  travel  is  done  over  a  road  net^vork  or,  more  specifically, 
on  a  number  of  different  roads  such  that  fev;  roads  serve 
travel  independently  of  others.   A  characteristic  of  modern 
travel  is  that  all  traffic  tends  naturally  to  the  maximum 


practicable   use   of   higher  type  facilities.   This  is  be- 
cause of  the  greater  time  savings,  comfort,  convenience,  and 
safety  inherent  in  facilities  providing  reasonably  direct 
routings  and  continuity  of  high  standards  over  a  relatively 
great  distance  (3) .   The  channeling  of  traffic  to  roads  of 
progressively  higher  standards  and  the  willingness  of  the 
public  to  do  so  explains  why  some  roads  are  more  important 
than  others  and  why  classification  is  needed. 

2.  Stability.   If  a  firm  basis  for  planning  and  fin- 
ancing a  highway  system  is  to  be  available,  it  is  necessary 
that  the  governmental  unit  have  some  assurance  that  changes 
in  the  system  v;ill  be  minimal  and  will  be  only  those  neces- 
sary to  meet  the  changing  travel  and  land  use  patterns  of 
the  state. 

3.  Economy.   Those  roads  serving  the  most  trip  pur- 
poses for  the  most  people  should  be  among  those  considered 
for  development  to  higher  standards  while  those  roads  on 
which  lov/er  standards  v;ill  suffice  may  also  be  identified. 
An  important  point  here  is  that  traffic  volume  indicates  the 
standards  of  construction  but  is  not  the  only  indicator  of 
functional  importance  (16). 

4.  Consistency.  This  implies  that  specific  methods 
and  procedures  for  classification  should  be  uniformly  ap- 
plied throughout  the  state. 

5.  Efficiency .   Efficient,  effective  management  of 
the  highway  program  by  all  involved  jurisdictions  must  be 


provided.   The  concept  of  jurisdictional  responsibility  is 
based  on  two  principles  (3): 

a)  The  geographic  limitations  of  a  government  unit 
must  be  recognized  in  measuring  its  ability  to 
plan  and  administer  effectively. 

b)  The  response  to  public  interest  must  be  recog- 
nized.  A  problem  of  statewide  interest  generates 
action  at  the  state  level  but  a  local  problem 
generates  only  local  action.   An  improperly 
classified  road  will  thus  generate  agitation 

for  action  at  the  affected  government  level  but 
no  action  at  the  jurisdictional  level. 

Advantages  and  Benefits 

Most  of  the  advantages  and  benefits  of  highway  classi- 
fication become  reality  only  if  the  classification  is 
properly  done. 

The  principal  advantages  of  classification  are: 

1.  It  provides  for  more  efficient  management  with  re- 
gard to  administrative  and  financial  responsibilities  (2,  4, 
5,  6,  7,  10,  11,  13,  14,  15,  16,  19,  25). 

2.  It  facilitates  the  fiscal  planning  of  long  range 
highway  development  and  improvement  programs  (2,  4,  5,  6,  7, 
9,  11,  13,  14,  16,  25,  26). 

3.  It  facilitates  the  planning  of  long  range  highway 
programs  (2,  4,  5,  6,  7,  11,  14,  15,  26,  26). 


4.  It  minimizes  conflicts  between  government  units 
since  each  mile  of  road  is  assigned  to  a  specific  jurisdic- 
tion (2,  5,  11,  16,  25)  . 

5.  It  establishes  a  basis  for  the  programming  of 
improvements  by  priority  or  relative  needs  (2,  5,  7,  11,  14, 
16,  26). 

6.  It  provides  a  basis  for  the  accumulation  of  cost 
data  by  v;hich  costs  of  the  highway  system  may  be  charged  to 
the  appropriate  beneficiaries,  thus  providing  for  more 
equitable  financing  (2,  6,  15,  16,  19). 

7.  It  enables  the  development  of  reasonable  minimum 
and  design  standards  for  each  system  (6,  9,  11,  14). 

8.  It  enables  the  establishment  of  an  efficient, 
specialized  staff  because  all  the  roads  of  a  particular 
system  require  the  same  degree  of  technical  conipetence  and 
ability  in  their  design,  maintenance,  and  operation  (5,  9, 
11,  15). 

9.  It  provides  for  logical,  integrated  highway  systems 
thus  promoting  the  efficiency  of  the  entire  highway  plant 
(5,  6,  11,  15,  25). 

10.  It  enables  the  application  and  comparison  of  cer- 
tain techniques,  experience,  and  results  both  within  systems 
and  betv;een  the  systems  of  various  states  (2,  16)  . 

11.  It  provides  for  an  easily  understood  picture  of 
road  needs  and  problems  by  the  general  public  and  by  highway 
administrators  (2,  11,  16). 


12.  It  permits  subclassif ication  of  any  system  for 
any  desired  purpose  (15), 

13.  It  provides  a  comparison  datum  by  which  road 
improvement  programs  and  progress  in  providing  an  adequate 
system  can  be  measured  (2,  16) . 

14.  It  provides  for  the  logical  recording  of  highway 
traffic  (16) . 

15.  It  facilitates  the  production  of  uniform  motor 
vehicle  accident  statistics  (2). 

16.  It  enables  a  more  complete  and  efficient  mainten- 
ance program  to  be  established  (2,  11). 

Classification  Criteria 

Some  of  the  references  studied  listed  a  set  of  criteria 
for  use  in  the  selection  of  highway  systems  while  in  other 
references,  the  criteria  are  implied  in  the  framework  of  a 
generalized  procedure  statem.ent.   Some  studies  provide  cri- 
teria for  the  selection  of  state,  county,  and  municipal 
systems  and  even  subclassif ications  of  each  of  these. 

However,  the  same  criteria  may  be  applied  to  the  selec- 
tion of  any  system  though  not  on  the  same  scale.   For 
instance,  one  often  used  criterion  is  the  provision  for  ser- 
vice to  points  of  traffic  interest.   This  may  be  a  city  or 
major  recreational  area  when  considering  the  state  system 
but  a  grain  elevator  for  the  county  system  and  a  shopping 
center  for  the  municipal  system. 

Some  of  the  criteria  employed  are: 


1.  Inter-center  service.   Population  centers  or  urban 
places  generate  traffic  in  different  ways  and  degrees.   The 
importance  of  a  place  determines  to  a  large  extent  the 
importance  of  roads  connecting  it  to  other  places.   Many 
classification  procedures  employ  this  concept  in  the  actual 
selection  of  highway  systems  (8,  9,  10,  11,  12,  13,  16,  17, 
19,  27).   Some  studies  refer  to  this  criterion  as  the  con- 
nection of  points  of  traffic  interest  (4,  7,  18). 

2.  Rural  access  or  balanced  area  service.   This  is  to 
ensure  that  all  areas  of  the  jurisdiction  are  reasonably 
served  by  highways.   This  is  often  done  by  trying  to  provide 
a  particular  class  of  highway,  say  a  state  highway,  within 

a  certain  distance,  say  5  or  6  miles,  for  most  of  the  people 
living  in  a  particular  area.   This  criterion  also  includes 
the  spacing  of  highways  so  that  they  are  closer  together  in 
areas  of  high  population  density  than  in  areas  of  low  popu- 
lation density  (5,  7,  9,  10,  15,  16,  17,  18). 

3.  Integrated,  continuous  systems.   The  most  important 
highway  system,  i.e.,  the  state  system,  must  be  made  up  of 
continuous  highways  such  that  an  integrated,  interconnected 
netv/ork  covering  the  entire  state  is  formed.   Lesser  systems, 
such  as  the  county  systems,  need  not  be  continuous  but  m.ust 
be  integrated  and  interconnected  with  the  state  system  in 
order  to  fulfill  their  function  as  feeder  roads  and  to  pro- 
vide an  efficient  traffic  circulation  system  (5,  7,  10,  13, 
16,  18,  27) . 
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4.  Traffic  considerations.   Occasionally,  v/hen  a 
specific  major  system  is  being  selected,  some  highways  which 
cannot  be  considered  for  inclusion  by  other  criteria  have 
high  traffic  volumes  or  perhaps  a  considerable  proportion  of 
non-local  traffic.   Highways  of  this  type  should  be  consid- 
ered for  inclusion  in  the  system  (4,  5,  7,  9,  13,  16,  17, 
18,  19). 

5.  Special  road  uses.   Roads  serving  important  interests 
should  be  considered  for  inclusion  in  a  major  system.   Exam- 
ples of  important  interests  or  uses  are  permanent  mines, 
railheads,  recreational  areas,  snail  tovms,  airports,  iso- 
lated industries,  national  defense  installations,  and  mili- 
tary bases  (5,  7,  9,  13,  17,  18,  19). 

6.  Utilization  of  existing  systems.   This  is  important 
with  respect  to  the  selection  of  one  of  two  alternate  but 
similar  routes.   The  route  on  a  presently  designated  system 
may  be  developed  and  maintained  to  higher  standards  and  should 
be  selected  (5,  7,  16,  17,  18). 

7.  Other  criteria  which  have  been  mentioned  in  some 
studies  include  topography  (5,  17),  major  post  and  school 
bus  routes  (4,  9,  16),  service  for  the  maximum  number  of 
trip  routings  with  the  least  mileage  (7,  18),  adherence  to 
legislatively  imposed  mileage  limit  controls  (7,  18),  pro- 
vision of  service  to  all  counties  and  to  all  county  seats 
(8,  10),  and  establishment  of  connections  with  major  routes 
in  adjacent  states  (9,  10,  13). 
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Elements  of  Classification 
One  of  the  more  important  questions  to  be  answeted  is 
the  number  of  systems  needed.   There  are  two  basic  considera- 
tions involved  in  answering  this  question: 

1.  Character  of  traffic  service  each  system  must 
provide . 

2.  Government  units  available  to  manage  each  system. 
For  classifications  purposes,  highway  traffic  consists 

of  trips  of  three  types:   those  of  statewide  interest;  those 
of  local  rural  interest;  and  those  of  local  urban  interest. 
Further  analysis  indicates  that  public  highways  must  satisfy 
trip  demands  falling  into  three  functional  categories. 

1.  Longer,  relatively  high  speed  trips  with  only  in- 
cidental land  service. 

2.  Shorter,  lower  speed  trips  within  a  limited  area 
v;ith  simultaneous  land  service. 

3.  Short,  low  speed  trips  where  the  emphasis  is  on 
land  service. 

These  considerations  give  rise  to  three  general  systems 
in  both  rural  and  urban  areas  as  shown  below  (18): 
Rural  Urban 

State  Highways  State  Highways 

County  Arterial  Highways  Arterial  Streets 

Local  Roads  Local  Streets 

The  analogy  of  a  highway  system  to  house  construction 
has  been  made  in  a  South  Dakota  study  (18) .  In  building  a 
house,  the  needs  are  a  few  main  beams,  some  stringers,  and 
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many  floor  boards.   Each  element  serves  its  purpose  in 
distributing  the  imposed  loads  with  the  main  beams  having 
the  greatest  capacity  to  carry  loads  and  the  floor  boards 
having  the  least.   A  house  built  of  all  main  beams  could  do 
the  job  but  at  very  high  cost.   All  floor  boards  could  not 
do  the  job. 

A  highway  system  is  needed  such  that  elements  of  in- 
creasing capacity  to  handle  traffic  are  available  to  dis- 
tribute traffic  and  to  do  so  economically.   Thus,  a  few 
main  roads  (state  highv;ays),  a  collector  system  (arterial 
highways  or  streets),  and  a  general  coverage  system  (local 
roads  or  streets)  are  needed. 

Government  units  available  to  manage  highway  systems 
include  the  Federal,  state,  county,  township,  and  municipal 
branches . 

Federal  classification  efforts  have  been  confined  to 
statutes  v;hich  have  directed  lesser  government  units  to 
classify  certain  highway  systems  in  order  to  be  eligible  for 
Federal-Aid.   These  efforts  have  resulted  in  the  establish- 
ment of  a  nationwide  system  of  Interstate  highways,  state 
primary  and  secondary  systems,  urban  extensions,  and  even 
county  arterial  systems  of  roads  eligible  for  Federal-Aid. 
However,  none  of  these  systems  are  managed  by  the  Federal 
government . 

The  classification  has  been  made  by  the  governmental 
units  involved  but  with  the  cooperation  and  approval  of  the 
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Federal  government  through  the  agency  of  the  Bureau  of  Pub- 
lic Roads. 

Federal  efforts  have  resulted  in  the  following  benefits 
(14)  : 

1.  An  integrated  network  of  Interstate  routes  serving 
all  states  and  a  major  intrastate  system. 

2.  Greater  uniformity  in  administrative  practices. 

3.  Improvement  of  all  major  systems  through  alloca- 
tion of  funds . 

4.  Adequate  standards  of  improvement  and  construction 
for  each  system. 

5.  Protection  of  the  highway  investment  by  requiring 
all  participating  jurisdictions  to  maintain  all  system 
mileage . 

The  individual  states,  however,  have  been  the  major 
unit  of  government  to  classify  highway  mileages.   All  states 
have  selected  a  primary  system  under  the  terms  of  the  1921 
Federal  Aid  Highway  Act  which  limited  the  system  to  7  per 
cent  of  total  state  mileage,  among  other  provisions. 

Some  states  also  have  designated  state  secondary  systems 
and  special  purpose  systems.   Many  states  have  other  system 
subclassif ications  for  administrative  or  other  reasons. 

The  state  of  VJashington,  for  instance,  has  subclassi- 
fied  the  state  highway  system  into  five  groups  (22): 

1.  Interstate  system. 

2.  Principal  state  systems.   Includes  highways  con- 
necting cities  of  20,000  or  more. 
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3.  Major  state  system.   Includes  roads  providing 
connections  for  cities  of  over  1000  or  which  serve  major 
tourist,  commercial,  or  industrial  traffic  movements. 

4.  Collector  state  system.   Includes  roads  serving 
other  populated  areas. 

5.  Other  state  system. 

The  principal  advantage  of  subclassif ication  is  that 
it  provides  a  sound  basis  for  the  determination  of  improve- 
ment and  construction  priorities. 

Need  for  Reclassification 
The  failure  to  properly  classify  roads  underlies  the 
basic  highway  problem  confronting  many  states  according  to 
the  Automotive  Safety  Foundation  (1).   Allocations  of  funds 
are  sometimes  based  on  political  or  financial  expediency 
rather  than  on  engineering  techniques  which  measure  the 
service  a  road  provides. 

Improper  classification  results  in  many  inequities. 
T^ong  these  are:   overtaxing  the  facilities,  manpower,  and 
funds  of  a  government  unit  to  provide  satisfactorily  for  each 
mile  of  highway  under  its  control;  disproportionate  distri- 
bution of  funds;  inequitable  distribution  of  costs  to  the 
highway  users;  and  misplaced  responsibility  by  charging  road 
officials  to  provide  for  services  in  which  they  have  little 
or  no  interest. 

Even  if  an  entire  road  system  were  properly  classified 
at  one  time,  there  are  several  reasons  why  this  will  not 
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always  be  the  case:   changes  in  land  use  and  population  re- 
sulting in  different  patterns  of  travel;  growth  in  popula- 
tion, vehicle  ownership,  and  vehicle  usage  resulting  in 
increased  traffic  and  greater  demand  for  more  and  better 
highways;  the  fast  approaching  completion  of  the  Interstate 
system  which  is  changing  traffic  patterns  as  motorists 
choose  better,  more  direct  routes;  changes  in  economy  and 
the  development  of  industry  and  natural  resources;  and  the 
construction  of  new  recreational  facilities  such  as  parks 
and  reservoirs. 

Thus,  a  periodic  review  of  highway  system  classifica- 
tion is  needed  to  resolve  existing  inequities  and  to  re- 
establish the  benefits  and  advantages  of  proper  highway 
classification. 

Improper  classification  of  a  road,  in  either  too  high 
or  too  low  a  class,  prevents  the  road  from  functioning 
efficiently.   A  local  road  may  be  quite  important  to  the 
community  it  serves.   It  gets  priority  with  respect  to  main- 
tenance, repair,  and  snow  removal.   If  this  road  were 
transferred  to  state  control,  its  importance  relative  to 
other  state  roads  places  it  near  the  end  of  the  list  for 
essential  services.   The  result  is  disruption  of  traffic 
and  land  service  and  bitterness  on  the  part  of  the  community, 
The  converse  is  equally  true.   A  truly  state  road  transfer- 
red to  county  control  will  not  receive  the  attention  it 
should.   As  a  result,  the  motoring  public  does  not  get  the 
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transportation  service  it  wants  or  needs  and  for  which  it 
has  paid. 

It  is  recognized  that  the  completion  of  the  Interstate 
highway  system  will  have  a  significant  effect  on  existing 
travel  patterns  not  only  in  Indiana  but  throughout  the 
Nation.   Indeed,  it  has  been  estimated  that  the  completed 
Interstate  system,  comprising  slightly  more  than  one  per 
cent  of  the  Nation's  total  road  and  street  mileage,  will 
carry  21  per  cent  of  all  traffic  (29). 

Hov/ever,  study  of  the  proposed  Interstate  system  shows 
that  in  many  instances,  Interstate  routes  closely  parallel 
or  coincide  with  existing  major  highways,  such  as  1-70  and 
US-40  across  Indiana.   xluch  of  the  traffic  on  the  existing 
route  will  simply  shift  over  to  the  Interstate  route. 

Other  cases  of  traffic  diversion  as  well  as  traffic 
generation  are  expected  because  drivers  will  attempt  to 
utilize  the  Interstate  routes  (and  others  built  to  similar 
standards)  to  the  greatest  extent  possible.   The  reasons  for 
this  are  many  (29): 

1.  Cost  reduction.   Studies  have  shown  in  souce  cases 
that  the  cost  of  a  freeway  is  balanced  out  by  motorist's 
savings  in  less  than  10  years. 

2.  Safety.   Use  of  the  Interstate  System  will  save  at 
least  5000  lives  a  year. 

3.  Freedom  in  travel.   Long  range  highway  travel  for 
any  purpose  by  any  mode  will  be  faster  with  greater  comfort 
and  less  strain. 
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4.   Reduction  of  urban  congestion.   Routes  through  and 
around  larger  cities  and  those  which  bypass  the  smaller 
cities  will  provide  relief  from  the  traffic  congestion  which 
today  clogs  many  urban  streets. 

The  concept  of  traffic  diversion  and  generation  is  not 

new.   It  was  pointed  out  in  the  historic  "Road  and  Canal 

Report"  by  Albert  Gallatin  in  1808  when  he  wrote  of  the 

effects  of  the  construction  or  improvement  of  routes  thusly 

(31)  : 

...The  general  gain  is  not  confined  to  the  difference 
between  the  expense  of  the  transportation  of  those 
articles  which  had  been  formerly  conveyed  by  that 
route,  but  many  which  were  brought  to  market  by 
other  channels  will  then  find  a  new  and  more  advan- 
tageous direction;  and  those  which  on  account  of 
their  distance  or  weight  could  not  be  transported 
in  any  manner  whatsoever,  will  acquire  a  value,  and 
become  a  clear  addition  to  the  national  wealth. 

It  is  often  pointed  out  that  there  have  been  large  in- 
creases in  population,  vehicles,  and  vehicle  miles  of  travel 
in  this  country  over  the  past  several  decades  and  that  this 
growth  is  expected  to  continue.   During  the  fifteen  year 
period  from  1948  to  1963,  the  number  of  vehicles  and  vehicle 
miles  almost  doubled  while  population  and  traffic-related 
fatalities  each  increased  by  about  one  third  (32). 

It  is  significant,  perhaps,  to  note  that  road  mileage 
during  this  same  interval  increased  by  a  factor  of  about 
1.04.   This  is  in  keeping  with  what  probably  is  National 
policy  in  respect  to  highway  transportation  facilities  as 
suggested  by  the  following  (29): 
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...most  construction  of  'nev/'  roads  actually  is  the 
replacement  of  betterment  of  existing  facilities. 
A  highway-improvement  program  therefore  is  not  designed 
to  achieve  'more'  highways  so  much  as  it  is  to  achieve 
'better'  or  'more  adequate'  ones. 


Purpose  of  this  Research 
Two  important  points  emerge  from  the  discussion  in  the 
previous  section: 

1.  Traffic  patterns  are  expected  to  change  so  that 
som.e  of  the  important  state  highways  of  today  will  not  be 
too  important  in  the  future  and  vice  versa. 

2.  Highway  improvements,  geared  to  provide  for  traffic 
requirements,  should  be  concentrated  on  the  existing  roads 
which  v;ill  be  important  state  highways  in  the  future. 

The  purpose  of  this  research  was  to  develop  a  method 
of  selecting  the  important  highways  of  the  future  in  Indiana, 
A  statewide  study  of  intercity  travel  desires  was  selected 
as  the  means  for  achieving  this  purpose.   The  present  State 
Highway  System  of  Indiana  as  shown  in  Figure  1  was  used  in 
this  analysis  with  the  Interstate  System  assumed  to  be  com- 
pleted and  fully  operational. 

The  scope  was  limited  in  that  a  complete  reclassifica- 
tion of  the  Indiana  state  highway  system  was  not  contem- 
plated because  only  intercity  travel  of  a  non-local  nature 
was  considered.   As  the  State  highway  system  primarily 
serves  intercity  travel,  however,  only  a  small  number  of 
miles  of  State  highway  serving  points  of  major  interest 
(such  as  reservoirs  and  parks)  was  not  included  in  the  study, 
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FIGURE    1.      THE     STATE      HIGHWAY      SYSTEM      OF      INDIANA 
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The  principal  aims  of  this  study  were  to  subclassify 
those  highways  of  primarily  statewide  interest  in  Indiana 
and  to  pinpoint  those  which  should  be  considered  for  recon- 
struction to  freeway  standards. 
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HIGHWAY  CLASSIFICATION  IN  INDIANA 

Historical  Development 

Road  building  in  Indiana  during  the  immediate  period 
following  statehood  in  1816  was  generally  confined  to  the 
National  Pike  or  Cumberland  Road  and  the  Michigan  Road  as 
shown  in  Figure  2.   Financing  of  these  projects  was  derived 
from  a  five  per  cent  tax  imposed  on  the  sale  of  Indiana  pub- 
lic lands.   Two  per  cent  was  to  be  retained  by  the  Federal 
government  for  the  construction  of  transportation  facilities 
to  and  through  Indiana  and  the  remaining  three  per  cent  was 
to  be  used  for  internal  improvements  within  Indiana  (31). 

The  National  Pike  was  built  by  the  Federal  government 
with  the  "2  per  cent  fund"  while  the  Michigan  Road  was  built 
by  the  State  government  with  the  "3  per  cent  fund."   Neither 
of  these  roads  were  very  good,  being  largely  impassable  much 
of  the  year,  but  they  did  open  Indiana  to  large   scale  pio- 
neer settlement. 

Many  other  road  and  canal  projects  were  initiated  by 
the  State  during  this  period  but  the  shortage  of  funds  and 
piecemeal  construction  resulted  in  almost  no  improvement 
in  the  transportation  problems  of  Indiana  to  1836.   Three 
conditions  led  Indiana  to  pass  an  internal  improvements  bill 
known  as  the  System  of  1836  (33): 
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1.  The  Federal  administration  at  this  time  was  opposed 
to  participation  in  internal  improvements, 

2.  Several  eastern  states  were  successfully  carrying  out 
internal  improvement  programs  of  their  own,  and 

3.  Public  disgust  with  the  transportation  problem   had 
made  it  a  political  issue. 

The  bill  provided  that  a  loan  of  up  to  $10,000,000  be 
made  for  a  period  of  25  years  at  an  interest  rate  of  not  over 
five  per  cent.   The  effects  of  this  bill  was  a  disaster  of 
immense  proportions  to  Indiana  (33).   Average  state  revenues 
at  this  time  were  $75,000  per  year  and  the  bill  made  no  pro- 
vision other  than  tolls  for  paying  the  $500,000  yearly  int- 
erest due  on  the  loan.   Every  road  and  canal  project  auth- 
orized by  the  bill  was  begun  almost  immediately,  none  was 
ever  completed,  and  virtually  no  revenue  was  ever  collected 
by  the  State.   Work  on  the  System  stopped  in  1839  and  by  1844 
the  State  debt  was  over  $15,000,000. 

This  debt  was  settled  later  with  great  favor  to  the 
State  but  resulted  in  a  complete  loss  of  its  credit.   Thus, 
the  new  State  Constitution  of  1851  imposed  strict  limitations 
on  the  amount  of  indebtedness  that  could  be  incurred  and  pro- 
hibited the  State  from  borrowing  for  public  improvements  (33). 

This  was  not  then  an  especially  significant  factor  in 
the  decline  of  road  construction  and  maintenance  because  the 
railroads  had  captured  virtually  the  entire  market  in  the 
transport  of  passengers  and  freight.   In  fact,  the  period 
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from  1850  to  1880  was  known  as  the  "Dark  Ages  of  the  Road" 
(31)  . 

Legislation  passed  in. 1877  authorized  the  several  coun- 
ties of  Indiana  to  issue  bonds  for  constructing  and  improv- 
ing roads  (34).   The  bonds  were  to  be  retired  by  assessments 
on  the  benefited  property.   Over  the  subsequent  40  years, 
many  miles  of  gravel  roads  were  built  and  maintained  by  the 
counties  as  a  result  of  this  legislation. 

Several  groups  of  cycling  enthusiasts  formed  the  League 
of  American  Wheelmen  in  1880  and  began  a  nationwide  campaign 
for  better  roads,  thus  initiating  the  Good  Roads  Movement  in 
many  states  (31).   Their  efforts  resulted  in  the  passage  by 
several  states  of  local  road-aid  laws,  led  by  New  Jersey  in 
1891,  and  the  establishment  of  the  Office  of  Road  Inquiry  by 
the  U.  S.  Congress  in  1893.   Other  Federal  agencies,  notably 
the  Post  Office  and  the  Forest  Service,  developed  an  interest 
in  good  roads  as  well. 

These  developments  coupled  v/ith  the  overwhelming  public 
acceptance  of  the  motor  vehicle  led  to  the  passage  of  the 
Federal-Aid  Road  Act  of  1916  which  constituted  the  basic 
Federal  highway  law  until  the  codification  in  1958  of  all 
highway  statutes  as  Title  23  of  the  United  States  Code  (35). 

The  Modern  Era 
The  1916  Federal-Aid  Road  Act  was  the  first  comprehen- 
sive step  toward  the  development  of  a  nationwide  system  of 
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interstate  highways.   Among  the  provisions  of  the  Act  were 
included  (31,  36): 

1.  Each  State  legislature  must  specifically  assent  to 
the  terms  of  the  Act  and  designate  a  State  agency  (highway 
department)  to  cooperate  with  the  Federal  government. 

2.  Projects  contemplated  by  the  State  agency  were  to 
be  approved  and  inspected  by  a  designated  Federal  agency. 

3.  Each  State  must  match  the  Federal  appropriation 
with  funds  of  its  own. 

4.  The  fund  apportionment  to  each  State  was  based  on 
its  area,  population,  and  miles  of  rural  mail  delivery  routes, 

5.  Federal-Aid  was  limited  to  construction  of  a  designa- 
ted  primary  highway  system  which  was  later  limited  to  7  per 
cent  of  the  existing  miles  of  road  within  each  State  by  the 
Federal-Aid  Highway  Act  of  1921. 

As  a  result  of  this  Act,  the  Indiana  State  Highway 
Commission  was  formed  by  the  State  legislature  in  1917  (37). 
Because  this  legislation  was  immediately  attacked  in  the 
courts,  it  was  repealed  and  replaced  by  the  1919  Highv;ay  Act 
of  Indiana  ( 34 ) . 

The  principal  accomplishment  under  the  terms  of  the 
1917  Act  was  the  designation  of  almost  990  miles  of  "main 
market  highways"  which  connected  the  main  market  centers  of 
Indiana.   These  roads  are  shown  in  Figure  3. 

Under  the  terms  of  the  1919  Act,  the  State  Highv;ay 
Commission  was  charged  with  the  laying  out  of  a  connected 
system  of  State  highways  which  would  reach  every  county  seat 


26 


Ml       C       H      I      G      A      N 


FIGURE    3.      MAIN       MARKET      HIGHWAYS    "     1917 
(  SOURCE  :     REFERENCE      37  ) 
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and  each  city  of  5000  or  more  inhabitants.   Thus,  a  3200  mile 
State  highway  system,  as  shown  in  Figure  4,  was  approved  and 
accepted  in  April,  1920  by  the  Commission  for  construction 
and  maintenance. 

This  system  amounted  to  about  5  per  cent  of  Indiana's 
road  miles  with  provision  for  an  increase  to  about  10  per 
cent  as  funds  became  available  to  properly  maintain  the 
additional  miles  of  road  (38). 

Throughout  the  succeeding  years,  the  bulk  of  State  high- 
way legislation  has  had  to  do  with  the  means  of  raising  rev- 
enues and  its  subsequent  apportionment  to  the  State  Highway 
Commission  and  to  the  counties  and  cities  of  the  State. 

The  Federal  government,  however,  passed  legislation  in 
1944  authorizing  Federal  Aid  for  a  system  of  secondary  high- 
ways as  well  as  urban  extensions  of  primary  and  secondary 
highvv'ays .   This  legislation  also  provided  for  the  selection 
of  a  national  system  of  interstate  highways  limited  to  40,000 
miles  (36).   Other  Federal  and  State  legislation  has  also 
resulted  in  various  subsystems  of  park,  forest,  and  insti- 
tutional roads.   However,  the  bulk  of  Federal  highway  legis- 
lation, just  as  for  Indiana,  has  been  devoted  to  financing 
considerations . 

Over  the  years,  the  Indiana  State  Highway  Commission 
and  the  system  of  roads  for  which  it  is  responsible  have 
enjoyed  substantial  growth.   The  total  length  of  the  State 
highway  system,  which  was  limited  to  a  maximum  of  12,000 


28 


M    I     C 


K         E  ^ 


FIGURE     4.      INDIANA      STATE      HIGHWAY      SYSTEM     -     1920 
(SOURCE:     REFERENCE      38) 
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miles  by  a  1937  Indiana  Statute  (34),  reached  approximately 
11,500  miles  in  1966  (39).   About  4,170  miles  of  this  total 
are  on  the  Federal-Aid  Primary  System  and  5,800  on  the 
Federal-Aid  Secondary  System.   In  addition,  there  are  12,970 
miles  of  Federal-Aid  Secondary  roads  under  county  control. 
A  total  of  1,115  miles  of  Interstate  highway  are  scheduled 
for  Indiana  with  more  than  half  open  to  traffic  as  of  1966. 

Some  idea  of  the  growth  of  the  highway  system  and  the 
reasons  for  its  development  can  be  obtained  by  comparing  some 
statistics  for  1920  and  1965.   The  miles  of  road  increased 
from  72,000  to  91,000,  the  State  highway  system  mileage  from 
3,200  to  11,160,  motor  vehicle  registrations  from  0.33  million 
to  2.4  million  ,  population  from  2.9  million  to  4.9  million, 
and  receipts  of  the  State  Highway  Commission  from  4.14  mil- 
lion to  $176.8  million  with  Federal-Aid  from  $0.49  million 
to  $88.6  million  (38,40). 

The  impact  of  Federal-Aid  participation  of  90  percent 
of  the  cost  of  the  Interstate  system  can  be  readily  visual- 
ized from  two  facts.   Total  Federal-Aid  to  Indiana  in  the  13 
year  span  from  1941  to  1953  amounted  to  $69.3  million  while 
Federal-Aid  in  1965  alone  amounted  to  $88.6  million  of  which 
$66.3  million  was  allocated  for  the  Interstate  system.   Up 
until  about  1960,  Federal- Aid  amounted  to  about  20  per  cent 
of  the  Commission's  receipts  but  since  then,  it  has  amounted 
to  about  50  per  cent  (40) . 
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Classification  in  Indiana 

The  Indiana  State  Highway  Commission  Law  of  1919  charged 
the  newly  created  Commission  to  lay  out  a  system  of  State 
highways  to  reach  each  and  every  county  seat  and  each  and 
every  city  of  over  50  00  population  and  to  provide  connections 
to  the  main  trunk  highways  of  adjacent  states  such  that  a 
continuous  system  of  improved  highways  would  be  formed.   The 
system  was  to  be  constructed,  reconstructed,  repaired,  and 
maintained  by  the  Commission  out  of  State  highway  funds  (38). 

This  same  Act  also  contained  the  specific  assent  of 
Indiana  to  the  provisions  of  the  1916  Federal-Aid  Road  Act 
and  all  subsequent  acts  amendatory  and  supplementory  thereto. 
The  State  Highway  Commission  was  authorized  to  cooperate 
with  the  Federal  government  under  any  Federal  law  in  any 
manner  necessary  to  secure  for  Indiana  the  proportion  of  any 
Federal  appropriation  which  may  be  made  in  the  future  (35). 

The  Indiana  legislature  has  seen  fit  to  give  the  Indiana 
State  Highway  Commission  much  of  the  authority  with  regard 
to  highv;ays  but  generally  has  not  provided  guidelines  con- 
cerning legislative  intent  (25).   For  instance,  the  Commission 
has  the  authority  to  designate  a  State  primary  system  but 
Indiana  law  includes  few  definitions,  standards,  or  criteria 
for  designating  the  system. 

The  Commission  has  the  following  authority  (25): 

1.   It  may  add  routes  to  the  system  if  funds  are  avail- 
able and  if,  in  the  opinion  of  the  Commission,  the  route  would 
make  a  desirable  addition  to  the  State  highway  system. 
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2.  It  may  take  over  roads  from  other  systems. 

3.  It  may  subclassify  the  State  highway  system  into 

two  or  more  classes  and  make  any  classification  changes  deemed 
necessary. 

4.  It  may  relocate  existing  highways  to  promote  public 
safety  and  convenience  or  if  the  relocation  will  provide  a 
more  serviceable  highway  system  than  is  possible  under  the 
existing  location. 

5.  It  may  delete  or  abandon  State  highways  as  necessary. 

6.  It  may  designate  municipal  connecting  links  and  add 
to,  relocate,  or  delete  them. 

In  some  cases,  the  legislature  has  imposed  some  condi- 
tions on  the  Commission's  authority  but  even  here,  the  con- 
ditions do  not  hold  if  the  Commission  feels  that  changes  are 
necessary  and  beneficial,  i.e.,  in  the  public  interest.   The 
legal  basis  for  much  of  the  Commission's  authority  stems 
from  the  fact  that  a  State  highway  is  one  so  designated  by 
the  Commission  (25). 

Indiana  statutes  currently  provide  for  only  two  rural 
highway  systems:   the  State  system  and  the  several  county 
systems.   Each  of  the  several  counties  has  the  authority  to 
make  changes  in  its  s/stem.   County  roads  include  all  rural 
roads  not  designated  in  the  State  system.   Generally,  when 
any  highway  is  deleted  from  the  State  system,  it  reverts  to 
the  counties  but  in  some  instances  may  revert  to  the  orig- 
inal landowner  (25). 
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CLASSIFICATION  STUDIES  IN  OTHER  STATES 

Pennsylvania 
Highway  classification  in  Pennsylvania  was  accomplished 
on  a  county  by  county  basis  (20) .   The  road  system  of  a  county 
was  examined  and  if  a  road  appeared  to  meet  the  definition 
of  one  of  four  proposed  systems,  it  was  so  classified.   Maps 
of  adjoining  counties  were  later  compared  and  inconsistenc- 
ies at  county  boundaries  eliminated.   The  classifications 
were  field  checked  and  public  hearings  held  in  each  county 
prior  to  legislative  adoption  of  the  systems. 

Illinois 

The  method  of  classification  proposed  in  Illinois  was 
based  on  the  concept  that  a  highway's  importance  may  be 
measured  by  its  predominant  usage.   The  assumption  was  then 
made  that  usage  of  a  highway  can  be  measured  by  the  economic 
importance  or  traffic  attraction  of  urban  places  connected 
by  the  highway  and  that  highways  connecting  urban  places 
of  similar  economic  importance  will  have  similar  functional 
usage . 

According  to  a  Louisiana  study  (12),  there  are  four 
basic  principles  which  explain  the  ability  of  some  urban 
places  to  attract  traffic  to  a  greater  or  lesser  degree  than 
other  urban  places.   These  principles  are: 
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1.  People  make  trips  so  that  the  larger  the  place,  the 
greater  the  number  of  trips  attracted  or  produced. 

2.  People  make  trips  to  see  other  people  or  to  buy 
and/or  sell.   Therefore,  urban  places  with  more  people, 
stores,  markets,  and  other  facilities  will  attract  more 
trips  than  those  with  less  such  resources. 

3.  The  larger  an  urban  place,  the  farther  people  will 
come.   This  is  because  the  greater  population  has  more  widely 
scattered  personal  relationships  and  also  has  greater  and 
more  specialized  trade,  service,  cultural,  and  recreational 
resources,  many  of  which  are  not  available  at  home. 

4.  Trips  make  traffic  such  that  highways  connecting 
large  trip  producing  and  attracting  urban  places  will  gen- 
erally carry  the  heaviest  traffic. 

In  Illinois,  the  economic  importance  of  an  urban  place 
was  measured  by  determining  its  economic  rating  (15).   Of 
the  many  factors  characteristtic  of  an  urban  place,  the  fol- 
lowing four  were  selected  to  establish  the  economic  rating: 
population  of  the  immediate  trade  area;  bank  resources  of 
the  immediate  trade  area;  newspaper  circulation  of  those 
published  in  the  immediate  trade  area;  and  the  assessed  val- 
uation of  the  trade  center.   A  Michigan  study,  performed 
independently  but  at  the  same  time,  included  the  additional 
factor  of  retail  sales  in  the  trade  center  (28) . 

The  rating  was  determined  as  follows.   A  factor,  say 
newspaper  circulation,  was  summed  over  every  urban  place  in- 
cluded in  the  study  to  obtain  a  state  total.   The  newspaper 
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circulation  of  a  particular  city  was  divided  by  the  state 
total  to  obtain  a  ratio.   The  ratios  resulting  from  the  con- 
sideration of  all  four  factors  were  added  and  divided  by  four 
yielding  the  economic  rating  of  each  urban  place. 

Factor  values  for  Chicago  were  not  included  so  as  not  to 
distort  the  rating  scheme  and  each  factor  was  assumed  to  be 
of  equal  value  in  computing  the  rating. 

The  econoiric  ratings  wei"e  plotted  on  a  uniform  scale  in 
decreasing  order.   Six  apparent  groupings  v:ere  readily  ob- 
served and  urban  places  falling  into  each  group  were  called 
metropolitan  centers;  regional  centers;  major  market  centers; 
market  centers  A;  market  centers  B;  and  minor  market  centers. 

Highways  were  then  selected  to  connect  those  places  of 
highest  economic  importance  with  each  other  and  with  simi- 
lar places  outside  the  State.   Highways  were  next  selected 
to  connect  places  second  in  importance  with  each  other  and 
with  those  of  greater  importance.   This  process  continued 
through  the  6  place  groups  such  that  a  basic  highv/ay  system 
was  formed.   This  method  of  ranking  cities  and  selecting 
a  basic  highway  system  will  be  referred  to  later  as  the 
Illinois  method. 

The  State  primary  system  was  composed  of  the  Interstate 
system  and  all  highways  selected  to  interconnect  cities 
through  those  identified  as  market  center  B.   However,  sev- 
eral other  factors  were  employed  to  assure  a  balanced, 
integrated  system.   Some  of  these  were:   high  volume  roads 
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performing  a  statewide  service;  roads  connecting  county  seats 
not  already  connected;  roads  needed  to  provide  a  primary 
highway  within  10  miles  distance;  roads  needed  to  serve 
recreational  facilities;  and  important  roads  in  the  Chicago 
area . 

South  Dakota 

In  a  South  Dakota  study,  cities  were  divided  into  8 
groups  on  the  basis  of  population,  retail  sales,  postal  re- 
ceipts, and  livestock  trading  (18).   Each  factor  was  measured 
for  the  city  and  not  the  immediate  trade  area. 

Aside  from  this,  the  Illinois  method  was  used  to  rank 
the  cities  and  to  select  a  basic  highway  system.   Supple- 
mentary routes  to  assure  balanced,  integrated  service  were 
selected  by  the  following  criteria:   provision  for  system 
continuity  and  interconnection;  accommodation  of  large  traf- 
fic flows;  provision  for  balanced  area  service;  accommodation 
of  special  road  uses;  and  the  observance  of  legislatively  im- 
posed mileage  limit  controls. 

Kentucky 
A  1962  Kentucky  study  reported  that  all  urban  places 
of  over  1000  inhabitants  were  ranked  by  population  (19). 
However,  all  those  of  over  2500  were  ranked  by  the  Illinois 
method  using  the  factors  of  population,  manufacturing  employ- 
ment, retail  sales,  and  bank  assets  for  each  individual 
city.   Urban  places  with  over  2500  persons  were  divided  into 
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7  groups  while  those  with  populations  between  1000  and  2500 
were  divided  into  5  groups.   Nearby  cities  in  adjacent  states 
were  similarly  grouped. 

The  cities  were  plotted  one  at  a  time  on  a  map  starting 
with  the  most  important  one.   As  each  city  was  plotted,  the 
most  direct,  feasible  route  to  connect  it  to  those  already 
plotted  was  selected.   All  roads  needed  to  connect  group  I 
cities  were  called  group  I  arterials  and  so  on  through  all 
12  groups.   A  thirteenth  group  of  roads  was  selected  to  ac- 
count for  routes  carrying  larger  than  average  traffic  and 
routes  connecting  large  traffic  generators  such  as  recre- 
ational and  military  areas. 

Manitoba 

All  cities  of  over  720  population  in  Manitoba  were 
divided  into  6  groups  on  the  basis  of  population  only  (16). 
A  basic  highway  system  was  then  selected  using  the  Illinois 
method. 

Supplementary  routes  were  selected  to  provide  for  rural 
access  service  and  an  integrated  network. 

Connecticut 

All  towns  in  Connecticut  were  divided  into  four  groups 
by  the  Illinois  method  using  the  following  factors:   popula- 
tion, employment,  grand  levy,  motor  vehicle  registration, 
and  sales  and  use  tax  receipts  (17). 

Roads  were  then  selected  which  served  the  predominant 
flow  of  traffic  between  any  two  towns.   The  roads  needed  to 
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connect  cities  in  groups  I  and  ii  became  the  State  primary 
system.   The  roads  needed  to  connect  the  remaining  cities 
became  the  State  secondary  system.   A  low  mileage  group  of 
State  special  service  roads  was  selected  to  provide  access 
to  recreational  areas,  institutions,  and  military  installa- 
tions . 

The  Connecticut  study  was  somewhat  unique  in  that  it 
included  a  study  of  existing  laws  and  developed  legislation 
to  implement  the  proposed  system.   The  legislation  became 
law  and  the  process  of  road  transfer  from  one  jurisdiction 
to  another  was  outlined  and  may  be  completed  at  this  time. 
A  special  arbitration  procedure  for  resolving  differences 
was  also  established. 

AASHO-NACO  Guide 

Part  III  of  the  AASHO-NACO  Guide  is  a  description  of 
the  technical  procedure  for  classifying  rural  highways  and 
supposedly  indicated  the  extent  to  v;hich  a  formalized  pro- 
cedure was  then  available  (3).   The  procedure  is  divided  into 
two  main  parts  and  a  brief  summary  of  each  follows: 

In  laying  out  an  arterial  network,  the  following  steps 
are  required: 

1.   Rank  population  centers  on  a  basis  which  indicates 
their  relative  radius  of  traffic  attraction.   Population 
alone  is  considered  a  sufficient  ranking  factor  due  gener- 
ally to  the  incomplete  availability  of  other  information, 
especially  for  the  smaller  urban  places. 
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2.  Plot  the  centers  graphically,  in  order  of  ranking, 
and  divide  them  into  6  to  8  similar  groups. 

3.  Repeat  steps  one  and  two  for  any  out-of-state  cen- 
ters v/hich  judgment  indicates  have  significant  traffic 
attraction . 

4.  Plot  each  group  of  centers  and  their  appropriate 
urban  limits  on  a  state  map. 

5.  Connect  the  largest  size  centers  by  the  most  direct, 
logical  routes.   Connect  the  next  largest  centers  and  con- 
tinue working  down  through  the  smaller  centers.   Good  judg- 
ment is  essential  in  determining  which  connections  should  be 
made  and  reference  to  a  traffic  flow  map  is  helpful  in 
reaching  decisions. 

5.   Log  the  routes  in  the  sequence  of  their  selection. 

7.  Determine  the  smallest  size  centers  to  be  connected 
by  noting  at  which  point  miles  of  road  are  being  added  with- 
out a  significant  increase  in  the  vehicle  miles  of  travel 
served  or  by  visually  noting  that  relatively  few  long  dis- 
tance trips  are  being  served. 

8.  Add  other  routes  as  required  including: 

a)  Service  to  other  major  traffic  generators 
(recreational  or  military) . 

b)  Significant  corridor  movements. 

c)  Service  to  all  areas  of  the  State. 

d)  Additions  needed  for  continuity. 

9.  Consider  inclusion  of  alternate  routes  where  one 
facility  cannot  handle  all  movement,  where  one  facility  is 
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a  parkway  from  which  commercial  vehicles  are  excluded,  or 
where  two  alternate  routes  are  separated  by  a  geographical 
barrier . 

In  laying  out  a  collector  network,  the  procedure  must 
be  generalized  to  an  even  greater  degree  since  information 
available  at  this  level  is  often  not  precise  and  seldom,  com- 
plete.  Generally,  it  should  be  performed  at  the  county 
level  and  information  obtained  on  the  following  factors: 
location  of  other  population  centers  not  on  the  arterial 
system;  location  of  heavier  than  average  traffic  flows;  lo- 
cation of  freeway  interchanges  and  river  crossings;  location 
of  important  traffic  generators;  and  rural  population  dens- 
ity and  land  use  distribution  throughout  the  county. 

A  collector  system  is  then  selected,  which  in  conjunc- 
tion with  the  arterial  system,  will  most  efficiently  meet 
the  highway  needs  of  the  county. 

Use  of  this  procedure  generally  results  in  an  arterial 
system  comprising  from  7  to  10  per  cent  of  rural  mileage 
while  the  extent  of  the  collector  system  ranges  from  20  to 
25  per  cent  of  rural  mileage. 

Ontario 
The  Ontario  Department  of  Highways  in  1957  developed 
two  classification  techniques  to  aid  in  defining  the  func- 
tion of  each  highway  (24).   These  techniques  were  Inter- 
center  Service  and  Rural  Access  Service. 
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Intercenter  service  is  based  on  the  concept  that  the 
importance  of  a  highway  may  be  measured  by  the  total  amount 
of  service  it  provides  to  interconnect  population  centers. 
Accordingly,  each  city  of  over  3500  population  was  placed 
in  one  of  6  groups  identified  in  decreasing  order  as  M,  A, 
B,  C,  D,  and  E. 

External  origin  -  and  -  destination  surveys  were  used 
to  determine  the  average  number  of  daily  through  trips  for 
each  type  of  city  connection.   For  example,  M  to  M  con- 
nections averaged  660  daily  trips,  M  to  A  connections  620,  and 
so  on.   The  average  numbers  of  daily  trips  were  found  to  con- 
gregate at  values  of  600,  300,  and  100.   Point  values  of  6, 
3,  and  1  were  assigned  to  highway  connections  between  each 
city  pair  making  up  one  of  the  three  characteristics  con- 
nections.  Judgment  was  used  to  decide  if  a  connection  should 
be  disregarded  when  the  distance  between  the  two  cities  was 
great . 

The  points  assigned  to  each  highway  section  were  summed 
and  this  number  used  to  classify  highways.   However,  since 
Southern  Ontario  has  88  per  cent  of  the  population  in  7  per 
cent  of  the  area,  different  point  range  classification  cri- 
teria were  used  in  Southern  and  Northern  Ontario.   These 
are  shown  in  the  following  table: 
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Highway  Class  Point  Range 

Southern  Ontartio 

Freeway  Over  20 

Major  Trunk  6  to  20 

Minor  Trunk  1  to  5 

Northern  Ontario 

Major  Trunk  Over  3 

Minor  Trunk  1  to  3 

A  measure  of  rural  access  service  was  needed  to  ensure 
that  persons  living  in  rural  areas  had  reasonable  access  to 
Trunk  highways.   The  principle  employed  was  that  areas  of 
similar  population  density  should  be  similarly  served. 

Areas  completely  bounded  by  Trunk  highways  and,  in  some 
cases,  by  shorelines  were  called  highway  cells.   Highway 
accessibility  was  measured  as  the  distance  in  miles  from 
the  most  remote  point  'v-ithin  the  cell  to  the  bounding  high- 
ways.  This  distance,  when  combined  with  the  cell's  popu- 
lation density,  provided  a  measure  of  under-  or  over-service 
to  particular  areas. 

Detailed  study  then  determined  whether  a  route  should 
be  added  or  removed  from  the  system  such  that  the  system 
would  retain  its  continuous,  integrated  character. 

VJashington 
The  State  of  Washington's  classification  procedure  was 
composed  of  two  parts.   The  first  was  an  analysis  designed 
to  yield  an  Index  of  State  interest  in  highways.   The  main 
purpose  of  the  Index  was  to  indicate  whether  or  not  the  road 
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section  could  be  considered  of  sufficient  State  interest 
to  be  on  the  State  highway  system  (21) , 

The  Index  was  a  composite  weighted  average  of  the  fol- 
lowing six  factors: 

1.   Intercity  travel  desire  factor  (ITDF) .    All  cities 
of  over  1000  population  in  Washington  and  v:ithin  20  miles 
of  its  borders  were  included  in  the  study.   A  basic  factor 
was  first  obtained  for  all  possible  city  pairs  by  taking 
the  square  root  of  the  product  of  the  populations  and  di- 
viding by  the  square  of  the  distance  between  them  by  the 
most  direct,  feasible  route.   The  basic  factors  assigned  to 
a  particular  road  section  v/ere  summed,  the  result  being 
called  the  intercity  travel  desire. 

The  assumption  was  then  made  that  if  the  minimum  aver- 
age annual  daily  traffic  (AAD)  on  a  road  section  could  be 
determined,  it  would  be  an  adequate  measure  of  through 
traffic  with  local  influences  essentially  eliminated.   If 
this  assumption  is  correct,  tnen  the  AAD  of  a  road  section 
should  be  highly  correlated  with  its  calculated  intercity 
travel  desire.   A  correlation  coefficient  of  0,9  3  was  ob- 
tained from  data  at  65  locations  on  7  cross-state  routes. 
The  correlation  was  performed  with  both  variables  in  log 
form  and  since  the  b-coef f icient  was  about  0.59,  the  re- 
lationship was  non-linear. 

The  ITDF  for  each  road  section  v;as  calculated  by  multi- 
plying the  basic  factor  assigned  to  a  road  section  by  the 
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length  of  the  trip.  When  this  had  been  done  for  all  basic 
factors  assigned  to  a  road  section,  the  sum  of  these  modi- 
fied factors  was  called  the  ITDF. 

The  multiplication  by  trip  length  was  performed  to  take 
into  account  the  fact  that  state  interest  increases  with 
trip  length. 

2.  Agricultural  usage  factor.   This  factor  provided  a 
measure  of  the  intensity  of  road  usage  generated  by  agricul- 
ture in  terms  of  the  value  of  farm  products  in  the  area 
served  by  the  road. 

The  area  served  by  each  road  was  first  defined  and 
the  land  assigned  a  productivity  value  as  determined  from 
data  of  the  U.  S.  Agriculture  Census.   The  land  area  in 
square  miles  was  multiplied  by  the  productivity  value.   The 
agricultural  usage  factor  v;as  this  product  divided  by  the 
length  of  the  road, 

3.  Logging  factor.   Forest  products  are  '^■Jashington '  s 
major  industry  and  are  especially  dependent  upon  highway 
transportation  inasmuch  as  over  half  the  cost  of  log  pro- 
duction is  a  transportation  charge.   The  logging  factor 
was  computed  as  the  number  of  million  board  feet  of  lumber 
hauled  over  a  road  section  during  a  year. 

4.  Mining  factor.   Mining  products  in  Washington  are 
quite  variable  both  in  magnitude  and  value.   Thus,  it  was 
decided  to  compute  the  mining  factor  as  the  total  payroll 
of  all  mines  served  by  a  road  section  divided  by  the  length 
of  the  road  section. 
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5.  Motor  freight  factor.   This  factor  was  taken  as  the 
average  freight  tonnage  per  day  hauled  over  a  specific  route. 
It  was  designed  to  take  into  account  long  distance  com- 
mercial trucking. 

6.  Recreation  factor.   Tourism  is  Washington's  fourth 
largest  industry  and  is  the  most  dependent  on  a  good  highway 
system.   A  basic  factor  was  computed  by  taking  the  square 
root  of  the  product  of  the  population  of  a  city  and  the 
number  of  visitors  to  a  major  recreational  area  less  1500 
divided  by  the  straight  line  distance  between  the  city  and 
the  area.   Only  cities  with  populations  of  over  2000  and 
recreational  areas  v^ith  over  1500  annual  visitors  were  con- 
sidered in  the  study. 

The  basic  factor  was  assigned  to  the  shortest,  most 
feasible  route  and  the  recreation  factor  was  taken  as  the 
sum  of  all  basic  factors  assigned  to  a  particular  road 
section. 

The  Index  of  State  interest  was  determined  from  the 
above  six  factors  in  tv;o  steps.   It  was  first  necessary  to 
express  each  factor  as  a  per  cent  because  of  the  diversity 
in  their  units  and  magnitude.   A  m.ean  maximum  value  was 
calculated  for  each  factor  as  the  arithmetic  average  of  the 
highest  values  observed  for  the  factor  on  any  100  miles  of 
road  in  the  State.   Each  factor  was  then  divided  by  the 
mean  maximum  value  for  all  road  sections. 

Each  factor  then  had  to  be  weighted  to  account  for 
their  unequal  contribution  to  State  interest.   The  existing 
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State  system  was  used  to  provide  the  proper  weights  which, 
when  finally  selected,  resulted  in  91  per  cent  conforinity 
to  the  existing  system. 

The  Index  was  used  to  measure  a  route's  eligibility 
for  consideration  as  a  State  highv/ay.   A  chart  ;vas  prepared 
showing  the  cumulative  miles  with  an  Index  equal  to  or 
higher  than  the  Index  being  plotted.   The  Index  at  which  the 
cumulative  mile  total  equalled  the  present  State  system  was 
20  which  was  assumed  to  be  the  cutoff  point  for  State 
interest. 

A  later  study  in  19  6  2  employed  the  Illinois  method  to 
classify  the  State  system  (excluding  the  Interstate  system) 
into  four  groups.   All  cities  of  over  1000  population  were 
ranked  by  population  and  successively  interconnected.   All 
highways  needed  to  interconnect  cities  of  over  20,000  popu- 
lation were  called  Principal  State  highways  while  those 
needed  to  interconnect  the  remaining  cities  were  called 
Major  State  highways. 

Collector  State  highways  were  those  needed  to  provide 
rural  access  service  and  service  to  recreational  areas. 
Other  State  highways  were  essentially  those  which  had  a  high 
enough  Index  to  be  included  in  the  State  system. 

This  classification  technique  was  designed  to  provide 
a  basis  for  priority  programming  of  improvements. 
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New  Mexico 

The  method  of  highway  classification  proposed  in  New 
Mexico  consists  of  two  steps  sinilar  to  those  employed  in 
Washington.   However,  the  method  of  applying  them  to  classi- 
fication was  different  (23). 

The  first  step  employed  the  Illinois  method  to  select 
a  basic  highway  system.   The  State  was  divided  into  immediate 
trade  areas  each  served  by  a  central  place.   Population, 
bank  resources,  newspaper  circulation,  and  sales  were  used 
to  rank  cities.   Five  city  groups  known  as  regional  centers, 
major  market  centers,  market  center  A,  market  center  3  and 
minor  market  centers  were  selected. 

A  basic  highway  system  was  selected  and  each  highway  was 

classified  by  the  type  of  city  connection  it  provided  as 

follows : 

Primary  Regional-Regional 

Major    -'4ajor 
Regional -Ma j or 

Secondary  A  !4arket-  ".arket 
Regional -Market 
Major    -Market 

Secondary  B  Minor  -Minor 
Regional -Mi nor 
Market   -Minor 

This  technique  is  called  the  Functional  Approach  and  it 
results  in  what  is  called  the  minimum  basic  road  system. 

The  second  step  involved  the  calculation  of  an  Index 
for  all  roads  to  provide  a  measure  of  State  interest.   The 
factors  considered  in  developing  the  numerical  Index  were: 
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intercity  travel  desire  factor  (including  in-state  and  out- 
of-state  phases) ,  recreational  factor  (including  in-state 
and  out-of-state  phases) ,  agricultural  factor,  mining  fac- 
tor, and  an  oil  and  gas  factor. 

The  criteria  leading  to  the  selection  of  these  factors 
were : 

1.  Relative  importance  to  the  economy  of  the  State. 

2.  Numerical  nature  of  the  pertinent  data. 

3.  Availability  and  reliability  of  data. 

4.  Adaptability  to  a  numerical  index. 

5.  Flexibility  of  computed  results. 

The  intercity  travel  desire  factor  was  computed  in 
essentially  the  same  manner  as  in  Washington  but  included 
an  out-of-state  phase  to  account  for  the  significant  use  of 
New  Mexico  highways  as  corridors .   The  recreational  factor 
was  derived  in  tv;o  phases  for  the  same  reason.   A  special 
study  v;as  performed  at  a  national  monument  and  the  recrea- 
tional factor  derived. 

The  agricultural,  mining,  and  oil  and  gas  factors  were 
computed  similarly  to  the  Washington  method  but  modified 
to  reflect  the  type  of  data  available  in  New  Mexico. 

The  factors  then  were  weighted  and  combined  to  give  the 
numerical  Index  for  each  road  section.   The  next  step  v/as 
to  determine  the  Index  cutoff  points  for  each  road  classi- 
fication so  that  all  roads  not  already  classified  could  be 
classified. 
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A  synthesis  of  the  Functional  Approach  v/as  devised 
based  upon  the  premise  that  the  road-user  services  and  rela- 
tive importance  implicit  in  the  Functional  Approach  categor- 
ies may  be  described  explicitly  by  the  relative  standing 
of  the  numerical  Indices  for  all  roads  in  the  minimum  basic 
system.   The  resulting  makeup  of  relative  standings  can 
then  be  extended  to  all  highway  sections  not  included  in 
the  minimum  basic  system  by  use  of  their  numerical  Indices. 

A  method  using  modal  classes  and  arithmetic  averages 
was  devised  to  select  the  most  representative  numerical 
Index  of  each  functional  category.   This  data  v/as  fitted 
to  a  negative  exponential  curve  and  six  functional  classes 
were  selected,  the  last  three  of  which  were  Other  1,  Other 
2,  and  Other  3. 

Comparisons  of  the  numerical  Indices  of  roads  not 
classified  by  the  Functional  Approach  with  the  Index  limits 
of  each  of  the  six  functional  classes  automatically  provided 
the  functional  classification  of  these  roads. 

This  method  resulted  in  the  recommendation  of  additions 
to,  and  deletions  from,  the  existing  system.   Hov/ever,  be- 
fore the  recommendations  could  be  implemented,  they  had  to 
be  considered  in  conjunction  with  the  systems  of  adjacent 
states,  jurisdictional  responsibility,  fund  apportionment 
policy,  the  financial  structure  of  the  state,  and  laws  al- 
ready in  effect  or  needed  to  implement  the  proposed  system. 
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STUDY  PROCEDURE 

The  principal  aims  of  this  study  xvere  to  subclassify 
those  highways  of  primarily  statewide  interest  in  Indiana 
and  to  pinpoint  those  which  should  be  considered  for  recon- 
struction to  freeway  standards.   The  method  selected  was  a 
statewide  study  of  intercity  travel  desire. 

This  method  was  selected  prim.arily  because  its  use  would 
result  in  the  assignment  of  a  numerical  factor  to  each  high- 
way section  in  the  State.   The  magnitude  of  the  factor 
would  be  a  measure  of  intercity  travel  desire  and,  hopefully, 
a  measure  of  both  the  relative  importance  and  anticipated 
traffic  volume  of  each  highway  section.   It  was  further  en- 
visioned that  some  function  of  the  importance  of  a  highway 
section  and  its  overall  condition  as  measured  by  a  suffic- 
iency rating  could  be  used  to  establish  a  set  of  relative 
priorities  for  construction,  reconstruction,  and  maintenance. 

It  was  decided  to  calculate  the  Intercity  Travel  Desire 
Factors  (ITDF)  in  the  framework  of  a  computer  oriented 
mathematical  model.   It  was  reasoned  that  this  method  would 
be  somewhat  more  objective  in  character  than  some  of  the 
methods  described  in  the  previous  chapter.   Even  the  least 
subjective  methods  described,  those  of  Washington  and  New 
Mexico,  depended  on  the  researcher's  judgment  with  respect  to 
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the  most  feasible  intercity  route.   It  was  the  intention  in 
this  study,  however,  to  determine  the  most  feasible  inter- 
city route  based  on  the  minimum  path  distance  through  the 
h  i  ghw  ay  ne  tavo  rk  . 

Distance  as  used  here  and  throughout  the  balance  of  this 
report  when  referring  to  minimum  paths  is  not  a  true  dis- 
tance measure.   This  v/as  due  to  the  fact  that  highway  sec- 
tions built  to  less  than  Interstate  standards  were    penalized 
by  assuming  a  lower  average  speed. 

The  average  speed  on  Interstate  highways  was  taken  as 
60  MPH  or  one  mile  per  minute.   The  length  in  miles  of  an 
Interstate  section  was  thus  numerically  equal  to  the  time 
in  minutes  required  to  traverse  the  section.    The  lengths 
of  all  other  highway  sections  were  multiplied  by  a  factor 
to  convert  distance  to  time  in  minutes.   The  factor  was  com- 
puted as  the  ratio  of  60  MPH  over  the  assumed  speed  for  the 
highway  section. 

Intercity  Travel  Desire  Factor 
The  Intercity  Travel  Desire  Factor  used  in  this  study 
is  based  on  the  gravity  concept  of  human  interaction.   The 
basic  postulate  of  the  gravity  concept  is  that  the  interac- 
tion between  cities  varies  directly  with  some  function  of 
the  populations  of  the  two  cities  and  inversely  with  some 
function  of  the  distance  between  them.   This  may  be  expressed 
mathematically  as: 
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ITDF.  .  -  f(P^,  P.)/f(d.  .) 

where  P^  and  Pj  are  the  populations  of  the  two  cities  and 
d^  .  is  the  distance  between  them, 

A  great  deal  of  effort  has  been  expended  in  developing 
the  gravity  concept  of  travel  synthesis  especially  in  the 
area  of  determining  the  functions  to  be  used  in  the  above 
expression.   A  detailed  summary  of  much  of  the  important 
work  in  this  area  may  be  i'ound  in  a  report  by  Tittemore  (51)  , 

The  form  of  the  expression  used  in  this  study  was  taken 
as : 

ITDF.j  =  (PiPj)^'^^/dij^ 

This  is  the  same  form  for  the  Intercity  Travel  Desire  Fac- 
tor as  was  used  in  the  studies  of  I'Jashington  (21)  and  New 
Mexico  (23)  as  reported  in  the  previous  chapter. 

In  the  Washington  study,  several  formulations  of  the 
travel  desire  factor  were  tried  and  each  was  correlated  with 
the  minimum  annual  average  daily  traffic  on  roads  connectina 
population  centers  (52).   The  form  of  the  factor  listed 
above  provided  the  best  results,  the  correlation  coefficient 
being  0.93  (21) , 

Several  limitations  regarding  the  use  of  the  gravity 
concept  have  been  pointed  out  (51).   The  principal  one  con- 
cerns the  saturation  of  a  city's  travel  desire.   An  upper 
limit  exists  on  the  number  of  trips  that  can  be  made  in  a 
given  time  period  by  each  person.   VJhen  this  limit  is  reached, 
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the  relationship  between  travel  desire,  population,  and  dis- 
tance is  no  longer  valid.   The  exponent  of  distance  has  been 
shown  to  vary  both  as  a  function  of  distance  and  trip  pur- 
pose.  Some  work  has  also  been  done  in  weighting  the  popu- 
lations of  cities  to  take  differences  in  income  and  car 
registration,  among  others,  into  account. 

Procedural  Steps 
The  procedure  followed  in  the  performance  of  this 
research  consisted  of  the  following  steps: 

1.  Coding  of  the  highway  network  and  the  location  of 
all  cities  to  be  used  in  measuring  Intercity  Travel  Desire 
Factors . 

2.  Determination  of  the  minimum  path  distance  between 
each  city  pair.   The  large  size  of  the  network  required  the 
development  of  a  tree  type  decomposition  algorithm  for 
finding  the  minimum  paths  in  large  networks. 

3.  Assignment  of  the  computed  travel  desire  factor 
between  two  cities  to  each  link  making  up  the  minimum  path 
between  the  two  cities.   The  travel  desire  factor  was  com- 
puted as  the  square  root  of  the  product  of  the  populations 
of  the  two  cities  divided  by  the  square  of  the  minimum  path 
distance  between  them. 

4.  Determination  of  the  adequacy  of  the  calculated 
Intercity  Travel  Desire  Factors  to  measure  traffic  by  means 
of  regression  analysis.   The  dependent  variable  was  the  sum 
of  the  average  daily  traffic  measurements  on  all  roads  going 
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in  and  out  of  a  city.   The  independent  variables  were  mea- 
sures of  the  Intercity  Travel  Desire  Factors. 

5.   Classification  of  the  state  highway  system  of 
Indiana  into  several  functional  classes.   Highways  recommended 
for  reconstruction  to  freeway  standards  were  also  selected. 

The  highway  network  coded  for  this  study  consisted  of 
the  entire  state  highway  system  of  Indiana,  the  major  routes 
in  the  bordering  states,  and  the  nationwide  system  of 
Interstate  highways.   This  network  was  quite  large  consisting 
of  more  than  1800  highway  nodes  and  2700  centroids  with  a 
total  of  over  6300  links. 

This  network  was  so  large  that  it  could  not  simultan- 
eously be  stored  in  the  available  computer.   Thus,  a  decom- 
position algorithm  capable  of  computing  minimum  paths  in 
large  networks  had  to  be  developed  in  order  to  carry  out  the 
objectives  of  this  research.   The  development  of  the  algo- 
rithm is  presented  in  detail  in  the  next  chapter. 
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DEVELOPMENT  OF  A  TREE  TYPE  DECOMPOSITION 
ALGORITHM  FOR  MINIMUM  PATHS  IN  LARGE  NETWORKS 


Introduction 
A  highway  system  can  be  considered  as  a  network  con- 
sisting of  a  set  of  nodes  N  representing  route  intersections 

and  links  L.  .  representing  the  portion  of  a  route  connecting 
1  j 

nj_  and  n^  .   The  highway  systen  can  then  be  described  by  a 

set  of  links  called  a  link  table.   Each  link  is  described  by 

its  two  end  nodes  and  the  associated  link  length  d-  .  v;hich 

may  be  measured  in  terms  of  time,  cost,  distance,  etc. 

Some  links  may  be  directed,  thus  providing  for  only  one 

way  traffic  while  others  may  be  undirected  and  provide  for 

two  way  traffic.   In  the  case  of  an  undirected  link,  two 

link  descriptions  are  necessary  in  the  link  table,  one  to 

indicate  movement  from  n.  to  n.  and  the  second,  movement 

from  n.  to  n .  .   It  is  not  necessary  that  d.  .  =  d...    Hov/- 
J      1  ID     Di 

ever,  the  network  used  in  this  study  consisted  entirely  of 

undirected  links  for  which  d. .  =  d- ■ . 

ID     Di 

A  path  through  the  network  leading  from  n.  to  n .  is 
defined  by  a  sequence  of  nodes  each  of  which  can  be  adjacent 
to  only  two  links,  i.e.,  an  incoming  and  an  outgoing  link. 
The  lengths  of  each  link  along  the  path  may  be  summed  so 
that  a  shortest  path  may  be  defined  as  the  path  whose  total 
link  lengths  is  a  minimum. 
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Minimum  path  algorithms  are  available  for  determining 

the  minimum  paths  oetv^een  all  node  pairs.   The  algorithms 

may  be  divided  into  two  classes,  matrix  and  tree  methods. 

Matrix  methods  simultaneously  determine  the  shortest  route 

betv/een  all  node  pairs  in  a  netv;ork..   In  a  paper  by  Mills 

(41) ,  he  states , 

The  cascade  method,  which  is  the  most  elegant  and 
efficient  of  the  matrix  methods,  requires  simul- 
taneous storage  for  2n2  variables,  where  n  is  the 
number  of  nodes  in  the  network. 

This  amounts  to  500,000  storage  locations  for  a  network  of 

500  nodes. 

The  cascade  method  is  described  in  a  paper  by  Farbey, 
Land,  and  xMurchland  (45).   Mills  also  states  that  the  time 
required  for  solution  with  the  cascade  method  increases  with 
the  cube  of  the  number  of  nodes  (41) ,   Because  of  the  rela- 
tively large  time  and  storage  requirem.ents  for  the  solution 
of  large  networks,  some  recent  work  has  been  devoted  both 
to  shortening  the  calculations  and  the  development  of  a  de- 
composition algorithm  including  that  of  Mills  (41)  and  Hu 
(46,47)  among  others. 

Tree  methods  are  used  to  find  the  shortest  route  from 
a  source  node  to  every  other  node,  the  process  being  repeated 
until  all  nodes  have  been  used  as  sources.   Tree  methods 
require  relatively  less  storage  capacity  than  matrix  methods 
but  supposedly  are  more  difficult  to  program  and  require  a 
great  amount  of  computation  (41). 
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It  was  originally  intended  that  the  tree  building  algo- 
rithm used  by  the  Bureau  of  Public  Roads  (43)  be  used  in 
this  study  to  determine  the  minimum  paths  between  centroids. 
However,  the  coded  network  was  so  large  that  it  could  not 
simultaneously  be  stored  in  the  computer.     Thus,  it  was 
necessary  to  decompose  the  network  in  order  to  arrive  at  a 
solution. 

The  purpose  of  this  section  is  to  present  the  basic 
minimum  path  algorithm  em.ploying  the  tree  method  and  then, 
the  development  of  a  decomposition  algorithm  to  handle  a 
large  network  by  breaking  it  up  into  several  manageable 
parts , 

Basic  Minimum  Path  Algorithm 
The  initial  development  of  the  tree  building  minimum 
path  algorithm  is  generally  attributed  to  E.  F.  Moore  (42). 
In  his  paper,  Moore  presented  four  algorithms  for  finding 
the  shortest  path  through  a  maze,  the  first  three  of  which 
considered  unit  links  v;hile  the  fourth  considered  the  impo- 
sition of  a  penalty  such  as  time,  cost,  or  distance  on  each 
link.   This  fourth  algorithm  was  later  modified  and  program- 
med as  the  basic  tree  building  process  used  in  the  Traffic 
Assignment  package  as  provided  by  the  U.  S.  Bureau  of  Public 
Roads  (43)  . 

The  essential  operation  of  the  algorithm  is  one  of 

comparison.   Consider  any  three  nodes  n. .  n,,  and  n  .   The 

I'D        k 

distance  d- •   is  considered  minimum  only  if  it  is  equal  to 
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or  smaller  than  the  sum  d. ,  +  d,  . .   The  length  of  a  one  link 
path  is  compared  with  the  length  of  a  path  through  an  inter- 
mediate node  nj^  and  the  original  path  is  changed  only  if  its 
length  is  the  larger  of  the  two  distances  (or  costs,  times, 
etc.) .   The  process  may  then  be  extended  to  the  comparison 
of  path  lengths  where  each  path  may  be  a  sequence  of  sev- 
eral nodes . 

Any  network  v;hich  has  a  finite  number  of  nodes  and 
links  must  also  have  a  finite  number  of  paths  connecting 
any  two  nodes  at  least  one  of  v/hich  is  a  minimum  path.   The 
problem  is  one  of  finding  an  efficient  computational  pro- 
cess which  can  find  this  minimum  path  directly  rather  than  by 
trial  and  error. 

Assume  a  network  exists  as  shown  in  Figure  5-a  and  that 
the  minimum  path  from  a  to  b  is  desired.   The  first  step  was 
to  find  the  link  out  of  node  a  which  had  the  smallest  length 
and  draw  an  arrov;  from  node  a  to  the  node  connected  by  this 
link.   In  this  example,  the  arrow  was  drawn  to  node  c  and  a 
'2'  written  alongside  node  c  as  shov;n  in  Figure  5-b. 

The  second  step  was  to  find  the  link  out  of  a  node  al- 
ready reached  (either  a  or  c)  to  a  node  not  yet  reached 
(either  d,  e,  or  h)  such  that  the  total  distance  back  to 
node  a  was  a  minimum.   At  this  point,  an  arrow  was  drawn 
from  node  c  to  node  d  and  a  '5'  written  alongside  node  d 

(v/here  '5'  is   the  sum   of  link  lengths  d    and  d   ).   The 

ac       cd 

v/ritten  number  is  referred  to  as  a  label.   This  process 
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FIGURE     5     -a,-b,-c.  A  SMALL  NETWORK   ILLUSTRATING 

ALTERNATE    MINIMUM   PATH  ALGORITHMS 
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continued  until  all  nodes  v/ere  reached  and  the  minimum  paths 

to  all  nodes  were  found  as  shown  in  Figure  5-b. 

At  any  point  in  this  process,  the  reached  nodes  may  be 
considered  as  a  set  S  where  each  node  is  labeled  with  a 
distance  v.,  jeS  (j  included  in  S)  which  is  the  minimum 
path  distance  back  to  the  starting  node.   All  unlabeled  nodes 

make  up  the  complement  set  S-   A  node  i  is  added  to  S  and 

_  * 

deleted  from  S  provided  that: 


V   =  min  (v   +  rain  d   )  [1] 

i    jcS    J    icS"   Ji 


This  process  is  time  consuming  and  relatively  inef- 
ficient for  computer  application  with  respect  to  the  fol- 
lov'ing  alternate  which  is  essentially  that  used  by  the 
Bureau  of  Public  Roads  (43).   The  first  step  of  the  alternate 
was  to  label  all  nodes  that  could  be  reached  from  the  start- 
ing node  with  the  associated  link  length.   This  differs 
from  the  first  process  in  that  the  labeled  distance  may  or 
may  not  be  the  minimum  distance.   The  second  step  H'as  to 
search  the  labeled  nodes  and  find  the  node  having  the  mini- 
mum label.   In  Figure  5-c,  the  labeled  nodes  at  this  point 
were  c,  d,  and  e  and  node  c  had  the  minimum  label.   Next, 
for  each  node  that  could  be  reached  from  c,  the  sums  of  the 
label  on  c  and  the  associated  link  length  were  obtained. 


Some  of  the  preliminary  detail  and  notation  in  this  section 
is  due  to  reference  44. 
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These  sums  were  then  compared  with  the  labels  on  the 
newly  reached  nodes,  in  this  instance  d  and  h.   If  the  node 
was  labeled,  the  label  was  replaced  only  if  the  sum  was 
less  than  the  existing  label,  the  case  for  node  d  in  Figure 
5-c.   If  the  node  was  unlabeled,  it  was  labeled  with  the  sum, 
the  case  for  node  h. 

The  process  then  reverted  to  the  second  step  to  find 
the  minimum  label  on  the  nodes,  d,  e,  and  h  and  continued 
until  no  further  label  changes  could  be  m.ade .   The  completed 
minimum  path  routings  by  this  process  are  shown  in  Figure 
5-c. 

The  arrows  are  reversed  in  Figure  5-c  because  the 
paths  are  traced  back  from  each  node  to  the  starting  node. 
In  tracing  the  paths,  the  only  restriction  is  that  the  dif- 
ference in  labels  of  two  adjacent  nodes  must  equal  the  link 
length . 

In  both  methods,  however,-  the  paths  and  distances  are 
identical  and  the  paths  from  all  nodes  to  the  starting  nodes 
were  found. 

In  the  second  method  it  is  convenient  to  label  all  nodes 

with  an  arbitrarily  large  number  at  the  start  except  that 

the  starting  node  is  assigned  a  zero  label.   Then  as  the 

process  continues,  all  nodes  for  which  the  labels  change  are 

placed  in  the  S,  all  other  nodes  in  the  complement  set  S- 

The  label  v.  of  a  node,  i,  is  changed  and  the  node 
1 


added  to  S  provided  that; 
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V .  >  min  V .  +  d  .  .  [21 

where  i  can  be  a  node  already  in  S.   The  node,  j,  is  removed 
from  S  when  all  nodes  that  can  be  reached  from  it  have  been 
checked  as  in  [2]  above.   The  process  is  completed  when  S 
is  empty  so  that  no  further  label  changes  are  possible. 

The  operation  indicated  by  [2]  may  be  expressed  in  a 
slightly  different  form  and  as  such  is  presented  in  the 
following  theorem. 

Theorem  I .   The  minimum  path  distance  to  node  i  (not  equal 
to  j)  must  be  equal  to  or  less  than  the  path  distance  by 
way  of  node  j,  i.e.: 

V.  --  V  .  +  d  .  .  [3] 

and  is  equal  to  the  minimum  of  all  sums  obtained  by  letting 
j  vary  over  all  nodes  directly  connected  to  i,  i.e.: 


V .  -  min  (v .  +  d . . )  [41 


The  sequence  of  nodes  defining  the  minimum  path  from 
node  i  to  the  starting  node  must  include  one  and  only  one 
node  directly  linked  to  node  i.   Thus,  if  the  minimum  paths 
to  all  nodes  j  directly  connected  to  node  i  are  known,  the 
minimum  path  for  node  i  is  defined  by  [4] . 

Suppose  the  minimum  path  between  any  two  nodes,  a  and 

b,  in  a  network  is  known.   At  any  intermediate  node  j  on 

the  path,  let  v  .  be  the  minimum  path  distance  to  a  and  v,  . 
^         '  aj  ^  bj 
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the  minimum  path  distance  to  b.   If  v  ,  is  the  minimum  path 
^  ab  '^^ 

distance  from  a  to  b ,  then  it  m.ust  follow  that: 

V  ,  =  V  .  +  V,  .  [  5  ] 

ab     aj     bj 

If  V  .  and  V,  .  are  defined  as  the  minimum  path  distances 
ai       bi  ^ 

to  nodes  a  and  b  respectively  from  any  node  i  in  the  network, 

their  sum  would  be  greater  than  v  ,  if  the  node  i  did  not  lie 
^  ab 

on  the  minimum  path. 

This  leads  to  the  following  theorem: 

Theorem  II.   If  a  set  of  nodes  T  is  selected  such  that  it  is 
impossible  to  form  any  path  from  node  a  to  node  b  without  at 
least  one  node  included  in  T,  the  minimum  path  from  a  to  b 
is  given  by: 


V  ,  =  min  (v  .  +  V,  .  )  [6] 


Lemma  I .   If  the  minimum  paths  from  all  nodes  (not  including 
b)  in  a  network  to  node  a  and  the  minimum  paths  from  all 
nodes  (not  including  a)  in  a  netv/ork  to  node  b  are  known, 
then  the  minimum  path  distance  from  a  to  b  is  given  by: 

V  ,  =  min  (v  .  +  v,  .  )  [7] 

ab         aj     bj 

where  j  varies  over  all  nodes  excluding  a  and  b.   Furthermore, 
the  minimum  path  from  a  to  b  may  be  traced  by  connecting  the 
sequence  of  nodes  j  for  which  [7]  holds. 

This  must  follow  from  [5]  because  node  j  can  be  any 
node  on  the  minimum  path  sequence  from  a  to  b .   This  also 
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serves  to  prove  Theorem  II  since  at  least  one  node  in  T 
must  be  on  the  minimum  path  sequence  of  nodes.   Lemma  I  may 
be  generalized  for  any  two  nodes  in  the  network  provided  that 
they  are  not  connected  by  a  direct  link. 

A  Partitioned  Network  of  Two  Parts 

Consider  a  network  which  has  been  separated  along  a 

set  of  border  nodes  T  such  that  it  is  impossible  to  go  from 

one  side  to  the  other  without  passing  through  one  of  the 

nodes  in  T.   Let  the  nodes  on  either  side  of  the  border  belong 

in  sets  A  and  B  respectively. 

In  the  network  AT  composed  of  the  sets  A  and  T,  the 

minimum  path  distances  from  each  node  in  T  to  all  nodes  in 

A  may  be  determined.   The  same  is  true  for  all  nodes  in  B 

when  considering  the  network  BT  composed  of  node  sets  B  and 

T. 

Let  w.   be  the  minimumi  path  distance  within  AT  from 

node  j  in  T  to  node  a  in  A  and  w.,  the  minimum  path  distance 

3b 

within  BT  from  the  node  j  in  T  to  node  b  in  B.   Then,  in  a 
manner  similar  to  Theorem  II,  the  minimum  path  distance  from 
a  to  b  should  be  approximately: 

V  ,  =  min  (w.   +  w.,  )  18] 

ab    .  ^    la     jb 

as  shown  in  Figure  6 . 

Equation  [8],  written  as  an  equality,  can  be  shown 
immediately  to  be  false  if  the  true  minimum  path  from  a  to 
b  is  as  shown  by  the  solid  line  in  Figure  7.   The  seeming 
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FIGURE    6.     A  PARTITIONED  NETWORK  OF  TWO  PARTS. 


Ti,i=l,5 

FIGURE  7.      ALTERNATIVE   PATHS    THROUGH  A  PARTITIONED 
NETWORK. 
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inconsistency  with  Theorem  II  develops  from  a  difference  in 

definition  of  v.   and  w.  .   In  Theorem  II,  v.   was  defined 
ja       ja  ja 

as  the  minimum  path  distance  from  node  j  in  T  to  node  a  using 
the  entire  netv/ork.   In  [8]  ,  however,  w,   is  defined  simi- 
larly  but  only  a  part  of  the  network  is  used.   As  a  result, 
some  nodes  that  lie  along  the  minimum  path  node  sequence  may 
not  be  available  so  that  only  an  apparent  minimum  path  is 
found  when  less  than  the  entire  network  is  used. 

In  Figure  7,  the  path  starts  at  b  and  goes  to  T^  by 
way  of  nodes  i  in  B,  from  T^  to  T.  by  way  of  nodes  i  in  A, 
from  T.  to  T^  by  way  of  nodes  i  in  B,  and  from  T_^  to  node  a 
by  way  of  nodes  i  in  A.   However,  if  the  true  minimum  path 
proceeds  from  b  to  T^-  to  T.  to  a  as  shown  in  Figure  7,  then 
equation  [8]  will  give  the  true  minimum  path  with  j  equal 
to  5  since  the  entire  m.inimum  path  from  T^-  to  a  is  included 
in  the  node  sets  A  and  T.   Thus,  a  problem  arises  when  the 
true  minimum  path  crosses  the  node  set  T  at  least  three 
times  when  the  origin  and  terminal  nodes  are  on  opposite 
sides  of  the  node  set  T. 

A  similar  problem  arises  when  considering  origin  and 
terminal  nodes  on  the  same  side  of  T  such  as  nodes  a  and  c 
of  Figure  5.   The  minimum  path  from  a  to  c  may  be  obtained 
directly  by  the  basic  minimum  path  algorithm  using  node  sets 
A  and  T  provided  that  no  portion  of  the  true  minimum  path 
utilizes  any  of  node  set  B.   This  is  illustrated  in  Figure  7 
where  the  path  is  from  c  to  T   to  T^  to  a.   If  the  T   to  T^ 
portion  of  the  path  is  contained  entirely  in  a,  no  problem 
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arises . 

Thus,  an  algorithm  devised  to  handle  a  partitioned 
network  must  be  capable  of  handling  both  these  situations. 

The  following  theorem  and  its  corollary  will  prove  of 
value. 

Theorem  III.   If  the  minimum  path  between  nodes  on  opposite 
sides  of  T  (such  as  nodes  a  and  b)  contains  at  most  two 
nodes  in  T,  the  minimum  path  may  be  found  directly  from 


v^v  =  min  (w.   +  w.,  )  [9] 

3b    j^T    :a     Db' 


With  three  or  more  nodes  in  T,  the  minimum  path  may  be  found 
directly  from  [9]  only  if  the  entire  portion  of  the  path 
along  the  nodes  in  T  lies  on  only  one  side  of  T. 

Corollary  I.  If  the  minimum,  path  between  nodes  on  the  same 
side  of  T  (such  as  a  and  c)  contains  at  most  one  node  in  T, 
the  minimum  path  may  be  found  by  the  basic  minimum  path 
algorithm  using  the  appropriate  node  set  and  T.  With  two  or 
more  nodes  in  T,  the  basic  minimum  path  algorithm  will  work 
only  if  the  entire  path  is  on  the  same  side  of  T  as  the  two 
nodes . 

Decomposition  Algorithm  for  A 
Partitioned  Network  of  Two  Parts 

The  proposed  algorithm  consists  of  several  steps  as 

outlined  below. 
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1.  Divide  the  network  into  two  parts  along  a  node  set 
T  such  that  it  is  impossible  to  pass  from  one  part  to  the 
other  without  going  through  at  least  one  node  in  T.   The  node 
sets  on  either  side  of  T  are  A  and  B  respectively. 

2.  Prepare  the  link  tables  for  node  sets  A  and  T  and 
for  B  and  T.   If  adjacent  nodes  in  T  are  directly  connected, 
the  link  descriptions  should  be  assigned  to  one  but  not  both 
of  the  resulting  link  tables. 

3.  Use  the  basic  minimum  path  algorithm  with  link 

table  AT  and  calculate  the  minimum  paths  to  all  nodes  in  AT 

from  all  nodes  in  T.   From  these  tree  tables,  generate  a 

two  dimensional  array,  w   ,  representing  the  distance  from 

each  node  in  T  to  all  nodes  in  A.   A  second  two-dimensional 

array,  w    ,  representing  the  distance  from  each  node  in  T 
t  ta 

to  all  nodes  in  T  by  way  of  nodes  in  AT,  should  be  generated 
as  well. 

4.  Repeat  step  3  with  link  table  BT,  generating  the 

distance  arrays  w^,  and  w^^,  . 
■'         tb       ttb 

5.  Use  the  basic  minimum  path  algorithm  with  link 
table  BT  and  calculate  the  minimum  path  tree  table  with  a 
particular  starting  node  in  BT . 

6.  Introduce  the  w     array  of  step  3  at  this  point. 

For  the  first  node  i  in  T,  sum  the  distance  w.,  to  the 

ik 

starting  node  k  in  BT  (from  the  tree  table  of  step  5)  and 

the  distance  w. .   to  node  i  in  T  (from  the  w^^   array) .   Com- 
xja  -^  tta      ^ 

pare  this  sum  with  w.,   (the  label  on  j)  and  if: 
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W  .,   >  W.  ,   +  W.  . 

jk     ik     ija 

replace  the  label  on  j  with  the  sum.  and  place  the  node  j  in 
the  set  S.   Continue  this  process  over  all  nodes  j  in  T. 
Repeat  the  entire  step  for  all  nodes  i  in  T. 

7.  If  the  set  S  is  empty,  go  to  step  9.   Otherwise, 
use  the  basic  ninimum  path  algorithm  v/ith  the  link  table  of 
step  5  to  make  any  changes  in  labels  on  nodes  in  BT  made 
necessary  by  label  changes  on  nodes  in  T  in  step  6.   Continue 
until  S  is  empty.   If,  during  the  operations  of  this  step, 

no  label  changes  were  made  for  any  node  in  T,  go  to  step  9. 

8.  Repeat  steps  6  and  7  until  conditions  are  such 
that  transfer  to  step  9  is  allowed. 

9.  Introduce  the  array  w    from  step  3.   For  the 
first  node  i  in  T,  sura  the  distance  w.,  from  the  starting 

XX 

node  k  to  i  (from  the  tree  table)  and  the  distance  w.   from 

m 

i  to  node  n  in  A  (from  the  w    array) .   Compare  this  sum 

t  a     ■* 

with  the  label  w,   on  n  (initially  set  to  an  arbitrarily 
kn  ^  -^ 

large  value)  and  if: 

w,   >   \j . ,     +  w . 
kn     ik     m 

replace  the  label  on  n  with  the  sum.   Continue  this  process 
over  all  nodes  n  in  A.   Repeat  the  entire  step  for  all  nodes 
i  in  T. 
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Discussion  of  the  Algorithm 

Due  to  the  fact  that  the  principal  reason  for  using 
this  algorithra  is  that  of  computer  storage  limitations,  some 
consideration  must  be  given  to  a  computer  oriented  solution 
of  the  algorithm.   One  important  point  in  this  regard  is 
that  in  large  networks,  the  distance  arrays  generated  in 
steps  3  and  4  may  easily  exceed  the  storage  capacity  of  a 
computer  by  themselves. 

The  tree  tables  generated  in  steps  3  and  4  should  be 
stored  on  tape  so  that  the  minimum  paths  may  later  be  traced. 
The  distance  arrays  must  also  be  stored  on  tape  and  in  such 
a  way  that  the  distances  to  all  nodes  nay  be  read  individu- 
ally for  each  node  i  in  T. 

Step  6  of  the  algorithm  is  valid  due  to  Theorem  I  v/hich 
implies  a  procedure  for  comparing  a  path  to  a  node  with  a 
path  by  way  of  another  node.   Each  entry  in  the  w     array 
may  be  considered  as  a  direct  connection  between  any  two 
nodes  in  T  so  as  to  conform  with  the  wording  of  Theorem  I . 

At  the  conclusion  of  step  5,  the  tree  table  contains 

the  distances  w.,  from  each  node  in  T  to  the  starting  node 
]b 

b  as  sho\\7n  in  the  right  hand  portion  of  Figure  6.   In  step 
6,  the  distance  w.,  from  T.  to  b  (refer  to  Figure  6)  is 
compared  to  the  distance  from  T.  to  b  by  way  of  Tj.  as  shown 
in  the  right  hand  portion  of  Figure  7. 

In  step  7,  the  operations  of  the  basic  minimum  path 
algorithm  eventually  makes  a  comparison  of  the  distance 
from  T^  to  b  with  the  distance  from  T-.  to  b  by  way  of  T. 
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and  T|-  as  shown  in  Figure  7.   Step  6  must  then  be  repeated 
so  that  a  comparison  of  say  T_  to  b  versus  T^  to  ^   ^Y  ^^^Y 
of  T  ,  T.,  and  T  may  be  made.   Eventually  all  possible 
comparisons  will  be  made  so  that  at  the  beginning  of  step  9, 
the  true  minimum  paths  from  all  nodes  in  BT  to  the  starting 
node  are  in  the  tree  table. 

Prior  to  step  9,  the  tree  table  should  be  stored  on 
tape  and  the  next  starting  node  read.   Repeat  steps  5  through 
8  for  each  starting  node. 

The  link  table  BT  will  not  be  needed  to  complete  the 
tree  table  so  that  the  storage  space  it  requires  is  avail- 
able for  the  operations  of  step  9. 

The  validity  of  step  9  is  due  to  Theorem  II  which 
states  that  a  minimum  path  between  nodes  on  opposite  sides 
of  T  cannot  exist  without  at  least  one  node  in  T.   The  mini- 
mum path  between  two  nodes  is  defined  as  the  minimum  obtained 
by  summing  the  respective  distances  from  each  of  the  two 
nodes  to  a  node  in  T.   At  the  conclusion  of  step  9,  the 
minimum  paths  to  all  nodes  in  the  network  from  the  starting 
node  have  been  found. 

If  all  the  nodes  in  BT  have  been  used  as  starting 
nodes,  then  all  paths  from  each  node  in  BT  to  each  node  in 
A  and  vice  versa  are  known.   Hence,  for  a  starting  node  in 
A  it  is  necessary  to  complete  only  steps  5  through  8  to  find 
the  minimuni  paths  to  all  other  nodes  in  A  using,  however, 
the  appropriate  link  table  and  distance  array,  i.e.,  AT  and 

^ttb- 
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Some  additional  considerations  concern  the  tracing 
of  the  sequence  of  nodes  corresponding  to  the  minimum  path. 
Theorem  III  and  its  corollary  indicate  that  if  more  than  one 
node  in  T  is  on  the  minimum  path  node  sequence,  problems 
arise  in  determining  directly  on  which  side  of  T  the  minimum 
path  lies.   These  problems  are  overcome  by  the  following 
procedure . 

There  should  be  no  links  in  the  link  table  BT  used  in 
step  3  which  provide  a  direct  connection  between  nodes  in  T. 
All  links  directly  connecting  nodes  in  T  should  be  in  the 
link  table  AT  but  with  a  dummy  node  providing  the  connec- 
tion.  Thus,  for  any  direct  connection,  four  link  descrip- 
tions as  opposed  to  two  v/ill  be  required.   The  dummy  node 
is  not  considered  as  belonging  to  T. 

These  steps  insure  the  direct  determination  of  the  side 
of  T  on  which  the  path  lies.   Consider  the  tree  table  from 
step  9  for  a  starting  node  in  BT.   If  a  path  trace  contains 
no  intermediate  nodes  betv/een  two  nodes  in  T,  that  portion 
of  the  path  includes  only  nodes  from  the  set  AT.   If  inter- 
mediate nodes  do  occur,  then  that  portion  of  the  path  in- 
cludes only  nodes  from  the  set  BT. 

For  a  starting  node  in  A,  the  converse  is  true.   The 
occurrence  of  intermediate  nodes  between  any  two  nodes  in 
T  indicates  that  portion  of  the  path  incl udes  only  nodes 
from  the  set  AT.   The  dummy  nodes  serve  to  provide  the 
intermediate  node  when  a  direct  connection  is  in  the  minimum 
path  node  sequence. 
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Extension  to  a  Partitioned  Network  of  Five  Parts 

Some  networks  are  so  large  that  they  must  be  parti- 
tioned into  more  than  two  parts.   Other  networks,  especially 
highway  netv/orks ,  lend  themselves  by  geographical  considera- 
tions to  multi-partitionings . 

One  example  is  the  road  and  street  network  of  a  large 
urban  area.   The  central  portion  of  the  network  has  a  high 
node  density  because  of  the  closely  spaced  street  grid 
pattern.   However,  as  distance  to  the  outlying  areas  in- 
creases, the  node  density  becomes  relatively  low  either  be- 
cause the  roadway  grid  pattern  is  less  closely  spaced  or 
because  the  detail  required  in  the  central  area  is  not  neces- 
sary . 

Another  example  is  that  of  a  statewide  highway  network 
where  the  node-link  configuration,  although  gross  is  nonethe- 
less considerable.   The  network  detail  required  in  adjacent 
states  is  relatively  sparse  compared  to  that  of  the  state 
in  question. 

In  each  of  these  examples,  the  central  portion  may  be 
considered  as  a  natural  partitioning  v\7ith  several  partitioned 
sectors  emanating  from  it  as  shown  in  Figure  8.   Here,  the 
area  is  partitioned  into  five  parts,  each  designated  by 
its  associated  area  node  sets  J,  M,  L,  K,  and  I.   The  border 
node  sets  A,  B,  C,  D,  E,  F,  G,  and  H  are  each  selected  such 
that  it  is  impossible  to  go  from  one  part  to  another  without 
passing  through  at  least  one  of  the  border  nodes  separating 
the  parts. 
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The  magnitude  of  the  computations  involved  in  producing 
a  tree  table  is  extensive  but  depends  to  a  large  degree  on 
the  information  desired  and  the  assumptions  one  is  willing 
to  make.   Generally,  the  nodes  of  a  highway  network  fall 
into  one  of  two  classes,  namely  nodes  which  may  be  route 
intersections  and  centroids  which  may  be  any  type  of  traffic 
generator  such  as  an  entire  city  or  a  residential  zone  in  an 
urban  area.   Minimum  paths  are  needed  from  a  centroid  to  all 
other  centroids  but  not  to  all  nodes.   This  reduces  consid- 
erably the  size  of  the  distance  arrays  needed  in  the  perform- 
ance of  the  algorithm. 

The  work  is  further  reduced  if  only  path  traces  which 
utilize  nodes  in  the  central  portion  are  needed.   An  addi- 
tional savings  is  made  if  one  is  willing  to  assume  that  the 
path  between  centroids  within  a  partitioned  area  will  never 
utilize  the  nodes  of  any  adjacent  area.   This  assumption 
will  generally  not  introduce  serious  errors  if  the  borders 
are  relatively  straight  as  in  Figure  8  and  if  the  centroids 
are  not  too  close  to  a  common  border. 

At  any  rate,  the  algorithm  presented  below  outlines 
the  procedure  for  the  determination  of  the  minimum  paths 
from  a  starting  node  to  all  other  nodes  in  the  entire  net- 
work.  It  may  then  be  modified  to  suit  the  needs  of  the 
prospective  user. 
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Decomposition  Algoritlun  for  A 
Partitioned  Network  of  Five  Parts 

The  decomposition  algorithm  as  outlined  here  for  five 

parts  can  be  modified  for  n  parts.   It  is  presented  in  the 

context  of  a  five  part  partitioning  because  it  was  in  this 

form  that  the  author  solved  the  Indiana  statewide  network 

which  is  presented  in  the  next  chapter. 

1.  Divide  the  network  into  five  parts  designated  as 
area  node  sets  J,  M,  L,  K,  and  I.   Select  the  border  node 
sets  A,  B,  C,  D,  E,  F,  G,  and  H  such  that  it  is  impossible 
to  pass  from  one  part  to  another  without  going  through  at 
least  one  node  included  in  the  border  set  separating  the 
two  parts. 

2.  Prepare  the  link  tables  JBDFE ,  MABC,  LCDG,  KGFH, 
and  IHEA  using  the  associated  node  sets  in  each  case.   All 
links  furnishing  a  direct  connection  between  nodes  in  any 
of  the  border  sets  should  be  broken  by  a  dummy  node.   All 
links  of  this  type  along  any  border  set  should  be  placed  in 
only  one  of  the  link  tables  adjacent  to  the  border  but  none 
should  be  placed  in  link  table  JBDFE. 

3.  Use  the  basic  minimum  path  algorithm  to  calculate 
the  tree  tables  from  all  border  nodes  'and  then  generate  the 
distance  arrays  needed.   The  required  distance  arrays  and 
the  link  tables  used  to  obtain  them  are  listed  in  Table  1. 
The  distance  array  w,  ,  represents  the  distance  from  each  node 
in  B  to  all  nodes  in  D  while  the  array  w,  ,  .  represents  the 
distance  from  each  node  in  B  to  all  nodes  in  B  by  way  of 


76 


Table  1.   Required  Distance  Arrays 


Link  Table  Distance  Arrays" 


JBDFE  bd,  bf,  be,  df,  de ,  fe 

bbj ,  ddj ,  f f j ,  eej 
bj,  dj,  fj,  ej 

iMABC  ab,  ac,  be 

aam,  bbm,  ccm 
am ,  bm ,  cm 

LCDG  cd,  eg,  dg 

ccl,  ddl,  ggl 
cl ,  dl ,  gl 

KGFH  gf,  gh,  fh 

ggk ,  ffk,  hhk 
gk ,  fk,  hk 

IHEA  he,  ha,  ea 

hhi,  eei,  aai 
hi,  ei,  ai 


Distance  arrays  listed  in  following  order: 
Line  1-Border  node  to  border  node 
Line  2-Border  node  within  border  node 
Line  3-Border  node  to  node 
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nodes  in  JBDFE ,   In  the  case  of  border  node  to  border  node 
distance  arrays  such  as  w,  , ,  the  reverse  array  w,,  is  also 
needed.   These  arrays  are  identical  in  regard  to  individual 
distances  between  any  two  nodes.   The  transpose  of  the  w,  , 
array  is  identically  the  w,,  array. 

4.  Use  the  basic  minimum  path  algorithm  with  link 
table  JBDFE  and  calculate  the  minimum  path  tree  table  with 
a  particular  starting  node  in  J. 

5.  Introduce  the  w.   array  of  step  3  at  this  point. 

For  the  first  node  i  in  B,  sum  the  distance  w.,  to  the 

'  ik 

starting  node  k  in  J  (from  the  tree  table  of  step  4)  and  the 

distance  w.  .  to  node  j  in  A  (from  the  v;,   array)  .   Compare 

this  sum  with  w.,   (the  label  on  j)  and  if: 
jk  -* 

w  .,  >  w  . ,  -)-  w  .  . 
jk     ik     J.J 

replace  the  label  on  j  with  the  sum.   The  initial  label  on 
each  of  the  nodes  of  border  sets  A,  C,  G,  and  H,  none  of 
which  were  reached  by  the  tree  table  of  step  4,  is  set  to 
some  arbitrarily  large  value.   Continue  this  process  over 
all  nodes  j  in  A  for  all  nodes  i  in  B.   Repeat  this  entire 
step  for  each  of  the  36  distance  arrays  listed  in  their 
calling  sequence  in  Table  2  for  a  node  in  set  J.   However, 
if  the  label  for  a  node  in  a  border  set  adjacent  to  area 
J,  i.e.,  border  sets  B,  D,  E,  or  F,  is  changed  during  the 
operations  of  this  step,  that  node  is  placed  in  the  set  S. 
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Table  2. 

Sequence  o 

f  Distance  Arrays. 

Node 

Node 

Node 

Node 

Node 

in 

in 

in 

in 

in 

Set 

Set 

Set 

Set 

Set 

J 

MABC 

LCDG 

KGFH 

IHEA 

ba 

bd 

db 

fd 

ef 

be 

cd 

cb 

gd 

hf 

dc 

be 

df 

fe 

eb 

dg 

ae 

gf 

he 

ab 

fg 

bf 

de 

fb 

ed 

fh 

df 

fe 

db 

fd 

eh 

ef 

be 

eb 

bd 

ea 

ed 

bf 

ed 

bf 

ac 

de 

fb 

de 

fb 

ca 

eg 

gh 

gc 

hg 

eg 

dg 

fh 

dc 

fg 

gc 

fg 

eh 

be 

dg 

gh 

ah 

ca 

ha 

ac 

hg 

eh 

ba 

ea 

be 

ha 

fh 

ea 

ba 

dc 

ah 

db 

bd 

df 

fe 

ab 

dc 

be 

dg 

fh 

cb 

eb 

fd 

ef 

be 

cd 

ea 

fg 

eh 

da 

gd 

fb 

ed 

bf 

de 

gf 

fd 

ef 

bd 

bf 

hf 

fe 

eb 

be 

db 

he 

gc 

hg 

cd 

gh 

ae 

gd 

hf 

cd 

gf 

bbm 

gf 

he 

cb 

gd 

bbl 

ha 

ae 

ah 

ca 

ffk 

he 

ab 

ae 

cb 

aam 

hg 

ah 

ac 

gc 

ccm 

gh 

ha 

ca 

eg 

ccl 

aai 

ccra 

gel 

aam 

ggi 

bbj 

ddj 

ffj 

eej 

ggk 

ccl 

ggk 

hhi 

hhk 

hhk 

ddl 

ffk 

ddl 

ffk 

hhi 

ddj 

ffj 

ddj 

ffj 

aai 

eei 

bbm 

eei 

bbm 

ee  j 

bbi 

eej 

bbj 

ffk 

eei 

bbm 

ddl 

ffj 

eej 

bbj 

ddj 

ggk 

hhi 

ccra 

ggk 

ggl 

hhk 

ccl 

ggl 

hhk 

aai 

aam. 

ccm 

• 

hhi 

aam 

aai 

ccl 
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6.  If  the  set  S  is  empty,  go  to  step  7.   Otherwise, 
use  the  basic  minimum  path  algorithm  with  the  link  table  of 
step  4  to  make  any  changes  in  labels  on  nodes  in  JBDFE  made 
necessary  by  label  changes  on  nodes  in  B,  D,  F,  or  E  in  step 
5. 

7.  Repeat  step  5,  and  step  6  if  necessary,  until 
such  time  that  not  a  single  label  change  is  made  during  the 
entire  operation  of  step  5. 

8.  Introduce  the  border  node  to  node  distance  arrays 

in  any  desired  sequence  such  as:   am,  bm,  cm,  cl,  dl ,  gl, 

gk ,  fk.  Ilk,  hi,  ei,  and  ai .   Consider  the  w    distance 
^    '         '         >         '         '  am 

array.   For  the  first  node  i  in  A,  sum  the  distance  w.,  from 

-'  '  ik 

the  starting  node  k  to  i  (from  the  tree  table)  and  the 

distance  w.   from  i  to  node  n  in  M  (from  the  w    array) . 
m  am 

Compare  this  sum  with  the  label  w,   on  n  and  if: 

w,    V  w  . ,  +  w  . 
kn     ik     in 

replace  the  label  on  n  with  the  sum.   Continue  this  process 
over  all  nodes  n  in  M  for  all  nodes  i  in  A.   Repeat  the 
entire  step  for  each  of  the  remaining  11  border  node  to  node 
distance  arrays. 

Discussion  of  the  Algorithm 

Most  of  the  remarks  made  earlier  in  the  discussion  of 
the  decomposition  algorithm  for  a  network  of  two  parts 
apply  here  as  well. 

All  of  the  tree  tables  calculated  in  step  3  must  be 
stored  on  tape  if  it  is  necessary  later  to  trace  the  minimum 
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paths.  The  distance  arrays  must  also  be  stored  on  tape  and 
perhaps  card  punched  which  will  facilitate  the  formation  of 
any  desired  sequence. 

The  reverse  or  transpose  arrays  need  not  be  stored 
separately.   This  is  due  to  the  fact  that  a  subscripted 
element  of  the  transpose  array  is  identical  to  the  element 
in  the  original  array  with  the  subscripts  in  reverse  order, 
i.e.  , 

T 
w  .  .   -  w  .  . 

Thus,  whenever  it  is  necessary  to  use  a  reverse  array,  the 
original  array  may  be  used  with  the  subscript  order  reversed, 

A  feature  of  the  algorithm  is  that  the  true  minimum 
paths  to  all  border  nodes  has  been  accomplished  with  the 
conclusion  of  step  7. 

The  basic  method  employed  in  achieving  this  result 
involves  giving  each  and  every  border  node  the  opportunity 
to  be  the  suceeding  node  in  a  minimuiri  path  node  sequence 
to  the  starting  node.   For  instance  in  step  5  when  the  w, 
array  is  introduced,  each  node  in  B  is  tested  to  determine 
which  will  provide  the  path  linkage  from  a  node  in  A  to  the 
starting  node.   As  a  matter  of  fact,  every  node  that  can 
provide  the  path  linkage  for  a  node  in  A  is  tested  including 
nodes  in  C,  E,  H,  and  A  as  well  as  B. 

This  procedure,  the  validity  of  which  is  based  on 
Theorem  I,  is  continued  until  eventually  each  node  in  A  has 


been  tested  to  determine  if  it  will  provide  a  path  linkage 
to  the  starting  node  from  nodes  in  B,  C,  E  and  H. 

Because  of  the  fact  that  each  border  node  not  reached 
in  the  tree  table  of  step  4  was  given  an  arbitrarily  large 
label,  a  node  not  yet  reached  can  never  provide  a  path 
linkage  for  another  node.   If  the  large  label  was  10,000, 
then  the  comparison  made  in  step  5  will  be: 

w.,  >  10,000  +  w.  . 

which  will  never  allow  a  label  change  on  the  node  j . 

For  this  reason,  the  distance  array  sequences  shown 
in  Table  2  need  not  be  strictly  adhered  to.   The  sequences 
shown  merely  provide  that  a  node  be  reached  before  it  is 
tested  and  that  nodes  are  tested  in  some  logical  order  so 
that  the  number  of  times  step  5  must  be  repeated  is  kept  to 
a  minimum.   However,  any  sequence  used  will  eventually  pro- 
vide the  same  final  result. 

The  reasons  that  step  5  must  be  repeated  until  no  label 
changes  are  made  can  be  visualized  in  Figure  9.   Here  a 
portion  of  the  network  at  some  stage  of  solution  is  shown. 
The  path  from  nodes  in  H,  F,  and  D  has  been  established.   If 
the  gf  distance  array  is  now  introduced  and  the  label  on  a 
node  in  F  changes  so  that  the  path  now  proceeds  from  the  node 
in  F  to  one  in  G,  the  label  on  the  node  in  H  is  in  error 
though  the  path  node  sequence  may  still  be  correct.   The  only 
way  to  correct  the  label  is  to  again  introduce  the  fh  dis- 
tance array  by  repeating  step  5. 


82 


The  border  node  within  border  node  distance  arrays 

such  as  ffk  are  last  in  the  sequences  of  Table  2  for  the 

reason  illustrated  in  Figure  10.   If  w^^,  is  introduced  first 

^  ffk 

and  then  w^^,   ,  the  path  from  points  1  to  3  will  be  as  shown 
by  the  solid  line.   However,  if  w  ,  is  introduced  first,  the 
path  will  proceed  from  points  1  to  3  as  shown  by  the  dotted 
line.   The  paths  in  both  cases  are  identical  but  in  the 
latter,  the  path  trace  procedure  is  somewhat  simplified.   For 
similar  reasons,  it  may  be  desirable  to  introduce  the  border 
node  to  node  distance  arrays  of  step  8  prior  to  the  border 
node  within  border  node  distance  arrays.   In  some  cases, 
however,  this  may  prove  impractical  due  to  computer  storage 
requirements . 

The  path  trace  procedure  required  here  is  somewhat  more 
complicated  than  that  described  for  a  two  part  network.   In 
that  case,  the  presence  or  absence  of  intermediate  nodes  be- 
tween two  border  nodes  indicated  on  which  side  of  the  border 
the  path  lay.   In  this  case,  that  procedure  will  work  only 
when  the  border  nodes  are  part  of  a  border  adjacent  to  the 
area  in  which  the  starting  node  is  included. 

If  the  path  trace  shows  two  consecutive  nodes  from  a 
another  border  set,  say  H,  then  the  side  of  the  border  on 
which  the  path  lies  is  determined  by  comparing  the  difference 
in  labels  of  the  two  nodes  in  the  tree  table  with  their 
entries  in  the  hhk  and  hhi  distance  arrays.   The  smaller  of 
these  entries  must  equal  the  difference  in  labels.   If 
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w,  ,  ,  is  the  smaller,  then  the  portion  of  the  path  in  ques- 
tion utilizes  nodes  in  K. 

The  computer  solution  of  the  algorithm  requires  a 
great  deal  of  temporary  tape  storage  to  get  from  one  step 
to  the  next.   The  distance  arrays  in  certain  instances  may 
consume  more  storage  area  than  that  required  for  solving 
the  entire  network  as  a  single  network. 

However,  there  are  several  reasons  which  make  the 
algorithm  useful  anc"  practical: 

1.  The  distance  arrays  may  be  read  individually  into 
the  same  storage  area  through  the  use  of  equivalence  state- 
ments . 

2.  The  distance  arrays  need  be  calculated  only  once 
and  then  used  in  the  calculation  of  all  required  tree 
tables . 

3.  Because  the  time  required  for  the  solution  of  the 
basic  minimxim  path  algorithm  increases  as  a  function  of  the 
number  of  nodes  in  the  network,  an  actual  time  savings 
develops  provided  that  the  number  of  tree  tables  required 
is  not  too  small. 

The  relationship  between  the  tim.e  required  to  compute 
a  tree  table  and  the  number  of  nodes  is  shown  in  Figure  11. 
If  t  is  the  time  in  seconds  and  N  is  the  number  of  nodes  in 
the  network,  the  relationship  is  given  by: 

t  -  0.55(N/100)^ 
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Based  on  this  relationship,  the  time  required  to  solve  the 
entire  highway  network  of  1837  nodes  (for  the  Indiana  high- 
way problem,  see  Table  3)  would  be  about  187  seconds.  The 
completion  of  the  tree  table  upon  the  incorporation  of  the 
centroid  table  would  require  some  additional  time  so  that 
an  estimate  for  the  complete  tree  table  would  be  about  190 
seconds . 

The  time  required  for  the  complete  tree  table  using 
the  decomposition  algorithm  was  about  87  seconds  for  Indiana 
based  centroids.   This  represented  a  time  savings  of  approx- 
imately 103  seconds  for  each  tree  table. 

A  total  of  about  4410  seconds  was  required  to  generate 
the  distance  arrays  necessary  for  solution  with  the  decom- 
position algorithm.   The  number  of  tree  tables  that  must  be 
solved  to  break  even  with  respect  to  time  is  4410/103  or 
about  43.   A  total  of  451  tree  tables,  one  for  each  Indiana 
based  centroid,  were  calculated  in  the  analysis  of  travel 
desire  factors. 
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DETERMINATION  OF  INTERCITY  TRAVEL  DESIRE  FACTORS 

Introduction 

The  procedure  used  in  determining  the  Intercity  Travel 
Desire  Factors  (ITDF)  consisted  first  of  finding  the  minimum 
path  route  betv;een  a  pair  of  cities.   The  ITDF  between  the 
pair  of  cities  was  computed  and  assigned  to  each  link  making 
up  the  path  provided  that  the  link  represented  some  portion 
of  a  highway  within  Indiana.   As  each  city  pair  was  investi- 
gated, a  cumulative  total  of  the  ITDF ' s  assigned  to  each 
link  was  maintained. 

This  analysis  was  performed  in  three  phases  with  the 
cumulative  total  of  ITDF  link  assignments  kept  separately 
for  each  phase.   The  phases  were: 

I.   All  city  pairs  formed  when  both  cities  were  more 
than  100  miles  from  Indiana's  border. 
II.   All  city  pairs  formed  with  at  least  one  city  in 
Indiana . 
III.   All  city  pairs  formed  with  at  least  one  city  within 
100  miles  of  Indiana's  border  but  neither  within 
Indiana . 
One  of  the  major  assumptions  made  in  this  analysis 
was  that  the  entire  41,000  miles  of  the  Interstate  System  was 
complete  and  open  to  traffic.   The  reason  behind  this 
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assumption  was  that  traffic  patterns  as  they  now  exist  will 
undergo  some  change  upon  completion  of  the  Interstate  system 
and  a  reclassified  state  highway  system  should  reflect 
these  anticipated  changes.   It  was  also  thought  that  some 
insight  might  be  gained  to  guide  the  selection  of  freeway 
standard  highway  sections  to  supplement  the  Interstate 
system. 

Network  Description 

A  total  of  five  partitioned  networks  plus  a  sixth 
network  (the  Interstate  system  outside  of  the  partitioned 
area)  were  coded  such  that  the  entire  network  consisted  of 
4582  nodes  of  which  2745  were  centroids  and  5345  links  of 
which  3229  were  centroid  links. 

The  node  and  link  sizes  of  the  six  networks  are  listed 
in  Table  3.   Each  of  the  networks  was  coded  in  a  slightly 
different  manner  as  explained  below. 

Interstate  Network 

The  first  network  coded  was  the  Interstate  network 
which  was  the  only  one  required  for  Phase  I .   The  assumption 
was  made  that  all  traffic  between  two  cities  where  both  were 
located  more  than  100  miles  from  Indiana's  border,  would 
utilize  only  Interstate  routes  in  crossing  Indiana. 

A  total  of  301  nodes  were  designated  at  intersections 
and  450  individual  links  described.   The  link  lengths  were 
determined  by  use  of  a  mechanical  map  measure  and  appropri- 
ate  scale  based  on  the  state  maps  of  the  Rand  McNally  Road 
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Table  3.   Size  of  Networks, 


Highway 

System 

Cities 

Area 

Links 

Nodes 

L; 

inks 

Centroids 

Interstate 

450 

301 

1 

,639 

1 

,639 

Indiana 

1,809 

1, 

,020 

691 

451 

Michigan 

186 

108 

178 

125 

Ohio 

320 

177 

244 

167 

Kentucky 

178 

101 

147 

114 

Illinois 

259 

163 

5  ^f 

330 

249 

G  CS 

Total 

3,116 

1, 

,837 

3 

,229 

2 

,745 

Atlas  (48).   The  actual  scaled  distance  was  used  as  the  link 
length  for  all  Interstate  links. 

Each  link  required  two  description  cards,  one  for  the 
forward  link  and  another  for  the  back  link.  The  cards  were 
coded  as  shown  below: 

Node  A     Node  B       Link  Length 
Node  B     Node  A       Link  Length 
The  numbering  of  the  nodes  extended  consecutively  from  1500 
to  1800. 

A  total  of  1639  cities,  each  with  a  population  of  5000 
or  more,  were  coded  as  centroids.   For  convenience  in  per- 
forming the  calculations  of  Phase  I,  the  centroid  links 
providing  the  connection  of  the  centroid  to  the  Interstate 
system  were  coded  differently  than  the  highway  links. 

The  Interstate  node  most  likely  to  serve  as  the  entry 
point  to  the  Interstate  system  was  determined  for  each  cen- 
troid.  The  distance  from  the  centroid  to  the  entry  node 
was  obtained  as  described  above  except  that  non-Interstate 
route  miles  were  first  multiplied  by  1.33  prior  to  the  addi- 
tion of  any  Interstate  route  miles.   The  conversion  factor 
of  1.33  was  used  to  change  miles  to  minutes  assuming  a  speed 
of  45  MPH  was  maintained  on  non-Interstate  routes  and  60 
MPH  on  Interstate  routes.   Thus,  in  effect,  all  distances 
used  in  all  of  the  networks  were  a  time  measure  in  minutes 
based  on  a  speed  index  of  one  m.ile  per  minute. 

The  description  of  each  centroid  link  included  three 
pieces  of  data:   The  most  likely  Interstate  entry  node,  the 
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link  length,  and  the  population  of  the  centroid.   All  popu- 
lation data  were  obtained  from  the  1960  U.  S.  Census  of 
Population  (49) . 

Indiana  Network 
A  total  of  1020  highway  nodes  and  1809  links  were 
needed  to  code  the  highway  system  of  Indiana.   The  Interstate 
highways  in  Indiana  were  coded  first  using  node  numbers 
from  2001  to  2199.   Link  lengths  were  scaled  directly  from 
the  1966  Indiana  State  Highway  Map  as  furnished  by  the 
Indiana  State  Highway  Commission. 

All  U.  S.  numbered  highways  in  Indiana  were  coded  using 
node  numbers  2201  to  2630.   Link  lengths  were  obtained  from 
the  Sufficiency  Rating  Report  for  Indiana  State  Highways 
(50).   Prior  to  coding,  the  link  lengths  were  converted  to 
a  time  measure  using  the  following  speeds  and  conversion 
factors  for  the  various  highway  typees: 

4-lane  rural         50  MPH       1.20 
2-lane  rural  4  5  MPH        1.3  3 

4-lane  urban  30  MPH        2.00 

2-lane  urban  15  MPH        4.00 

The  node  locations  coded  up  to  this  point  are  illus- 
trated in  Figure  12. 

Next,  all  Indiana  numbered  highways  were  coded  using 
node  numbers  3001  to  3337.   The  link  lengths  were  treated 
just  as  for  U.  S.  numbered  highways.   Detailed  illustrations 
of  the  complete  Indiana  highway  network  are  shown  in  Figures 
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13  through  20.   The  link  table  is  included  in  Appendix  A, 
Table  Al. 

A  total  of  451  cities  and  towns  in  Indiana  were  coded 
as  centroids  using  centroid  numbers  656  to  1106.   The  break- 
down of  these  cities  by  population  was: 

over  5000  70  cities 

1000-5000  175  cities 

100-1000  206  cities 

Centroid  links  were  coded  similarly  to  the  highway 
links  with  a  forward  link  specifying  the  centroid  number, 
the  highway  node  number,  and  link  length  while  the  back 
link  listed  the  highway  node  number,  centroid  number,  and 
link  length.   A  centroid  located  in  juxtaposition  with  a 
node  was  coded  only  to  that  node  with  the  link  length  taken 
as  a  function  of  population  as  follows: 

over  5000  1.0  miles 

1000-5000  0.5  miles 

100-1000  0.2  miles 

A  centroid  located  along  a  link  between  two  nodes  was 
coded  to  both  nodes  such  that  the  sum  of  the  two  centroid 
link  lengths  equaled  the  highway  link  length. 

The  populations  of  two  cities  within  five  miles  of 
each  other  were  summed  and  treated  as  a  single  city,  the 
location  of  which  was  selected  at  some  intermediate  point. 
The  centroid  table  as  used  in  this  study  is  included 
in  Appendix  A,  Table  A2 . 
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FIGURE    13.       INDIANA    HIGHWAY   NETWORK  - 

SECTION   901-  NORTHWEST   CORNER 
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FIGURE   15.       INDIANA    HIGHWAY    NETWORK    - 

SECTION    903  -   NORTHEAST    CENTRAL 
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FIGURE  17.       INDIANA    HIGHWAY   NETWORK    - 

SECTION  905  -  SOUTHEAST  CORNER 
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FIGURE  18.        INDIANA    HIGHWAY    NETWORK    - 

SECTION    906  -   SOUTHWEST    CENTRAL 
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FIGURE    20.         INDIANA    HIGHWAY     NETWORK   - 
SECTION    908    -    INDIANAPOLIS 
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Michigan  Network 

The  Michigan,  Ohio,  Kentucky,  and  Illinois  highway 
networks  were  each  coded  similarly.   The  networks  extended 
a  distance  of  100  miles  from  Indiana's  borders  so  that  some 
portions  of  Tennessee,  Missouri,  and  Wisconsin  were  also 
included. 

All  cities  and  towns  of  over  1000  population  in  this 
delimited  area  were  coded  as  centroids.   Several  cities  in 
Ohio  and  Illinois  outside  of  the  delimited  area  were  in- 
cluded here  because  there  was  no  most  likely  entry  node 
available  in  the  Interstate  network. 

The  node  and  centroid  number  sequences  as  well  as  a 
population  breakdown  of  the  centroids  are  listed  in  Table  4 
for  each  border  state. 

The  highv/ay  and  centroid  link  lengths  were  determined 
exactly  as  for  the  Indiana  network  except  that  distances 
were  scaled  from  appropriate  state  highway  maps. 

Decomposed  Network 

The  highway  network  used  in  this  study  was  decomposed 
into  five  parts,  which  represented  Indiana  and  portions  of 
Michigan,  Ohio,  Kentucky,  and  Illinois  as  shown  in  Figure 
21.   A  sixth  network  was  that  portion  of  the  Interstate  net- 
work not  included  in  the  other  five  parts. 

Several  assumptions  and  simplifications  were  made  in 
the  calculation  of  minimum  paths  using  the  decomposition 
algorithm. 
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Table  4.   Node  Numbers  in  Border  States. 


r,^    ^  -KT  J   M   1-  Centroid 

State        Node  Numbers 


Numbers   Over  5,000   1,000-5,000 


Michigan  4,001-4,092  1-125  39  86 

Ohio  4,093-4,219  126-292  65  102 

Kentucky  4,220-4,303  293-406  34  80 

Illinois  4,304-4,409  407-655  98  151 
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FIGURE  21.     DELIMITED  AREA  OF  THE  PARTITIONED  NETWORKS. 
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1.  Travel  between  two  cities  within  a  state  were 
assumed  to  utilize  only  highways  within  that  state.   This 
appears  to  be  reasonable  because  the  minimum  path  would 
probably  not  cross  the  state  boundary  unless  the  cities 
were  close  to  the  border,  the  border  itself  was  not  very 
straight,  or  a  high  type  facility  was  available  close  to  the 
border  in  the  adjacent  state.   The  borders  as  shown  in 
Figure  21  are  relatively  straight,  but  the  east-west  Inter- 
state highways  closely  paralleling  the  borders  of  Indiana 

at  the  north  and  south  ends  of  the  State  may  result  in 
border  crossings  for  travel  between  cities  in  both  Michigan 
and  Kentucky.   These  interactions  were  judged  to  be  not 
significant,  however. 

2.  The  minimum  paths  from  one  centroid  to  all  other 
centroids  within  the  same  state  were  not  required  for  any 
state  except  Indiana.   This  greatly  simplified  the  perform- 
ance of  the  decom.position  algorithm  because  it  virtually 
eliminated  step  4  (use  of  the  link  table  to  build  minimum 
path  trees  to  all  other  nodes  from  a  source  node)  for  all 
states  but  Indiana,  requiring  only  the  border  node  to  cen- 
troid distance  array  to  reach  each  border  node.   This  assump- 
tion is  justified  for  two  reasons:   first,  travel  factors 
for  links  in  Indiana  only  were  required,  and  secondly,  such 
travel  factors  would  not  be  affected  if  the  first  assumption 
is  valid. 

3.  Travel  between  two  cities,  one  of  which  was  in 
Indiana,  was  assumed  to  cross  Indiana's  borders  only  once, 
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i.e.,  there  would  be  no  crisscrossing  of  a  border.   The 
same  reasoning  is  applied  here  as  for  the  first  assumption. 

These  assumptions  were  not  made  primarily  to  reduce 
the  magnitude  of  the  computations  but  were  made  necessary 
due  to  physical  limitations  involving  the  computer.   The 
principal  limitation  was  the  lack  of  sufficient  tapes  to 
maintain  all  the  necessary  intermediate  data.   Another  con- 
cerned the  conserving  of  time.   At  any  rate,  any  errors 
introduced  by  these  assumptions  are  considered  negligible. 

This  judgment  was  reinforced  by  subsequent  analysis. 
When  the  decomposition  algorithm  v;as  applied  to  determine 
paths  between  cities  in  different  border  states,  assumption 
3  was  not  made  because  it  was  feasible  to  introduce  the 
distance  arrays  required  to  check  for  the  occurrence  of 
crisscrossing.   It  occurred  in  only  several  instances  and 
in  every  case  was  characterized  by  the  fact  that  one  of  the 
cities  was  very  close  to  the  border  and  the  crossings  were 
made  to  utilize  a  section  of  Interstate  highway. 

Additional  comments  and  discussion  on  the  use  of  the 
decomposed  network  and  algorithm  will  be  made  in  the  follow- 
ing sections  on  travel  desire. 

Travel  Desire-Phase  I 
Phase  I  of  the  determination  of  Intercity  Travel 
Desire  Factors  included  all  city  pairs  formed  where  both 
cities  were  located  more  than  100  miles  from  Indiana's 
border.   The  Interstate  network  consisting  of  301  nodes  and 
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450  links  v/as  used  in  this  phase.   All  city  pairs  formed 
from  a  total  of  1639  centroids  were  tested  to  determine  if 
the  minimum  path  between  them  utilized  any  Indiana  Inter- 
state routes.   If  so,  the  factor  was  calculated  and  assigned 
to  those  routes  included  in  the  minimum  path, 

A  modified  version  of  Small  Trees,  a  computer  program 
listed  in  Appendix  Bl,  was  used  to  determine  all  required 
minimum  paths.   A  total  of  10  Interstate  system  nodes  and 
23  links  were  included  in  Indiana.   The  first  step  was  the 
formation  of  a  distance  array  FINDR  which  included  the  mini- 
mum path  distances  from  each  of  the  10  Interstate  nodes 
within  Indiana  to  each  of  the  301  nodes  in  the  netv/ork. 

The  remainder  of  the  procedure  was  confined  to  those 
Interstate  nodes  to  which  centroids  were  linked  beginning 
with  the  lowest  Interstate  node  num.ber  proceeding  to  the 
highest. 

The  tree  table  was  first  calculated  for  each  node 
and  the  paths  to  it  traced  from  all  higher  numbered  nodes. 
If  the  path  trace  contained  at  least  one  of  the  10  Inter- 
state nodes  located  in  Indiana,  the  path  trace  was  printed. 
The  route  through  Indiana  was  coded  by  a  sequence  of  numbers 
from  1  to  23,  each  of  which  designated  one  of  the  23  Inter- 
state links  in  Indiana.   Each  individual  route  was  assigned 
a  code  numiber  from  1  to  2  9  which  represented  each  of  2  9 
different  routes  used  in  passing  through  Indiana. 

Thus,  for  each  different  Indiana  route  employed  in 
reaching  the  source  node,  a  set  of  at  least  two  data  cards 
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was  generated.   The  first  contained  the  source  node  number, 
the  route  number,  and  the  total  of  higher  numbered  nodes 
which  used  the  designated  route  to  reach  the  source  node. 
The  second  card  (and  others  as  required)  listed  the  node 
numbers  corresponding  to  the  total  given  on  the  first  card. 

Each  of  the  29  routes  was  coded  similarly.   The  first 
card  listed  the  route  number,  the  key  node,  and  total  links 
on  the  route  while  the  second  listed  the  link  numbers.   The 
key  node  was  any  one  of  the  ten  Interstate  nodes  within 
Indiana  which  was  included  in  the  route. 

After  each  node  had  been  used  as  a  source  node  and  all 
of  the  node-node  interactions  utilizing  Indiana  routes 
properly  coded,  the  Intercity  Travel  Desire  Factors  were 
computed. 

The  distance  d.  .  v;as  the  sum  of  the  four  distances 
obtained  from  the  input  data.   Two  were  the  centroid  to 
Interstate  entry  node  distances  obtained  from  the  centroid 
table.   The  other  two  distances  were  the  entry  node  to  key 
node  distances  obtained  from  the  FINDR  array.   The  factor 
was  computed  and  then  assigned  to  each  link  on  the  route 
connecting  the  two  controids.   A  cumulative  total  of  these 
factors  was  maintained  for  each  link. 

A  total  of  346,734  centroid  interactions  were  computed 
for  this  phase  of  the  analysis.   The  final  total  of  Inter- 
city Travel  Desire  Factors  for  each  of  the  192  Indiana  high- 
way links  comprising  the  29  Indiana  Interstate  routes  are 
listed  in  Appendix  A,  Table  Al  under  the  heading.  Phase  I. 
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Travel  Desire-Phase  II 
Phase  II  of  the  determination  of  Intercity  Travel 
Desire  Factors  included  all  city  pairs  formed  with  at  least 
one  city  within  Indiana.   The  area  node  and  border  node  set 
code  names  used  in  the  description  of  Phase  II  are  shown 
in  Figure  21.   The  number  of  nodes  in  each  border  set  are 
in  parentheses  alongside  the  border  set  code  name  while  the 
node  sets  used  in  generating  the  necessary  minimum  path  trees 
are  shown  similarly  within  each  area. 

Figure  21  also  shov7s  10  Interstate  network  nodes 
selected  in  such  a  way  that  the  delimited  area  is  bounded 
by  them.   Thus,  the  minimum  path  from  any  Interstate  node 
outside  the  delimited  area  to  any  node  within  the  delimited 
area  must  pass  through  at  least  one  of  these  ten  nodes, 
hereafter  referred  to  as  boundary  nodes.   All  Interstate 
nodes  within  the  boundary  nodes  are  included  in  the  appropri- 
ate area  node  sets  of  the  partitioned  netv/orks . 

The  general  procedure  used  to  compute  the  Intercity 
Travel  Desire  Factors  was  as  follows: 

1.  Select  a  source  centroid  in  Indiana  and  using  the 
JBCDE  link  table  (listed  in  Appendix  A,  Table  Al) ,  calculate 
the  minimum  path  tree  table  to  all  nodes  and  centroids 
within  Indiana  employing  the  computer  program.  Small  Trees. 

2,  Use  this  tree  table  and  the  required  distance 
arrays  to  generate  an  expanded  minimum  path  tree  table  for 
all  other  centroids  in  the  network  employing  the  computer 
program,  Expanded  Trees. 
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3.  With  the  expanded  tree  table  and  the  centroid 
population  arrays,  compute  the  travel  desire  factor  between 
a  centroid  and  the  source  centroid,  trace  the  minimum  path 
node  sequence,  and  assign  the  factor  to  each  Indiana  link  on 
the  path,  employing  the  computer  program.  Trades. 

4.  Update  the  cumulative  sum  of  travel  desire  fac- 
tors calculated  for  previous  source  centroids  by  adding 
those  calculated  for  the  current  source  centroid. 

Phase  II  was  terminated  when  each  Indiana  centroid  had 
been  used  as  a  source  centroid. 

Distance  Arrays 
The  distance  arrays  required  for  use  of  the  decomposi- 
tion algorithm  in  Expanded  Trees  were  generated  first.   These 
arrays  as  described  below  were  generated  by  a  modified  ver- 
sion of  Small  Trees  using  the  appropriate  area  link  and  cen- 
troid tables.   Only  the  border  nodes  were  used  as  source 
centroids  in  this  step. 

The  required  distance  arrayss  were: 

DISTMC-   The  distances  from  all  M  (Michigan)  centroids 

to  the  border  nodes  in  sets  B  and  H. 
DISTLC-   The  distances  from  all  L  centroids  to  the 

border  nodes  in  sets  H,  C,  and  F. 
DISTKC-   The  distances  fromi  all  K  centroids  to  the 

border  nodes  in  sets  F,  D,  and  G. 
DISTIC-   The  distance  from  all  I  centroids  to  the 
border  nodes  in  sets  G  and  E. 
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DMBLLL-   The  distances  from  each  node  in  H  to  all 

nodes  in  H,  C,  and  F  by  v/ay  of  nodes  in  L.   The 
array  name,  DMBLLL,  may  be  interpreted  as  follows: 
The  distance  (D)  from  all  Michigan  nodes  bordering 
Ohio  (MBL)  to  all  nodes  bordering  Ohio  (L)  by  way 
of  nodes  in  Ohio  (L) .   The  companion  array,  DLBMMM, 
was  generated  but  discarded  because  it  was  suspected 
and  easily  verified  that  all  minimum  paths  from 
nodes  in  H  to  nodes  in  B  were  by  way  of  nodes  in 
areas  L  and  J  due  to  the  presence  of  the  Ohio  and 
Indiana  Turnpikes. 

DKBLLL-   The  distances  from  each  node  in  F  to  all 
nodes  in  F,  C,  and  H  by  way  of  nodes  in  L. 

DLBKKK-   The  distances  from  each  node  in  F  to  all 
nodes  in  F,  D,  and  G  by  way  of  nodes  in  K. 

DIBKKK-   The  distances  from  each  node  in  G  to  all 
nodes  in  G,  D,  and  F  by  way  of  nodes  in  K. 

DKBIII-  the  distance  from  each  node  in  G  to  all  nodes  in 
G  and  E  by  way  of  nodes  in  I . 

DISTMN-   The  distances  from  each  Interstate  node  (N) 
included  in  the  Michigan  area  node  set,  MBH  ,  in- 
cluding boundary  nodes,  if  any,  to  the  border 
nodes  in  B  and  H. 

DISTLN,  DISTKN,  DISTIN-   Each  of  these  is  interpreted 
similarly  to  DISTMN  but  for  each  of  the  other 
three  border  states. 
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NNATND,  NACEND,  DISTNC-   This  is  a  link  table  generated 
by  using  each  of  the  boundary  nodes  as  a  source 
centroid  with  the  Interstate  network.   NNATND  is 
one  of  the  boundary  nodes.   NACEND  is  a  centroid 
linked,  Interstate  node  connected  to  a  boundary 
node  at  a  distance,  DISTNC.   Note  that  this  is 
not  a  distance  array  in  the  same  sense  as  those 
described  above. 

Small  Trees 

Small  Trees  is  a  computer  program  listed  in  Appendix 
Bl .   It  was  used  to  solve  the  basic  minimum  path  algorithm 
with  the  Indiana  network  link  table  for  a  source  centroid 
within  Indiana.   The  details  concerning  the  basic  minimum 
path  algorithm  were  given  in  a  previous  chapter. 

The  mechanics  of  Small  Trees  are  outlined  briefly  in 
Appendix  Bl. 

Expanded  Trees 

Expanded  Trees  is  a  computer  program  listed  in  Appen- 
dix B2.   It  was  used  to  extend  the  tree  table  calculated 
by  Small  Trees  to  include  the  minimum  path  distances  to  all 
other  centroids  in  the  network.   The  details  concerning  the 
decomposition  algorithm  for  large  networks  were  given  in 
the  previous  chapter. 

The  mechanics  of  Expanded  Trees  are  outlined  briefly 
in  Appendix  B2. 
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Trades 

Trades  is  a  computer  program  listed  in  Appendix  B3 . 
It  was  used  to  calculate  the  Intercity  Travel  Desire  Factors 
between  the  source  centroid  and  all  other  centroids  and 
assign  the  factors  to  each  link  on  the  routes  connecting 
each  centroid  pair. 

Trades  as  originally  programmed  exceeded  available 
storage  by  about  3000  locations.   The  Interstate  centroid 
array  required  4917  locations  for  the  three  data  elements 
of  each  of  1639  centroids.   Because  many  centroids  utilized 
the  same  Interstate  entry  node,  it  was  reasoned  that  the 
sum  of  centroid  populations  and  an  average  distance  to  the 
Interstate  entry  node  could  be  substituted  for  the  group  of 
centroids,  thus  reducing  the  storage  requirements.   This 
procedure  was  subject  only  to  the  provision  that  any  error 
introduced  was  negligible. 

The  form  of  the  travel  desire  factor  used  in  this 
study  was  the  product  of  the  square  roots  of  populations  of 
two  centroids  divided  by  the  square  of  the  distance  between 
them : 

1/2     2 
ITDF.  .  =  (P  .  P  .)  '     /d,  . 
ij      ID       IJ 

The  desired  simplification  would  be  of  the  form: 

^  Y      1/2         -  2 

^  ITDF.  .  -  (P.   ''       P.)^  /(d.   +  d)'' 

11      1       D        m 

J  J 

where  n  is  the  Interstate  entry  node  and  i,  the  source 
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centroid.   It  was  found  that  the  accuracy  of  the  simplifica- 
tion depended  on  three  things : 

1.  The  average  distance  must  be  based  on  the  square 
roots  of  population  rather  than  on  population, 

2.  The  less  variable  the  d  .  distances  were,  the  more 

nj 

accurate  the  simplification,  and 

3.  The  larger  the  value  of  the  constant  distance  d.  ; 
the  more  accurate  the  simplification. 

Therefore,  the  following  procedure  was  used  in  provid- 
ing the  simplification.   The  average  distance  was  computed 
as  : 


d,  v^P,+d_vP„  + 
j     nl     nl     n2     n2 

^  "  ^nj 


and  the  travel  desire  as 


I    ITDF.  .  -  (P.)^/^  [  I     (P  .)^/^/(d.   +  d)^] 
V      13      1        V    nj     'in 


In  order  to  reduce  the  effects  of  the  variability  in  the 
d  . ,  d  was  calculated  for  each  of  five  distance  ranges  for 
each  Interstate  entry  node.   The  ranges  v/ere  0  to  100  miles, 
100  to  300  miles,  300  to  600  miles,  600  to  1000  miles,  and 
over  1000  miles. 

Because  there  were  only  86  Interstate  entry  nodes  and 
only  a  few  \vith  distances  in  all  five  ranges,  a  total  of  181 
Interstate  centroid  cards  requiring  543  storage  locations  were 
needed. 
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This  procedure  was  verified  with  the  actual  data  used 
in  the  study  and  the  error  introduced  was  insignif iciant . 

The  mechanics  of  Trades  are  outlined  briefly  in 
Appendix  B3. 

Update  Trades 

Because  of  time  limitations  imposed  on  the  length  of 
jobs  run  on  the  computer,  the  three  programs  --  Small  Trees, 
Expanded  Trees,  and  Trades  --  could  be  run  in  succession  for 
only  about  35  to  40  source  centroids  at  a  time.   Therefore, 
after  each  run  it  was  required  that  the  output  of  Trades 
be  added  to  that  of  all  previous  Trades  and  the  result 
stored  on  another  tape. 

A  total  of  209,393  centroid  interactions  were  computed 
for  this  phase  of  the  analysis.   The  final  totals  of  Inter- 
city Travel  Desire  Factors  for  each  of  the  2510  Indiana  cen- 
troid and  highway  links  are  listed  in  Appendix  A,  Tables 
Al  and  A2  under  the  heading.  Phase  II. 

Travel  Desire-Phase  III 
Phase  III  of  the  determination  of  Intercity  Travel 
Desire  Factors  included  all  city  pairs  formed  v/ith  at  least 
one  city  within  100  miles  of  Indiana's  borders  but  neither 
within  Indiana.   The  centroids  used  in  this  portion  of  the 
analysis  were  restricted  to  those  with  a  population  of  5000 
or  more.   The  area  node  and  border  node  set  code  names  used 
in  the  description  of  Phase  III  are  shown  in  Figure  21. 
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The  general  procedure  used  to  compute  the  Intercity 
Travel  Desire  Factors  was  as  follows: 

1.  Select  a  source  centroid  from  one  of  the  border 
states  and,  using  the  required  distance  arrays,  calculate 
the  minimum  path  tree  table  to  all  nodes  and  centroids  em- 
ploying the  computer  program.  Border  Trees. 

2.  Determine  the  portion  of  the  path  consisting  of 
Indiana  highway  links  and  the  centroid  pairs  connected  by 
each  different  path. 

3.  With  the  data  of  step  2  and  the  centroid  popula- 
tion arrays,  compute  the  travel  desire  factor  between  a  cen- 
troid and  the  source  centroid  and  assign  the  factor  to  each 
Indiana  highway  link  on  the  path  employing  the  computer 
program  Tradet. 

Phase  III  was  terminated  when  each  border  state  cen- 
troid had  been  used  as  a  source  centroid. 

Distance  Arrays 

The  distance  arrays  required  for  use  of  the  decomposi- 
tion algorithm  in  Border  Trees  were  generated  first.   These 
arrays  as  decribed  below  were  generated  by  a  modified  ver- 
sion of  Small  Trees  using  the  appropriate  area  link  and 
centroid  tables.   Only  border  nodes  were  used  as  source 
centroids  in  the  step  and  only  those  centroids  with  a  popu- 
lation of  over  5000  were  considered. 

Several  distance  arrays  required  for  Border  Trees  were 
also  required  for  Expanded  Trees.   These  were  defined  in  the 
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previous  section.  Travel  Desire-Phase  II  and  are  simply 
listed  here:   DISTMC,  DISTLC,  DISTKC,  DISTIC,  DMBLLL,  DKBLLL, 
DLBKKK,  DIBKKK,  DKBIII,  DISTMN,  DISTLN,  DISTKN,  DISTIN,  and  the 
link  table  providing  connections  froin  the  boundary  nodes  to  cen- 
troid  linked.  Interstate  nodes,  NNATND,  NACEND,  and  DISTNC. 
Other  distance  arrays  required  were: 

DISTML-   The  distance  from  each  node  in  B  to  all  nodes 
in  C  by  way  of  nodes  in  J,  i.e.,  the  distance 
from  the  Michigan-Indiana  border  nodes  to  the 
Ohio-Indiana  border  nodes  as  calculated  using  the 
JBCDE  link  table. 
DISTiMK-   The  distance  from  each  node  in  B  to  all  nodes 

in  D  by  way  of  nodes  in  J. 
DISTMI-   The  distance  from  each  node  in  B  to  all  nodes 

in  E  by  way  of  nodes  in  J. 
DISTLK-   The  distance  from  each  node  in  C  to  all  nodes 

in  D  by  way  of  nodes  in  J. 
DISTLI-   The  distance  from  each  node  in  C  to  all  nodes 

in  E  by  way  of  nodes  in  J. 
DISTPCI-   The  distance  from  all  nodes  in  D  to  each  node 

in  E  by  way  of  nodes  in  J. 
DMBJMJ-   The  distance  from  each  node  in  B  to  all  nodes 
in  B  by  way  of  nodes  in  either  M  or  J ,  i.e.,  the 
distance  from  each  Michigan-Indiana  border  node  to 
all  other  Tlichigan-Indiana  border  nodes  where  the 
entry  in  the  distance  array  was  the  minimum  of  two 
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possible  paths,  the  first  entirely  in  Michigan 

and  the  second  entirely  in  Indiana.   In  effect, 

two  distance  arrays  were  combined  into  one  by  this 

procedure . 
DLBJLJ-  The  distance  from  each  node  in  C  to  all  nodes 

in  C  by  way  of  nodes  in  either  L  or  J . 
DKBJKJ-  The  distance  from  each  node  in  D  to  all  nodes 

in  D  by  way  of  nodes  in  either  K  or  J. 
DIBJIJ-  The  distance  from  each  node  in  E  to  all  nodes 

in  E  by  way  of  nodes  in  either  I  or  J. 

Border  Trees 

Border  Trees  is  a  computer  program  listed  in  Appendix 
B4.   It  was  used  to  solve  the  decomposition  algorithm  for 
determ.ining  minimum  paths  in  large  networks. 

The  mechanics  of  Border  Trees  are  outlined  briefly 
in  Appendix  B4. 

Tradet 

Tradet  is  a  modified  version  of  Trades,  a  computer 
program  listed  in  Appendix  B3.   It  was  used  to  calculate  the 
Intercity  Travel  Desire  Factors  between  each  centroid  pair 
and  assign  the  factors  to  each  link  on  the  Indiana  portion 
of  the  routes  connecting  them. 

The  output  of  Border  Trees  was  coded  in  a  manner  simi- 
lar to  that  used  for  Travel  Desire-Phase  I.  For  each  source 
centroid,  the  distance  array  to  all  other  centroids  was  card 
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punched.   The  Indiana  portion  of  each  route  connecting  a 
centroid  to  the  source  centroid  was  assigned  a  route  number 
and  all  centroids  utilizing  that  route  were  so  coded. 

A  total  of  87  different  routes  utilizing  Indiana  high- 
ways were  needed  to  account  for  all  centroid  interactions 
where  the  source  centroid  was  located  in  Michigan.   For 
source  centroids  in  Ohio,  Kentucky,  and  Illinois,  the  number 
of  different  routes  was  91,  58,  and  46  respectively. 

The  mechanics  of  Tradet  and  Trades  are  similar  but 
with  several  exceptions.   The  path  tracer  and  the  centroid 
pairs  utilizing  each  path  trace  (route)  v;ere  read  in  for 
Tradet  but  had  to  be  determined  from  the  tree  table  in 
Trades.   It  was  possible  to  compute  all  the  travel  desire 
factors  in  one  computer  run  with  Tradet  but  ten  were 
required  to  complete  Trades. 

A  total  of  33,814  centroid  interactions  were  computed 
for  this  phase  of  the  analysis.   The  final  totals  of  Inter- 
city Travel  Desire  Factors  for  each  Indiana  link  involved 
are  listed  in  Appendix  A,  Table  Al  under  the  heading, 
Phase  III. 

Summary 
The  methodology  employed  in  computing  the  Intercity 
Travel  Desire  Factors  representative  of  Indiana  highv/ays  has 
been  presented  in  this  chapter.   The  travel  desire  factors 
are  listed  in  Appendix  A  opposite  the  link  each  represents. 
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In  Phases  I  and  III  of  this  analysis,  the  centroid 
pairs  and  the  route  connecting  each  pair  were  coded  and 
card  punched.   This  was  done  with  the  intention  that  other 
forms  of  a  travel  desire  factor  besides  that  used  could  be 
tried.   However,  this  was  not  feasible  in  Phase  II  because 
virtually  every  centroid  pair  was  connected  by  a  different 
route.   An  analysis  of  other  forms  of  a  travel  desire  factor 
could  have  been  performed  if  sufficient  tapes  v^?ere  available 
to  store  the  expanded  tree  table  for  each  Indiana  source 
centroid. 

The  number  of  centroid  interactions  computed  in  Phases 
II  and  III  as  reported  above  are  in  error  because  of  the 
grouping  of  Interstate  centroids  made  necessary  to  reduce 
computer  storage  requirements  as  described  in  the  subsection. 
Trades.   If  this  is  corrected,  an  estimate  of  the  number  of 
centroid  interactions  computed  in  this  analysis  is: 

Phase  I  346,734 

Phase  II  311,453 

Phase  III  130,000 

for  a  total  of  approximately  790,000  interactions. 
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ANALYSIS  OF  INTERCITY  TRAVEL  DESIRE  FACTORS 

Introduction 

The  purpose  of  this  chapter  is  to  demonstrate  the 
adequacy  of  the  assumed  ITDF  model  to  synthesize  travel. 
Current  and  past  traffic  volume  data  could  not  be  used  direc- 
tly as  the  system  analyzed  was  the  current  Indiana  State 
Highway  System  with  the  Interstate  System  complete.   The 
ITDF  model  used  in  this  study  was  the  product  of  the  square 
roots  of  the  populations  of  a  city  pair  divided  by  the  square 
of  the  minimum  path  distance  between  the  cities. 

The  number  resulting  from  this  calculation  for  a  city 
pair  was  then  assigned  to  each  link  making  up  the  minimum 
path  sequence.   A  cumulative  total  was  kept  for  each  link  so 
that,  after  all  city  pairs  had  been  considered,  each  link 
was  represented  by  an  Intercity  Travel  Desire  Factor,  here- 
after called  the  factor. 

Because  the  factors  were  computed  in  three  Phases,  some 
links  had  more  than  one  factor  associated  with  them.   In 
Phase  I,  the  assumption  was  made  that  all  centroid  interac- 
tions would  be  assigned  to  Interstate  routes  only;  hence, 
only  these  had  a  non-zero  Phase  I  factor.   This  assumption 
was  not  made  in  Phase  III  but,  nevertheless,  the  greater 
proportion  of  the  Phase  III  factors  were  for  Interstate 
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routes.   This  development  indicated  that  the  assumption  of 
only  Interstate  route  assignments  in  Phase  I  appeared 
reasonable.   Furthermore,  more  than  three-fourths  of  the 
Indiana  state  highway  links  had  only  factors  resulting  from 
Phase  II. 

Because  of  the  way  in  which  the  network  was  coded, 
each  Indiana  centroid  had  a  Phase  II  factor  associated  with 
it  which  represented  the  magnitude  of  its  interaction  with 
all  other  centroids,  i.e.,  when  the  centroid  served  as  an 
origin  or  destination  for  interaction  with  other  centroids. 
This  meant  that  the  factors  associated  with  the  links  coming 
into  a  centroid  could  be  separated  into  what  may  be  termed 
a  centroid  attraction  factor  and  a  thru  factor. 

All  of  the  factors  are  listed  in  the  Link  Table  and 
Centroid  Table  of  Appendix  A. 

The  analysis  of  the  factors  was  performed  in  two  parts, 
a  centroid  and  a  link  analysis. 

Centroid  Analysis 
It  has  prev.'.cusly  been  stated  that  the  completion  of 
the  Interstate  highway  system  and  its  subsequent  utilization 
by  traffic  will  result  in  changes  in  the  traffic  patterns 
now  existing  on  the  state  highway  system.   Because  of  the 
fact  that  the  Interstate  system  was  assumed  to  be  completed 
in  the  performance  of  this  study,  the  factors  should  reflect 
travel  patterns  not  as  they  currently  exist  but  as  they  would 
exist  if  the  Interstate  system  were  com.plete  or  as  they  will 
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exist  in  approximately  1972  (the  scheduled  completion  date 
of  the  Interstate  system)  provided  that  growth  is  more  or 
less  uniform.   For  this  reason,  it  was  decided  that  the 
establishment  of  a  relationship  between  link  volumes,  as 
they  currently  exist,  and  link  factors  would  not  provide  an 
adequate  means  of  demonstrating  the  adequacy  of  the  ITDF 
model  to  synthesize  travel  although  a  link  analysis  to  es- 
tablish such  a  relationship  was  necessary  if  the  factors 
were  to  serve  their  function  as  a  highway  classification  tool. 

However,  it  was  reasoned  that  even  though  the  traffic 
patterns  would  be  changed,  the  relative  magnitude  of  traffic 
attracted  by  individual  centroids  would  not  be  affected  by 
assuming  a  completed  Interstate  system.   This  implies  that 
the  traffic  attracted  to  a  centroid  should  be  essentially 
the  same  but  that  it  may  enter  the  centroid  on  different 
links.   Also,  it  must  be  recognized  that  the  amount  of  traf- 
fic which  currently  passes  thru  a  city  might  be  significantly 
altered  if  the  Interstate  system  were  complete. 

The  impact  of  this  latter  effect  on  the  total  traffic 
entering  or  leaving  a  city  decreases  as  city  size  increases 
as  shown  in  Table  5  from  Matson,  Smith,  and  Kurd  (53) . 
Table  5  gives  the  per  cent  of  traffic  entering  a  city  which 
may  be  bypassed  (traffic  not  desiring  to  stop  in  the  city) 
around  the  city  as  a  function  of  city  size.   The  percentage 
of  bypassable  traffic  for  cities  from.  10,000  to  300,000 
population  is  relatively  uniform  appearing  to  average  about 
20  per  cent.   Because  of  the  high  proportion  of  traffic 
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Table  5.   Proportion  of  Approaching  External  Traffic  which 
May  Be  Bypassed  Around  A  City. 


Cities  by 
Population  Size 


Per  Cent 
Bypassable  Traffic 


Less  than  2,500 

2,500  -  10,000 

10,000  -  25,000 

25,000  -  50,000 

50,000  -  100,000 

100,000  -  300,000 

300,000  -  500,000 

500,000  -  1,000,000 


50.7 
43.3 
21.9 
21.0 
16.2 
18.2 
7.2 
4.2 
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having  a  terminus  in  these  cities,  it  appeared  reasonable 
to  expect  that  a  significant  relationship  should  exist  be- 
tween the  sum  of  the  computed  factors  and  the  sum  of  the 
traffic  volumes  for  the  links  entering  these  cities  even  if 
significant  changes  in  thru  traffic  would  occur  with  a  com- 
pleted Interstate  system,  provided  that  the  ITDF  model  was 
adequate.   In  other  words,  a  high  percentage  of  the  variabil- 
ity of  traffic  volume  should  be  explained  by  a  regression  of 
the  sum  of  the  traffic  volumes  for  the  links  serving  the 
centroid  on  the  Phase  II  centroid  attraction  factor. 

In  order  to  test  this  reasoning,  several  regression 
analyses  were  performed  utilizing  the  Stepwise  Regression 
program,  BMD-2R,  as  developed  by  the  Health  Sciences  Comput- 
ing Facility  at  UCLA  (43) . 

The  dependent  variable  Y  was  taken  as  the  sum  of  the 
Annual  Average  Daily  Traffic  (AADT)  estimates  on  all  high- 
ways entering  a  city  as  measured  just  beyond  the  city 
limits.  This  data  was  collected  from  the  1962  Traffic  Map 
as  furnished  by  the  Indiana  State  Highway  Commission  (55). 
The  independent  variables  were  various  measures  of  the  com- 
puted factors  associated  with  the  links  entering  each  cen- 
troid. 

Data  for  a  total  of  390  Indiana  centroids  were  coded 
for  use  in  this  analysis.   These  data  sets  were  divided  into 
three  centroid  population  groups.   The  first  group  consisted 
of  59  centroids,  each  with  a  population  of  over  5000,  the 
second  of  167  centroids,  each  with  a  population  between  1000 
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and  5000,  and  the  third  of  154  centroids,  each  with  a 
population  between  100  and  1000.   A  regression  analysis  was 
performed  on  each  of  these  three  data  sets  as  well  as  on 
the  total  of  390  data  sets. 

The  division  of  data  sets  by  population  was  performed 
to  test  the  reasoning  that  the  adequacy  of  the  ITDF  model 
would  decrease  as  population  decreased  as  measured  by  the 
per  cent  of  variability  of  the  dependent  variable  explained 
by  the  independent  variables.   The  per  cent  of  variability 
explained  is  commonly  referred  to  as  the  coefficient  of 

determination  or  the  square  of  the  correlation  coefficient 

2 
(R  )  in  regression  analysis. 

2 
It  was  also  expected  that  an  increase  in  R   would 

occur  when  the  thru  factors  associated  with  the  links 
serving  a  centroid  were  included  in  the  analysis.   However, 
it  was  also  recognized  that  a  serious  problem  concerning 
the  effects  of  local  traffic  existed  for  the  smaller  cities 
in  addition  to  the  thru  traffic  effects.   The  situation  as 
it  generally  exists  for  small  cities  is  that  traffic  volumes 
on  roads  entering  the  city  are  correspondingly  small.   A 
good  deal  of  this  traffic  may  be  bound  for  the  city  but  be 
local  in  character.   The  am.ount  of  local  traffic  is  a  func- 
tion of  the  population  density  in  the  surrounding  area  and 
the  degree  of  agricultural  and  other  economic  development. 
These  factors  have  wide  variation  throughout  the  State  so 
that  a  similar  variation  should  be  expected  in  the  amount 
of  local  traffic  as  well. 
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2 
Each  of  these  effects  serve  to  force  a  decrease  in  R 

as  city  size  decreases  because  the  amount  of  thru  and  local 
traffic  assumes  a  much  greater  significance  with  respect 
to  the  total  traffic  entering  the  city  than  does  the  inter- 
city traffic  attracted  by  the  city.   This  would  be  expected 
to  occur  even  if  the  ITDF  model  perfectly  explained  inter- 
city travel  desire. 

The  results  of  the  regression  analyses  performed  on 

2 
the  centroid  data  are  given  in  Table  6.   The  values  of  R 

are  shown  to  decrease  from  0.829  to  0.218  to  0.116  for  large, 

medium,  and  small  city  sizes  respectively  when  the  centroid 

attraction  factor  is  the  only  independent  variable.   However, 

2 
the  values  of  R   are  increased  to  0.879,  0.479,  and  0.419 

respectively  when  the  thru  factor  measures  of  Phases  I,  II, 

and  III  also  serve  as  independent  variables.   These  results 

are  in  agreement  with  the  discussion  given  above. 

2 
The  values  of  R   for  the  data  taken  as  a  whole  range 

from  0.816  when  the  centroid  attraction  factor  is  the  only 

independent  variable  to  0.872  when  the  three  thru  factor 

measures  are  also  included. 

Based  on  these  results,  it  appears  reasonable  to  con- 
clude that  the  assumed  ITDF  model  is  an  adequate  measure  of 
intercity  travel  desire. 

A  perplexing  aspect  of  the  models  developed  in  this 
analysis,  however,  is  the  negative  sign  of  some  of  the  coef- 
ficients for  the  Phase  I  and  Phase  III  thru  factors,  i.e., 
X   and  X  .   (See  Table  6)   The  implication  is  that  traffic 


128 


M 

0) 
.H 
XI 

(0 

-^^  ; 
M 
rO 
> 

-P 
C 
Q) 
T3 
C 

0) 


-P 

C  0) 

0)  rH 

G  ITJ 

0)  -H 
O,  M 

Q  > 


•^ 

(J\ 

1   1  CM 

US 

1   1  <N 

iH 

1  1  vo 

r-H 

1    1   r-l 

(N 

n  \0  1J3   03 

I  I     U5  in  i-l  CTl 

II  .... 

I       I     (N  O  -^   O 


o  fo  Ln  •»»  in  -T 


I 


I  I  I  I  in  r~ 
I  I  I  I  •  . 
I     I     I     I   ^  .-I 


CO   O  ■^   <D 


(N  vu  1-t  IX)  ^  in 


I     U3   iH 

I 


r~  r~  n  in 
I     I   <N  n  yj  00 


I     (N  O  ID  o 


M    M    M 

M    M 

M 

.  W   .  W   •  W 

.  W  •  W 

•  w 

l«-l    4-1    14-1 

m   iM 

14-1 

g-i  .in  .  M-l  . 

M-l   .  IM   . 

14-1   . 

0)  X)  <D  Tl  0)  T3 

(1)  T)  0)  'O 

0)  01 

O  -P  0  -P  0  4-1 

O  -P  O  -P 

O  -P 

u       u 

•  w    •  w 

14-1  14-1 


O  +)   O  -P   O  -P 


CJMUWUW       uwucouw       uwuwucn 


I     fM  O 

I 


I    o  o 
I 


I       I     CM   O  -^  O 


00  i-H  r-~  i-H  r^  .H 


O  -P   O  -P   O  -P 
U  M  L)  en  u  w 


Q) 

o  c 

o 

n3  -H 

tr>  u 

•H    U 


M   -P  M 

0)  -P  M 

4-)     (0  O 

C  +J 

0)  TJ  O 

■H  (0 

to    O  "H 

M  u 

c:  -p  3 

H    C  M 

-I  0)  x; 

O  4J 

c 

OHM 

M  H 

w 

E-t 
Q 


flj    0) 
0)    Ul 

x:  x: 


.  u 

w  O 

M  4J 

0  u 

-P  nj 

U  >4-l 

(0 

4-1  3 
M 

3  x: 

M  -P 

x: 

4J  M 

M 

H  H 

Q)  0) 

01  U) 

x:  x: 

CU  04 


4-1   14H   M-l    14-1    IM 

o  o  o  o  o 


tn  u)  M  t/i  Lo 
II   II    II    II 


^  (N  n  ■q. 
?H  X  X  X   X 


129 


decreases  as  the  thru  factors  associated  v/ith  a  centroid 
increases.   This  makes  no  sense  at  all. 

An  explanation  for  this  occurrence  is  an  extension  of 
the  thru  traffic  effect  as  described  previously.   Thru  traf- 
fic makes  up  a  relatively  small  proportion  of  the  external 
traffic  entering  a  city  and  it  may  be  assumed  that  the  type 
of  thru  traffic  supposedly  measured  by  the  Phase  I  and  Phase 
III  factors  constitutes  a  small  proportion  of  the  thru  traf- 
fic.  Thus,  other  traffic  variations  such  as  in  local  traf- 
fic probably  overwhelm  any  variability  in  traffic  volume 
that  may  be  explained  by  the  Phase  I  and  Phase  III  factors. 
This  was  borne  out  by  examination  of  the  standard  errors 
which  revealed  that  many  of  the  coefficients  of  X-.  and  X. 
are  not  significantly  different  from  zero. 

The  high  value  of  some  of  the  Phase  I  and  Phase  III 
factors  leads  to  the  conclusion  that  they  are  overestimated 
relative  to  Phase  II  factors.   In  other  words,  the  ITDF 
model  yields  too  high  an  intercity  travel  desire  for  the 
very  long  interactions  characteristic  of  Phase  I  and  Phase 
III.   Many  of  these  long  trips  are  either  not  made  at  all 
or  are  made  by  another  mode  such  as  plane  or  train.   Thus, 
for  long  distance  trips  of  say  300  to  400  miles  or  more,  the 
exponent  of  distance  in  the  ITDF  model  should  be  some  value 
greater  than  two. 
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Link  Analysis 

A  link  analysis  was  performed  to  develop  a  relation- 
ship between  link  volumes  and  link  factors,   A  regression 
analysis  of  the  minimum  AADT  associated  with  a  link  on  the 
Phase  II  link  factor  was  performed  to  develop  this  relation- 
ship. 

The  minimum  AADT  on  a  link  (the  lowest  AADT  of  a  sec- 
tion of  the  link)  was  used  as  the  dependent  variable  because 
the  minimum  value  would  more  closely  reflect  intercity 
travel  as  shown  in  Figure  22.   In  this  idealized  relationship, 
it  can  be  seen  that  local  traffic  increases  as  a  centroid 
or  route  intersection  is  approached  but  intercity  traffic 
remains  constant.   Thus,  the  minimum  traffic  volume  more 
closely  represents  what  the  ITDF  model  measures. 

Only  the  Phase  II  link  factor  was  used  in  the  link 
analysis  for  tv;o  reasons;  namely  over  three-fourths  of  the 
links  had  zero  values  for  Phase  I  and  Phase  III  factors  and 
the  results  of  the  centroid  analysis  showed  that  the  numer- 
ical values  of  the  factors  from  the  three  phases  were  not 
additive.   The  establishment  of  a  relationship  between  the 
three  factors  is  a  research  topic  in  itself. 

The  link  data  used  in  this  analysis  were  selected  to 
meet  one  of  two  criteria:   the  link  was  located  in  an  area 
remote  from  Interstate  highway  rights-of-way  or  the  link  was 
located  near  portions  of  the  Interstate  system  completed 
prior  to  1962.   These  criteria  were  imposed  in  an  attempt  to 
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insure  that  travel  patterns  in  the  areas  of  the  selected 
links  would  not  be  significantly  affected  by  completion  of 
the  entire  Interstate  System. 

After  the  link  data  had  been  collected  and  plotted  on 
a  scatter  diagram,  it  was  decided  to  eliminate  those  data 
sets  having  a  Phase  II  link  factor  of  less  than  50  and  to 
try  a  functional  relationship  of  the  form: 

Y  =  a  +  b  Log(X) 

where  Y  is  the  minimum  link  AADT  and  X  is  the  Phase  II  link 
factor. 

A  total  of  126  data  sets  were  used  in  this  analysis. 
The  regression  model  developed  was: 


Y  =  -  8977  +  5523  Log, q (X) 

2 

which  had  an  R   of  0.919.   A  plot  of  this  equation  is  shown 

in  Figure  23. 

2 
The  variance  of  the  estimate,  s^  ,  of  a  predicted  Y 

Y 

for  a  given  X  is  estimated  by: 

S.2  =  s  2   u  +  ^^  ''"^  -  ^^'  ^ 


Y      ^  "    Z  (X  -  X)^ 

For  X  equal  to  the  mean,  the  variance  of  the  estimate  is 
given  by: 

S.2  =  551  412  [1  +  J:-  +  (2.377  -  2.377)^ 
s^     ^^31,^1^  li  +  ^26  +        25.43      ^ 
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and  the  standard  error  of  the  estimate,  s  ,  equals  745.5. 

Y 
The  standard  error  of  the  slope  was  147.3  so  that 

the  95  per  cent  confidence  interval  extended  from  5,234  to 
5,812. 

Based  on  this  analysis,  it  was  decided  that  the  regres- 
sion model  is  an  adequate  means  of  predicting  the  traffic 
volume  associated  with  the  Phase  II  link  factor.   Further- 
more, these  results  reinforce  the  conclusion  that  the  ITDF 
model  used  is  adequate  for  the  synthesis  of  travel  demand. 
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CONCLUSIONS  AND  APPLICATION  OF  RESULTS 

Conclusions 
Conclusions  which  were  reached  during  the  course  of 
this  research  included  the  following: 

1.  The  tree  type  decomposition  algorithm  is  a  prac- 
tical method  for  computing  minim.um  path  trees  for  large  net- 
works and  can  result  in  a  substantial  reduction  in  computer 
time  if  the  number  of  tree  tables  required  is  not  too  small. 

2.  The  Intercity  Travel  Desire  Model  used  in  this 
research  is  an  adequate  miethod  of  synthesizing  travel  demand. 

3.  The  Phase  I  and  Phase  III  portions  of  the  determin- 
ation of  intercity  travel  desire  reflected  thru  traffic 
demand  but  not  in  a  manner  which  permitted  summation  with 
each  other  or  with  the  factors  of  Phase  II.    Furthermore, 
these  Phases  need  not  be  determined  in  any  future  study  as 
the  factors  found  from  Phase  II  were  adequate  for  classifica- 
tion purposes. 

4.  For  cities  outside  of  the  delimited  study  area, 
only  those  with  a  population  of  more  than  10,000  or  20,000 
should  be  considered  in  calculating  the  travel  desire  factors 
Factors  for  the  smaller  cities  are  virtually  negligible  be- 
cause of  the  long  distances  involved. 


136 


5.   The  relationship: 

Y  =  -  8997  +  5523  Log, „  (X) 

where  Y  is  the  AADT  and  X  is  the  Phase  II  link 
factor,  was  found  to  be  an  acceptable  highway 
classification  tool  for  Phase  II  link  factors  of 
65  or  greater.   The  limiting  value  of  65  represents 
a  traffic  volume  of  about  1,000  vehicles  per  day. 
Some  representative  volumes  and  the  factors  found  to 
be  associated  with  them  are: 

Volume  Factor 

1,000  65 

2,000  97 

3,000  148 

4,000  225 

6,000  522 

10,000  2,700 

Application  of  Results 

It  is  proposed  that  the  State  Highway  System  of  Indi- 
ana be  subclassif ied  into  four  designated  systems.   These 
four  subsystems  are: 

1.   Principal  State  Highway  System.   This  system  should 
be  composed  of  the  presently  designated  Interstate  highway 
system  and  those  highv/ays  which,  on  the  basis  of  their 
Phase  II  Intercity  Travel  Desire  Factor  and  other  planning 
criteria  should  be  reconstructed  to  freeway  standards  by  1982, 
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The  design  year  was  selected  as  1982  because  the  AADT-link 
factor  relationship  was  based  on  1962  volume  data  thus  pro- 
viding a  20  year  planning  interval. 

2.  Primary  State  Highway  System.   This  system  should 
be  composed  of  the  additional  highways  required  to  provide 
for  the  interconnection  with  the  Principal  State  Highway 
System  of  all  Indiana  cities  over  5000  population  and  should 
also  include  those  roads  having  a  Phase  II  link  factor  of 

at  least  125. 

3.  Secondary  State  Highway  System.   This  system  should 
be  composed  of  the  additional  highways  required  to  provide 
for  the  interconnection  with  the  previously  designated  sub- 
systems of  all  still  unconnected  county  seats  and  should 
also  include  those  highways  having  a  Phase  II  link  factor  of 
at  least  50. 

4.  Collector  State  Highway  System.   This  system  should 
include  the  remainder  of  the  present  State  Highway  System 

of  Indiana  not  already  in  a  previously  designated  subsystem. 

The  proposed  Principal  State  Highway  System  for  1972 
and  1982  is  shown  in  Figure  24.   The  routes  proposed  as  sup- 
plementary to  the  Interstate  system  are  shown  in  two  groups. 
The  first  consists  of  freeways  deemed  necessary  by  1972  and 
the  second,  freeways  deemed  necessary  by  1982. 

In  general,  when  the  AADT  of  a  two-lane,  two-way  high- 
way lies  somewhere  about  7,000  to  8,000  vehicles  per  day, 
serious  consideration  should  be  given  to  making  it  a 
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multilane  facility  with  the  degree  of  access  control  deter- 
mined by  individual  study.   If  the  design  hour  volume  (DHV) 
is  taken  at  about  15  per  cent  of  the  AADT,  its  magnitude  is 
1,050  to  1,200  vehicles  per  hour  for  AADT ' s  of  7,000  to 
8,000  respectively. 

According  to  the  "Highway  Capacity  Manual"  (56),  the 
service  volume  for  level  of  service  C  approaches  1400  pas- 
senger cars  per  hour  under  ideal  conditions  while  the 
service  volume  for  level  of  service  B  approaches  900  pas- 
senger cars  per  hour  under  ideal  conditions.   Level  of 
service  B  is  associated  with  the  design  of  rural  highways 
while  level  of  service  C  merely  implies  satisfactory  oper- 
ating conditions. 

The  service  volume  at  level  of  service  C  on  a  two-lane 
highway  with  adequate  lane  and  shoulder  v;idth,  no  passing 
sight  distance  or  alinement  restrictions,  and  only  10  per 
cent  trucks  is  about  1,200  vehicles  per  hour.   Thus,  demand 
volumes  of  about  1,100  to  1,200  vehicles  per  hour  virtually 
dictate  multilane  design  even  if  ideal  conditions  can  be 
designed  into  a  two-lane  highway  and  if  a  desirable  level 
of  service  is  to  be  attained. 

The  selection  of  the  Phase  II  link  factor  limits  used 
to  select  this  subsystem  of  highways  was  based  on  a  pro- 
jected volume  of  7500  vehicles  per  day  assuming  an  annual 
average  growth  rate  of  four  per  cent.   Highways  in  Indiana 
which  will  be  carrying  this  volume  of  traffic  by  1972  are 
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those  which  had  a  link  factor  of  350  or  more  as  calculated 
in  this  study.   Those  which  had  a  link  factor  of  180  to  350 
will  require  multilane  design  between  1972  and  1982. 

The  four  subsystems  of  the  State  Highway  System  of 
Indiana  as  suggested  by  this  research  for  1972  are  shov;n  in 
Figure  25. 

The  Primary  State  Highway  System  includes  those  high- 
ways with  link  factors  between  125  and  350.   The  lower  value 
is  representative  of  a  1962  traffic  volume  of  about  2,600 
vehicles  per  day  so  that,  with  a  four  per  cent  annual  growth 
rate,  the  1972  volume  would  be  about  4,000  vehicles  per  day. 

This  projected  volume  of  4,000  was  used  because  it 
represents  the  point  at  which  high  standards  for  two-lane 
rural  highway  design  are  often  recommended  by  state  highway 
departments.   This  proposed  system  for  1972  also  includes 
those  highways  with  link  factors  between  180  and  350  which 
are  suggested  for  transfer  to  the  Principal  System  between 
1972  and  1982. 

The  Primary  System  was  also  selected  to  ensure  the 
interconnectivity  of  cities  over  5,000  with  the  Principal 
System.   Similar  nearby  cities  in  adjacent  states  as  well 
as  the  major  routes  of  adjacent  states  were  also  considered. 

The  Secondary  State  Highway  System  was  selected  to  pro- 
vide for  the  interconnection  of  all  county  seats  and  also 
to  provide  a  general  coverage  of  all  areas  of  Indiana.   A 
minimum  link  factor  of  50  was  used  because  its  use  appeared 
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to  provide  good  overall  service  to  most  areas  of  the  State 
and  included  most  highways  with  a  1972  volume  of  approxi- 
mately 1,000  vehicles  per  day  or  more.   Most  of  the  small 
Indiana  cities  would  be  served  by  this  system. 

With  regard  to  the  selection  of  each  of  the  four 
subsystems,  it  is  important  to  note  that  the  factor  limits 
of  each  subsystem  were  not  strictly  adhered  to  in  all  cases. 
Each  subsystem  was  selected  so  that  it  was  integrated  with 
previously  designated- subsystems  and  so  that  there  were  no 
isolated  sections  unconnected  to  either  an  equal  or  higher 
system. 

Thus,  some  flexibility  in  selecting  the  subsystems  was 
necessary  to  achieve  an  integrated,  interconnected  highway 
system  serving  all  areas  of  the  State. 

The  remaining  highways  in  the  current  State  highv;ay 
system,  after  selection  of  the  three  subsystems  discussed 
above,  includes  some  highway  sections  which  should  be  con- 
sidered for  deletion  from  the  State  highway  system.   Their 
low,  in  some  cases  zero,  link  factor  indicates  that  the 
service  they  provide  is  local  in  character  and  that  the 
counties,  rather  than  the  State,  should  be  responsible  for 
them.   All  the  remaining  highways,  however,  are  shown  in 
Figure  25  as  belonging  in  the  Collector  System. 
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RECOMMENDATIONS  FOR  FURTHER  RESEARCH 

Based  on  the  preceding  analyses  and  conclusions,  the 
following  further  research  is  recommended: 

1.  The  optimum  number  of  partitionings  of  large  net- 
works and  the  influence  of  the  relative  size  of  each  parti- 
tioning should  be  determined  for  use  of  the  tree  type  decom- 
position algorithm  with  regard  to  both  the  accessibility  of 
the  required  data  and  the  savings  of  computer  time. 

2.  The  relative  efficiency  of  both  computer  storage 
and  time  requirements  should  be  determined  for  the  tree  type 
decomposition  algorithm  with  respect  to  matrix  type  decomposi- 
tion algorithms  and  also  with  respect  to  the  basic  minimum 
path  algorithm  when  programmed  for  high  storage  capacity 
computers . 

3.  The  results  of  this  study  should  be  updated  and 
the  synthesis  procedures  refined  by  means  of  a  comprehensive 
traffic  survey  covering  the  entire  State  of  Indiana  or  at 
least  a  portion  of  it.   This  traffic  survey  should  be  con- 
ducted in  1970  when  data  on  population  will  be  gathered  by 
the  U.  S.  Bureau  of  Census.   The  traffic  survey  should  include 
sample  data  on  volumes,  out-of-state  vehicles,  truck  traffic, 
roadside  origin-destination  surveys  along  rural  roadways  and 
at  cordon  lines  around  cities,  and  other  desired  information. 
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These  data  will  then  enable  conduct  of  the  following  needed 
research  on  intercity  traffic  desire  models: 

a.  Determination  of  the  effect  of  distance  between 
cities  on  intercity  traffic  volumes. 

b.  Determination  of  the  best  intercity  traffic  desire 
model  for  estimation  of  intercity  traffic  volumes. 

4.  Land  use  or  other  type  models  designed  to  explain 
the  amount  of  locally  generated  traffic  entering  urban  areas 
and  that  which  exists  on  rural  highways  should  be  developed. 
These  models  should  then  be  combined  with  an  intercity  travel 
desire  model  to  give  a  more  comprehensive  picture  of  the 
relative  importance  of  highways. 

5.  The  link  factors  developed  in  the  Phase  II  portion 
of  this  study  should  be  combined  with  the  Sufficiency  Ratings 
of  the  various  highway  sections  to  calculate  the  priority  of 
highway  sections  for  improvement  or  reconstruction  to  adequate 
standards  or  for  the  preparation  of  plans  for  construction  of 
new  facilities  to  replace  existing  highways. 

5.   An  economic  model  should  be  devised  to  properly 
allocate  available  funds  to  the  four  State  highway  system 
subclasses  proposed  in  this  study. 

7.   A  study  is  recommended  to  determine  the  necessary 
additions  to  the  proposed  Primary  and  Secondary  State  High- 
way Subsystems  in  Indiana  to  accommodate  such  special  pur- 
poses as  service  to  airports,  institutions,  recreational 
areas,  and  military  installations. 
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'14' 

17. A 

70  .  ' 

0.0 

781  ■; 

'06q 

10.4 

'^  7.  1 

0.0 

3240 

'25  1 

1  '  .  9 

1  .3 

".0 

7'^i'4 

'OAq 

=-  .  8 

8  8 ,  q 

0,0 

^242 

'24' 

'.4 

7'=.'  .0 

0.0 

2'07 

'1  'R 

10.  A 

7  p .  q 

0.0 

3243 

3244 

10.9 

7.0 

'^,^ 

2486 

'07  7 

6.7 

115.7 

^  A 

2315 

'24  5 

0.  7 

47.9 

0.0 

3'16 

"17 

7.7 

0.0 

0,0 

2  421 

''Op 

7  7.9 

7.4 

c.o 

2606 

7A1  8 

6.q 

0.0 

0,0 

7  044 

7^71 

.  q 

18.4 

c.o 

'142 

'1^2 

1  -.7 

09.8 

0,0 

2  7  44 

3'i  g 

17.' 

40  .q 

0.0 

'1=^2 

llr-,K 

18.7 

'0.0 

0,0 

^■^70 

"  7i, 

'.4 

^7.7 

<^.l> 

747A 

'1  At; 

6.6 

4  6,4 

0,0 

74^4 

'  -  2A 

77.  A 

^  p 

C.r- 

7S2R 

7A01 

'  .  ' 

40.6 

0,0 

2452 

^  196 

A.  ' 

71  .4 

".0 

74Aq 

'174 

8.0 

41  .7 

0,0 

25P7 

'03  7 

1^.  1 

2'7.7 

0.0 

31  77 

'1  2i 

7.7 

770,8 

0.0 

2'^a7 

761  ' 

4.6 

'37.7 

c.o 

3117 

'17' 

6.P 

1  7,P 

0,0 

2612 

761? 

l.B 

0.0 

0.0 

'noq 

'11-' 

1  7.Q 

0.0 

0,0 

2''4'- 

'75'' 

17.7 

14.7 

:  ,0 

'ooq 

'  1  no 

7.  1 

0.0 

0,0 

'231 

ii'" 

o.  ■; 

■^7.4 

c.o 

7R'7 

'  1  00 

1  7.8 

o.o  ■ 

0,0 

3202 

3  20  3 

25.9 

0.0 

c.o 

27R7 

'0'3 

23.  1 

CO 

CO 

'201 

'  70  7 

14.0 

o.o 

r  n 

'1Q7 

4101 

1  '  .0 

7  7  .  A 

0,0 

7  '^9 

'77'^ 

2  7.0 

71  .P 

0.0 

2  '  0  7 

7  1  P  7 

1  3.4 

7  .9  ,  A 

0.0 

2'P6 

7'DQ 

16.  P 

1:5." 

0.0 

71  'q 

7707 

6.4 

7  A  ,  A 

n  ,  0 

'264 

T  ,->  A  P 

Q  .  A 

-A. 6 

0.  0 

7  1  'O 

7  7  on 

18. A 

1  A  ■  .  i. 

'1  "  .  4 
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^274 

^276 

4.6 

157.5 

0.0 

2300 

1120 

1'=..4 

113,0 

0,0 

1?76 

^27n 

lfl.8 

0.0 

0.0 

2017 

7111 

18,0 

=  3,5 

0.0 

?031 

T  179 

20.1 

fll  .0 

0.0 

2012 

1193 

,2 

289.9 

0.0 

?n^l 

7-J3B 

.  ^ 

10  1,7 

0,  0 

ai  01 

1 1  04 

I''.  ■^ 

17,"; 

0.0 

7^?fl 

•'214 

l-^."^ 

101  .-■ 

0.0 

74'^o 

T]04 

•^ .  4 

1  7  ,  t 

0.0 

?171 

■'417 

1R.4 

71^.7 

74  .  "= 

744^ 

1111 

1  7  .  T 

4.1 

0.0 

?4lfl 

3293 

17.7 

11.7 

0.0 

2140 

1711 

1  1  .P 

10.^ 

0.0 

3?q3 

9  7QQ 

15.0 

0.1 

0.0 

3210 

1211 

1.6 

48.1 

0.0 

?^7'^ 

3^17 

5.  "^ 

.p 

C  ,0 

7106 

17  10 

11.° 

74.'^ 

0.0 

3^1? 

13)1 

6.7 

5.4 

o.p 

71O6 

1717 

70.0 

8.4 

0.0 

?Aiq 

131^ 

■^.0 

110.7 

A  7  .  f. 

2711 

17^2 

16.1 

17.2 

0.0 

20^5 

2419 

1.  0 

190.5 

67.6 

3139 

1  140 

1.9 

606.6 

0.0 

?nu.'=> 

331  R 

17.  1 

1  87.0 

0.0 

20QR 

1  1  10 

^  .^ 

165,1 

0.0 

?^7q 

^31  R 

1  7.6 

155.1 

0.0 

200R 

7700 

10. R 

717.1 

0.0 

■J?TQ 

41-SQ 

16.0 

.tc.7 

0.0 

7  70^ 

7TOC 

1  7  .  Q 

7^.7 

0.0 

■^7^^ 

7.?-^L 

A.  B 

0^.5 

0,0 

711R 

1104 

10.7 

47, R 

0.0 

?^4P 

^234 

7.6 

77.0 

0.0 

24  19 

1104 

5.6 

P8,o 

0.0 

IN|0    OTHFR  HIGHWAV^ 

3?38 

3239 

6.9 

0.0 

0.0 

2579 

75R7 

24.0 

44,5 

0.0 

2060 

^239 

1.0 

0.0 

0.0 

2570 

2619 

1.8 

0,0 

0.0 

5060 

2205 

4.0 

A  ^  0 

r.  0 

7  100 

7747 

11.7 

o.n 

0,0 

7007 

^189 

15. R 

60.7 

0.0 

7410 

117R 

11.4 

10,6 

0.0 

?62'=. 

TIRO 

11.0 

17.6 

0,0 

2426 

13T^ 

16.1 

0,0 

0.0 

7377 

T  10  = 

l=.i 

77  .  D 

^"",5 

2426 

1174 

17,4 

16,1 

0.0 

3306 

•Jin 

1  1  .4 

141.5 

67,*^ 

mo 

1111 

1  P.  1 

14.5 

0,0 

2377 

33n 

16. P 

0.0 

0.0 

3165 

1168 

7.'^ 

15,5 

0,0 

2347 

3279 

R.O 

7.1 

CO 

2224 

116P 

n  .4 

10,6 

0,0 

7347 

3276 

17.3 

1  14.9 

0.0 

2774 

1776 

7.1 

11.'=^ 

0.0 

377S 

1276 

17.1 

54.6 

0,0 

3776 

1777 

H  ,  7 

7.1 

0,0 

7?P0 

177=, 

10.  7 

15.0 

'^.  0 

7104 

1710 

<..  1 

0,0 

0  .0 

215R 

2617 

6.1 

n.o 

0,0 

2176 

3112 

4.1 

4,6 

0,0 

^74R 

3240 

11.0 

71.6 

0,0 

24'.R 

11  97 

1  A,7 

44,7 

0,0 

2^16 

3247 

11.6 

77.7 

0.0 

2457 

^107 

70,  0 

10,^ 

0,0 

7021 

7317 

7.  1 

0,0 

0.0 

2359 

1160 

P,8 

15,5 

0.0 

2029 

3335 

R.  P 

0.0 

0.0 

2159 

116? 

1  0,4 

61,7 

0,0 

7029 

737c; 

1  .'"' 

^  1  R  .4 

^.0 

271'=. 

1167 

4,4 

°P  ,  7 

0,0 

7^7"^ 

121R 

1.0 

47A.0 

0,0 

1100 

-JTO  7 

'  ■=•  .  0 

17  ,  R 

0,0 

3711 

121B 

14. R 

•^47. 7 

C.O 

7177 

177  1 

6,  1 

1  ,7 

0,0 

329fl 

1 100 

14.  T 

7^  .  =; 

0,  0 

371  A 

1791 

4,"' 

0,0 

0,0 

?43n 

i3on 

74.7 

101  .0 

0.0 

1717 

1716 

0,0 

0,0 

0  ,0 

7446 

1257 

?7.7 

0.0 

0.0 

3261 

1764 

1  6, A 

59,4 

0,0 

2416 

32R9 

10.7 

54.4 

0.0 

2383 

3263 

10,7 

t4,  R 

0,0 

3288 

3289 

6.2 

41  .7 

0.0 

3211 

1215 

10,4 

1.7 

0,0 

2451 

3200 

24.  1 

6.0 

0.0 

2120 

1214 

10,6 

7  1  7  ,  Q 

0,0 

2  066 

7076 

11.6 

1  R'^.9 

p.? 

2';07 

1070 

14,0 

R7.6 

0,0 

2066 

7710 

7.0 

10  7.6 

-'.  0 

3070 

40RI 

1  7,0 

pi^.P 

0,0 

7210 

7710 

0.  B 

00.0 

A,  0 

2  76A 

1  1  1R 

1,7 

44,0 

0  ,0 

7701 

5210 

7.7 

Rl  .^ 

0.0 

2A70 

1 1  IP 

17,6 

0,0 

0,0 

2201 

3247 

11.2 

45.0 

0.0 

2417 

2  44? 

36.1 

R1.5 

0,0 

2531 

1105 

P.  1 

10.1 

C.O 

247'=. 

3000 

9.1 

1Q5.2 

0,0 

3105 

3107 

0.4 

33.0 

C.O 

2078 

1184 

17.0 

59.6 

0,0 

2533 

3064 

R.6 

0.0 

0.0 

2420 

1107 

70.4 

1  7  .A 

0,0 

7  "^^T 

1)01 

R  .  ~' 

^0.1 

0.0 

2141 

17AR 

°  .  7 

c  A.5 

0,0 

^10' 

1 1  00 

R  .  7 

51.0 

0.0 

74^0 

1  7^=, 

1  7,0 

"'4  .  0 

0.0 

3109 

3114 

10.4 

74.6 

0.0 

23R1 

2415 

14,1 

'^4  .  4 

0,0 

31  la 

1117 

R.O 

07.4 

0,  0 

loni 

1001 

7,6 

0,0 

0,0 
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^117 

•il?0 

^.4 

7r..i 

0.0 

?'^2U 

7S7C) 

4.2 

77.0 

".0 

^0fl4 

&0?S 

!7.n 

^^.7 

r.o 

2592 

^084 

9.1 

4^.4 

r.o 

?S9? 

^.■;7ft 

11.1 

4.4 

0.0 

7^71. 

^07*, 

8.0 

4.4 

o.n 

2^05 

2^.7T 

9.0 

77.3 

0.0 

Tlli,? 

■»0i.^ 

6.7 

fin,c, 

0.  o 

2560 

^OitT 

4.6 

16'='.1 

r.o 

2560 

2  587 

in.  P. 

5^.6 

o.n 

I 

^^n     OTHFR  HIT. 

HWA  YC 

?0C1 

217P 

17.9 

0.0 

CO 

212S 

2627 

14.0 

n.o 

0.0 

^16'^ 

2  IfiQ 

14.6 

fil  .Q 

0.0 

716A, 

2627 

14.2 

0.0 

0  .  O 

?167 

26?o 

.  0 

n.o 

0.  " 

216P 

2627 

16. P 

o.P 

.4 

2  160 

7627 

4.6 

0.0 

0.0 

2169 

262n 

6.1 

n.o 

0.0 

2170 

2179 

14.1 

0.0 

CO 

21^2 

2626 

7.4 

^5^.3 

0.0 

2172 

2627 

6.1/ 

0.0 

CO 

2A71  ^7«fl  21.4 

'^?!iP     3290  10.3 

2510  3072  1';.4 

2576  7=^90  3.0 

71^7  7307  6.4 

209ft  75P9  n.6 

2^07  41S«  47.0 

2312  4216  78. O 

3187  4169  16.0 

3314  4773  1,n 


21?^  76^0 

7174  7676 

7174  7ft30 

7176  7678 

21 77  7ft78  1 O.q 

21  70  7A56    1.1^ 

2187  7188  l=i.8 

2187  2628   7.8 

21°8  2678  I'^-.O 

2676  7678   -[.-i 

2h21  2f^?^     1^.8 


^.6 

1  .9 

4.6 

1  1  .0 


7?. 7 

74.6 

^1.2 

790,9 

445,7 

0.0 

'^6.7 

TO  ,  1 

14.0 

37,7 

10.' 

87]  ,4 

17,7 

48,2 

1  '•'^.1 

700  ,  ? 

0.0 

0.0 

0.0 

713,3 

0.0 


0,0 
0.0 
O.O 

0.0 
75.'^ 
0.0 
0.0 
0.0 
0.0 
0.0 

0.0 
0,0 

0,0 

0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
0.0 
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Mawr         MODF  NO'-iP  AMTF     ITHF 


rPMTon I n 
MAMF 


HWV 
Moric 


r^MT  nl  c,T  OHACF   J  I 

NODF  AMTr     ITTF 


NEW  BUFFALO. Ml   ?felH 


IR.?       CROWN    POINT 


AN'^ONI&.OH 
AN'^ONI  A  .OH 
IIMIOM    riTV,'T-l 
HARRt  SnN»nH 
HARRISON«nH 
CLFVFS.OH 
CLFVF<^,nH 
NlFW     PAD  I  c;  , OH 
MCi,i    DAP  I  c;  ,riH 
LOHANSrOOT 
''nrHF^TFP 
RFMlNriTON 

rfn?;<;fi_afp 

CtAPY 
GARY 
LOWFLL 
inwFl.  L 
rrnAD    i^v:tr 
fFHAP    LAKF 


^1  ?q 
i.1  '^'^ 
1 1  70 
;'0  15 
7016 
1  749 
7454 
7  14  1 
7^07 
74f<1 
777R 
74n-i 
74ns 
71  "^4 
71  S6 

7  1'^A 
7'^iS7 


1  4 
1  64 
1  64 
1  66 
1B5 
IR^ 
1R6 
1  fl6 
744 
744 
6":,  6 
6'=  ' 
66R 
66Cj 
66^' 
66  r 
661 
661 

667 


4  .n 
1  "^  .n 

<;  ,n 
1  .f^ 
"^  .0 

■^  .n 
7  .n 
fi  ,n 

n  ,0 


1  (-1  ,n 

1  n.n 


1  ri.R 
n  ,6 


47.  1 

n.n 
^  hi .1 

RR.q 
n.n 
o.n 

R0.9 

R4.  1 
•=  ^^4.7 
7  4  s  .  1 
17  7.'=. 
^  7R  ,  O 
71^^.4 
"^■=•6.7 

74.^. 

7=  .4 
117.1 
7 '^0.7 


^.T,  JOHN 
■^T,  JOHN 
LA  PORTF 
"fHIGAM  riTV 
K  I  Nfi'^FOPn  HT"^. 
KING^FOPn  HT';. 
DL YMOUT  H 
rULVFP 
APfnO'^ 

MORorro 
rvOnoLAMn 

krFNTLANn 

k-  F  N'  T  L  A  N  0 

HF^ROM 

VAlPAPA  I  CO 

KOilT<: 

rHFSTFRT'^M 

CHF';tfRTOM 

FOANirt  <;\/ilL'^ 


2412       66^ 
7T67       664 


7^6° 
7497 
760R 
7490 
7491 
256? 

7?Rn 
?-^67 
7S45 
7240 

774] 
74]  1 
?ri?6 
^01  1 
2149 
2160 
"■404 
77  56 


664 
66'=' 
666 
66  7 
667 
668 
669 
670 
671 
677 
67? 
67? 
674 
676 
676 
677 
677 
A70 


1  .0 

7.5 

4.? 

.5 


1  .n 
1  .0 
?.? 
4.6 


6  .  1 
1  .1 


474.0 

69,  i, 

77.  1 

770.7 

1 006.4 

5R.  7 

R7.4 

?71  .5 

115.9 

14R.n 

114. ■=> 

114. R 

116.9 

79. R 

1-70.0 

57?. 8 

96.9 

72.4 

45T.  7 

147.'^ 

41  ."^ 


FRa\|rF<;\/I  lLF 

SOUTH     PFNIO 

NORTH    LIRFPTY 

WALKFRTON 

NFw    CARL  I  c-LF 

^l  F  1.1    r  A  D  L  I  c  i_  r 

L  vr  irir 

LYDirr. 

KNOX 

NORTH    Jiin<;o\i 

.•^ONON 

MONT ITFt  LO 

MOM     ;  ("FLLO 

nF»/OTTF 

ppPMCTTCi/T  lie 

PlIPNFTT'^VILLF 

FORT  WAYNF 

^]FW  HAVFN 

MONROEVILLF 

MOMROFVI LLF 


7  7  67 

77R7 

?n?i 

76^? 
7  6  '^  4 
76*^4 
761  6 
74P  7 

?ni6 

7755 
775? 
7  7  ■=  4 
7409 
■'4R7 
7  5  44 

7517 

77?0 
7557 


6  "79 

6an 
6R1 
602 
6R7 
60^ 

AQ4 

6R4 
60  = 
6P6 
6R7 
6OR 
6n  Q 
6fl9 

690 
691 

697 

69? 
69? 


4.1 


.7 

1  •  R 


1  =  .6 

1  .^  .  R 

1  .0 

,  7 

14.-' 
6.0 


49,6 

1 67?. 6 

1  ?  ^  ,  R 

1  RR  .7 

60.4 
1  '■  6  .  ? 

40,= 

R  6  ,  ? 
1  R  7  ,  4 
1'  0  .  4 
1  1  1  .R 

R  R  ,  R 
171.7 

9  n  ,  4 

4A  ,  P 
4^.4 
17  19.5 
7M  .1 
3R.6 
R9,q 


CYPACll^F 
<;VR  ATI  ic  F 
WILF IPO 
MI LFOPp 

M  I  LF^R^^ 

P  I  pprcjOM 

PI  cor F TOM 

WARSAW 

■.M  h\n•\l^     I  A"  F 

W!  NON'*     L'*'  F 

L  A  O  R  A  N  o  F 

POUPPON 

P  C  F  M  F  M 

DF  Pi  I 

At   P  1  r.M 

KFNioAl  LVILLF 

I_IGnNIFP 

O'^CFpL  A 

05CF0LA 

ANGOLA 


7  =  79 

?07? 

9^? 

?  =  R0 
7  =  57 
T067 
755R 
7667 
255R 
2594 
2560 
26P2 
2  =  4^ 

TOP  = 

7674 
7571 
7  =  7H 
7579 

7'='10 


704 
704 
705 

70  = 
706 
70<^ 
706 
707 
70R 
70R 
709 

710 
711 
717 

71  -i 
71  4 
7  1  ■-. 
"'1  6 
716 
71  7 


R  ,0 

R.7 

19.1 

4.7 
7.7 

1  =  ,  n 
,  = 


.  7 
,  7 


1  .4 

16.0 


110.5 

4=  .  4 

1  ■;  .0 

47.  1 

7  9.0 

1  77,  f, 

1  f,.^ 

47  )  .8 

97.7 

107.2 

157.4 

167.4 

197,9 

6  4  4.7 

170,= 

?  7  ?  .  = 

217.6 

60.7 

44.  W 

7  7  R  ,  7 
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WATFRLnn 

fli  miiPN 

r.  A  p  (?  c-  T  T 

GARRf^TT 

BUTLFR 

FLk.HART 

DUNLAP 

DUNLAP 

GO'^HFM 

NAPPANFF 

WA<ARmi;a 

WAK  API  ic  A 

HUNTINGTON 

WOLCOTTVILL'^ 

LAKF  JAMFS 

LAKF  JAMFS 

OFCATUR 

MONTPFLIFR 

HAPTFORD  riT^ 

f-'lJMCIF 

ALBANY 

YORKTOWN 

YORKTOWN 

YORKTnwN 

r^AL'^VILLF 

HALFVILLF 

FATON 

FATON 

MARION 

FAIRMOlINT 

CAIRVOUNT 

MPLA^!Cl 

tIPLANn 

r,a<;  CITY 

CARMFL 

ARCADIA 

ARCADIA 

ARCAOI A 

ARCADIA 

CICFRQ 

cir=^Rn 

NORL'^SVILLF 

VfFcTFIFLD 

SHIRLFY 

5HIRLEY 

SHIRLFY 


?sr.7 

?57? 

25?6 

?')?"' 
■JOTS 

75-^7 

?1?6 

^084 
2501 
^114 
^121 
246'- 
312'i, 
735 
2104 
2465 

2  104 
246"^ 
24(S7 
3121 
3092 
2107 
^1  ^fi 
-^  1  77 
31  2^ 
2471 
^1  3R 

743 
o*,0 
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^T^T 

1  10? 
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15.5 
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■^?Qf. 
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46,4 
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1103 

71,6 
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HA7FL  TOWN 

7137 

1  08^ 
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OTWPLL 

T^OO 
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IQ.i 
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^3  03 
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1^  1  .  7 

HATFIELD 

1  lO-^ 

1  104 

7,9 

0,4 

PATOKA 

7337 
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6.1 

71  .T 

HATF IFLD 

2^7] 

1  104 

8,9 

19.5 

PALMYRA 

7431 

1085 

.7 

37.5 

HATE  I FLD 
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1  104 

4,7 
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VFDOPA 

7  44  ■> 

ion  6 

6.0 

10.1 

RICHLAND 

7T'l 

1  105 

7.7 

18,4 
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^333 
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6.5 

35  .Q 

P I OHL AND 

7417 

1  105 

10.9 

19.  7 

N^OMROF  riTY 

3  30? 

1  087 

16.1 

1  1  .« 

5T,  MPIMRAO 

7470 

1  106 

0,7 

51.0 
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APPENDIX  B 

COMPUTER  PROGRAMS  FOR  DETERMINATION  OF 
INTERCITY  TRAVEL  DESIRE  FACTORS 


Bl.   Small  Trees 

The  mechanics  of  Small  Trees  will  be  outlined  briefly. 

Four  storage  areas.  A,  B,  C,  and  K  were  set  up  where 
storage  areas  A  and  K  contained  the  link  and  centroid  tables 
respectively.   A  total  of  3,618  link  description  cards  v/ere 
required  to  simulate  the  1,809  links  making  up  Indiana's 
highway  system.   The  i    link  card  consisted  of  lANODA. 
(the  A-node  in  storage  A),  IBNODA.,  and  LABELA.  (the  distance 
between  the  two  nodes) .   The  sequence  of  link  cards  was  sub- 
ject to  only  one  restriction.   All  cards  having  an  equal 
lANODA  had  to  be  read  into  storage  consecutively.   However, 
the  actual  numerical  order  of  lANODA  card  groups  was  not 
important. 

The  centroid  table,  consisting  of  687  link  cards,  in- 
cluded only  the  back  link  descriptions,  i.e.,  the  link  from 
a  highway  node  to  a  centroid.   The  sequence  of  these  cards 
was  not  important.   The  forward  link  cards  were  read  in 
with  each  source  centroid  and  were  needed  to  generate  the 
initial  connection  of  the  source  centroid  V\7ith  the  highway 
network.   Each  centroid  card  consisted  of  lANODK,  IBNODK, 
and  DLINK. 
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Storage  area  B  contained  the  tree  table  with  three 
columns  designated  as  lANODB ,  IBNODB,  and  LABELB .   The 
lANODB  column  was  a  listing  of  each  centroid  and  node  number 
in  the  network.   Centroid  numbers  were  generated  as  a  con- 
secutive numerical  sequence  from  656  to  1,106  while  the 
node  numbers  were  generated  from  the  lANODA  column  of  the 
link  table.   A  total  of  1,477  centroids  and  nodes  were  thus 
generated. 

IBNODB  is  the  next  node  on  the  minimum  path  node  se- 
quence from  lANODB  to  the  source  centroid  while  LABELB  is 
the  distance  from  lANODB  to  the  source  centroid.   In  the 
section  on  the  Basic  Minimum  Path  Algorithm  of   a   previous 
chapter,  it  was  stated  that  the  label  on  node  i  is  changed 
and  the  node  added  to  a  set  S  provided  that: 

V .  >  min  V .  +  d . . 
^    jcS   ^     ^^ 

or  in  the  context  of  this  section: 

LABELB .  >  min  LABELB .  +  LABELA . . 

If  the  label  on  node  i  is  changed  then  the  following  change 
must  also  be  made: 

IBNODB .   =   lANODB . 
1  D 

Storage  area  C  contains  the  nodes  in  the  set  S  designated 
as  lANODC  and  the  label  on  each  node  in  S,  LABELC. 
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A  significant  amount  of  computer  time  in  calculating 
a  tree  table  is  consumed  in  searching  the  several  node  arrays 
to  find  the  location  of  certain  node  numbers.   In  order  to 
reduce  the  magnitude  of  search  time,  the  following  substi- 
tutions were  made  after  the  lANODB  array  of  the  tree  table 
had  been  generated.   The  first  node  in  the  lANODA  array  of 
the  link  table  was  found  in  the  lANODB  array  and  the  sub- 
script of  the  node  substituted  for  the  node  number  in  the 
lANODA  array.   For  example,  if  the  first  node  number  in  the 
lANODA  array  was  2,100  and  was  located  in  the  557th  location 
of  the  lANODB  array,  557  was  substituted  for  2,100  in  the 
lANODA  array.   This  same  type  of  substitution  was  made  for 
each  node  in  the  lANODA,  IBNODA,  lANODK,  and  IBNODK  arrays. 
This  step  reduced  the  necessary  search  time  by  about  one 
half.   A  similar  substitution  was  not  made  in  the  tree  table 
node  arrays  because  the  original  node  numbers  were  needed  in 
Expanded  Trees . 

A  final  initialization  step  was  that  of  setting  both 
the  lANODC  AND  LABELC  arrays  to  zero. 

The  source  centroid  and  the  forward  centroid  links  were 
then  read  in  and  the  appropriate  changes  made  in  the  tree 
table  and  the  reached  nodes  placed  in  C  storage. 

Subroutine  RANK  searched  for  the  minimum  LABELC  dis- 
tance.  The  comparison: 

LABELB .  >  min  LABELC.  +  LABELA . . 
^    jeS        3  ^^ 
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was  made  with  label  and  node  number  changes  performed  as  re- 
quired by  the  algorithm.   The  tree  table  was  completed  when 
set  S  (or  the  lANODC  array)  v/as  empty.   The  tree  table  was 
then  written  on  tape  for  use  in  Expanded  Trees,  the  procedure 
then  continuing  for  the  balance  of  the  source  centroids . 
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r^MALL  TRr  f  «; 


f.FxFrij 
SI RjnR 
SIBFTr 


TF 


I  M  JOR 


SMAP 


NDI  AN 

PFWTMH  8 

niMFNSION  I  ANnnA(  36la)  ,  IMNnnA(  ^ftl8),LA(^FLA(^6  1«), 
1  I  AMOnPl  14  77  )  ,  lONODRf  1  4  ■  7  )  ,  L  A  RFI.P  (  1  A  7  7  )  ,  [  A  NOnr  (  ISO  )  ♦ 
1LABFLC(1?>0),  I  ANO0K(  687  )  ,  IBN0DK(68  n.nLINK{'^fl7) 

COMMON    iANonc»LARFLr,v;rNon«i^RAM 

RFAL  LABP'LA,LAwFLR»LARFLr,LAHrLn,LAP,FLS 
PFAn('^,in)LINicQ,NIX,lvMy,kTN;on 

10  FORMAT ( A  I ^ ) 

PFAni  =^ ,  1 1  )  (  I  ANnnR(  i  )  ,  t  =  1  .6  ) 

1 1  FORMAT (6  I S ) 

RE  AD  (  5,30)  (  1  ANOL'IC(  I  )  ,  IBNOD<  (  1  )  .DL  I  NK  (  I  )  ♦  I  =  1  ♦  6P  7  ) 
RFAn(S,30)  (  I ANODAf  I  )  , IPNODA (  I  )  .L ARFl A (  I  )  ,  I  =  1  ,  L  I  NK?  ) 
30  FORMAT  (  ?  I  S  ,F=..0  ) 
J  =  7 

r>n  1  "^  T  =  1  ,4^^  1 
IANOnB(J)=IXNIX+T 

15  J=J+1 

I ANODR( J )  =  1 ANOOA (  1  ) 
no  1?  T  =  ]  ,L  [Mtf  <^ 

IF{  I  ANODA (  I  )  .FO.  1 ANOnp ( J  )  )nO  in  12 
J  =  J+1 

I  ANOn)R{  J  )  =  I  A^40^A  (  I  ) 
1?  r  ON  T  INI  IF 
Mnnr  c  =  J 
WRITF  (6.  16  )  1  ANODAILIMK-;)  ,  I  R  NOD  A  (LINK^;  )  ,  I  A  NOOB  (  NODF=.  )  .NODFS.L  IN<<^ 

16  FORMAT ( IHO.5 17) 

IF(  IANO0R{N00ES  )  .NF.4388 )G0  TO  ^00 
I  F(  IRNOOA  (L  INK";  )  .NF.P157  )GO  TO  ^On 
I  F(NOnP'^.NP.  1477  )r,0  TO  300 
f^O  f,0     1  =  1  ,lf  TNOn 
r  A^|n^^  (  n  =1^ 
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no  A?  1  =  1.1  INI^'^ 
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^-O  TO  4  7 
44  rOMTTMIlF 
4?  COMTIMMF 

no  46  T  =  1  .L  I  N<=: 

no  LP.     J-  1  ,NOOFC 

IF  (  IRNOOAf  I)  .N^.  lANO'-ri  (  J)  )r-.  ^   j   .^g 

I  RNOriA  (  I)  =  J 
GO  TO  4  6 
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rONTI  Nl  '~ 
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I  F  (  [  AMOuK  (  I)  .NF.  I  ANOnR  (  <  )  )  r,o  TO  502 

I ANOHK ( J ) =K 
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no    sof,   K  =  1  .NnnP"^. 
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hOU    CnNTINIIf^' 
lOCO     PFfln(  S  ,71))  IHOMf  »KRI '^P 
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T  ^|  o  r  V  =  I 

ro    TO    1  1  "^ 

TOO   roMTiMfir 

11  r    LA"PLR(!f\ioi^Y)=o.o 
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B2 .   Expanded  Trees 

The  mechanics  of  Expanded  Trees  will  be  outlined 
briefly. 

The  first  tree  table  on  the  tape  generated  by  Small 
Trees  was  read  in.   Because  the  first  centroid  number  of 
this  tree  table  was  656,  the  arrays  of  the  tree  table  were 
read  into  the  expanded  tree  table  storage  area  beginning 
with  the  656th  location.   The  first  655  locations  were  al- 
located for  the  centroids  of  the  four  border  states  and  the 
last  99  for  the  Interstate  nodes  including  boundary  and 
centroid  linked  nodes  for  a  total  of  2,225  centroids  and 
nodes  in  the  tree  table. 

All  of  the  border  node  sets  and  distance  arrays  were 
then  read  in  except  for  the  centroid  distance  arrays  DISTMC, 
DISTLC,  DISTKC,  and  DISTIC.   These  were  so  large  they  were 
stored  on  tape  and  read  into  equivalenced  storage  areas  as 
needed. 

All  node  numbers  were  then  changed  to  the  subscript  or 
location  of  that  node  number  in  the  lANODB  array  of  the  tree 
table  just  as  in  Small  Trees. 

After  these  and  other  initialization  steps  had  been 
completed,  the  first  tree  table  produced  by  Small  Trees  was 
re-read.   The  minimum  path  routes  from  border  node  to  border 
node  were  then  determined.   For  instance,  the  path  from  a 
node  in  H  to  a  node  in  K,  C,  or  F  was  determined  as  follows 
(refer  to  Figure  21) .   If  i  is  the  node  in  H  and  j  a  node 
in  H,  C,  or  F,  the  minimum  path  is  determined  as: 
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LABELS.  =  min  (LABELB .  +  DMBLLL . .) 
1  J  ID 

This  procedure  was  performed  in  the  program  by  a  series 
of  Do-Loops.   The  centroid  distance  arrays  were  then  read 
one  at  a  time  and  the  minimum  path  for  each  determined  by 
the  same  procedure. 

The  completed,  expanded  tree  table  was  then  written  on 
tape  and  the  next  small  tree  table  read. 
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EXPANDFD  TREES 

trXECUIF         IBJOB 
SIBJOR 

SIRFTC    FXTREE       SMAP 
ftrwiiNiD    8 
RFKirgO    9 
KEWIMD     10 

CI  MENS  ICN    lANODB (2225), ICN0DB(222  5) ,LABELB(222  5),DISTMC(l8,125)t 
1DISTLC(10,167)  ,niSTKC(  13,ll'i),DISTIC(25,2^9),DISTNC(280), 
INNATNni 280)  ,MACENO(  280) ,OMHLLL( 6,  10) ,MBLN0D(6) , DKOLLL ( 2 , 30 ) , 
ICLPKKK  (2,13  ),KDLN0D(2)  ,DIBKKK(3,13),DKBIII(3,25),  IBKNODO)  , 
1MM\'0UE(  18),LLN0CE(30),KKNCDE(  13)  »  I  I  NODE  (25)  .NODE  (30),JBRCEN(6), 
lNNATi\D{5|  .DISTMiNC),  18)  ,LNATN0(4)  .DISTLNC,  30),KNATND(2)  , 
lf:lSTKN(2,13),INATND(  12),DISTIN(  12,25) 
FQUIVALENCE( DISTIC(  1,  1  ) ,DISTKC(  1,1) ,OISTLC( I , 1 ) ,D I STMC (1 , 1 ) ) 
REAL    LABELB 

RE AD( I C, 52) I  HOME, NODES 
5  2    Fnl<^'AT  (  215) 

READ(10,51)(  IANOCB(  I ),  IBNODB( I ),LABELB( I), 1=650,2126) 
51     FOKMAT(  10(  2K,F5.0)  ) 
REKIimD     10 
IC     12     1=1,655 
12     IANnCB(  I)=  I 

RrAn(5,15)(IAM0DB( I), 1=2127,2225) 
15    FORMAT  (  201^.) 

RFAC(5,  15) (MMNODE(  I  ), 1  =  1,  18) 
REAC(5,l5)(LLNnCE{ I), 1=1,30) 
K£Ar)(5,15)(KKNQDE(I),I  =  l,13) 
READ(5,15)(IIMOOE{l),I=l,25) 
CC  20  1=1,18 
CO    21    J=1107,2225 

IF(  MMNODE(  I  ).NE.IANODB{  J)  )r,C    TO    21 
NMNOOEi I )=J 
GO    TU    20 
2  1    CONTINUE 
20    CONTINUE 

CU    21     1=1,30 
CO    23    J=1107,2225 

IF{ LLNODF(  I  ).NE.IANODB( J)  )G0    TO    23 
LLNOuEi  I  )=J 
CO    TO    27 
23    CONTINUE 
22    CONTINUE 

CO    2^     1  =  1,  13 

DO    25    J=1107,2225 

IF( KKNCnE( I ).NE.IANODB( J) )G0    TO    25 

KKNOUF (  I  )  =  J 

CO     TO    2^ 

25  CONTINUE 
2^    CCNTINUF 

CO    26     1=1,25 
CO    27    J=1107,2225 

IF( IINCDEI I ).\E.IANODB( J) )G0    TO    27 
IINODEI  I)=J 
GO    TO    26 
27    CONTINUE 

26  CONTINUE 

L'O     130     1  =  1,6 
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RrAn(5 

FORMAT 
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CC     132 

R[;At)(5 

RT  An(5 

RTAC(5 

CO     134 

READ(5 

KFAD(5 

UFAC(5 

CO     30 

CO     3  1 

IF(  MliL 

fTSLNQD 

GO     T3 

ClKJTlN 

CONTIN 

cn    32 

CO     3  3 
IF(  KBL 
KBLNOC 
CO    TU 
CONT  IM 
COigTIN 
CO     34 
CO    3  5 
IF(  I  UK 

iPKrjuD 

GO     TO 
CCMTIN 
C(!f>JTIN 
CO    4U 
READ(5 
REAC(5 
CO     41 
K  c  A  r,  (  5 
Kt=  Ar(  5 
CO    4^ 
HFAC(5 
XFAD(5 
CO    4  3 
HFAD15 
t<EADl5 
CO    44 
CO    45 
IF( MMA 
r'NATND 
GO     TO 
CONTIN 
CONTIN 
CO     46 
CO    4  7 
IF(  LNA 
LNATND 
GO     TO 
CONTIN 
CONTIN 
cn    46 
cn    49 
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.54 
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(20 
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.54 

1  = 
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1  =  1 
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(  I) 
30 
UE 
UF 
1  =  1 
J=l 
NOD 
(  I) 
3? 
UF 
UE 
1  =  1 
J=l 
NOD 
(  I  ) 
34 
UE 
UE 
1  =  1 
.53 
,54 
1  =  1 
,53 
.5^ 
1=1 
,53 
.54 
1  =  1 
.53 
.54 
1  =  1 
J  =  l 
TND 
(  I  ) 
44 
UE 
UE 
1  =  1 
J  =  l 
TND 
(  I) 
46 
UE 
UE 
1  =  1 
J=l 
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)  (DISTLN(  I .J)  ,J=1,30) 

,2 

)KNATNO( I ) 

)( niSTKN(  I, J)  ,J  =  1,13) 
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IF  (  KNATND(  I  )  .^4F.  I  ANODB(  JMGn    TO    <f9 
KNATND( I )=J 
GO     TO    ^fl 
CONT  INUn 
CHNTINUC 
CO    60     1  =  1,  12 
CO    hi    J=1107,2225 

IF(  InATND( I  I.ME. lANODBJ J)  )G0    TO    61 
INAT'nID(  I  )=J 
GO    TD    60 
CONIT  INUF 
CONTINUE 

KCAD(5,57)  (NNATND(  I  ),NACENO(  I),DISTNC(I),I  =  1,280) 
FORMAT (2I5,F5.0) 
CO    62     1=1,280 
CO     63    J=1107,2225 

IF( NNATNDl I ).NB. lANOCBI J) )G0    TO    63 
NNATNOI I )=J 
GO    TO    62 
CllNiT  INUE 
CONTINUE 
CO    6^     1=1,280 
CO    6b    J=1107,2225 

IF(NACENID(  I  ).NE.  IAN0DB(  J)  )G0    TO    65 
NACEi^D(  I  )=J 
CO    TU    64 
CONT  INUE 
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FORMAT (614) 
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KFAPC 10,52 )  I  HO MF, NODES 
CO     101      1=1,2225 
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LABFLBI  I  1=9999. 
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JTKX=649+J 

It)NOCB(JTRY)  =  I  BNODBI  JTRX) 
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lANOCB ( JTRX  )=JTRX 
ITCHY=0 
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CO     112     J=l,30 
IXC=LLNODE( J ) 

TOT AL=LABELB( IXC)+DMBLLL(  I  ,J) 
IF( TOTAL. OF. LABELBI IXB) )Ga     TO     112 
inNOLB (  IXB)=  lANOCBI  IXC) 
LAHELB(  IX6)=T0TAL 
ITCHY=ITCHY+1 
CONTINUE 

IF(  ITCHY. NE .0)G0    ID    201 
lTCHy=0 
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.0)GO    TO    203 

5 

I  ) 

18 

J) 

B( IXC)+DISTMN( I, J) 

.LABFLB( IXR ) )G0  TO  120 

=  IANODR(  IXC  ) 
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REACIS 
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LABFLB! I 
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1?5 
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NODE( J) 
LABFLB(  I 
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ISTMC( J,I ) 
IXB))GO  TO 
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XC  )+0 
BFLB( 
NOCB( 
TAL 

MNIX 


(  I  ) 

TLC (  I,  J),  J  =  1,MNIX) 


ISTLCI J,  I) 
IXB)) GO  TO 
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XC  )*D 
BELB( 
NODB( 
TAL 

MNIX 


(  I  ) 

TKCt  I ,J)  ,J=1,MN1X» 


XC)4-D 
BELB  ' 
NODB( 
TAL 


I STKC( J,  I  ) 
IXB  )  )G0  TO 
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150 
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ISO 


152 


15A 


301 
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300 


»I)ATA 


rPMTlNUF 

-I  ■■.•   (t^.S^  )MIX,MNIX 

DO     liZ     I  =  1,NU 

READ(8,53)NnDE ( I ) 

KEAD(8,55)(0ISTIC( I,J) ,J=l,MNIX» 

CO     IbA     1=1,2^9 

IXr.=  l-t-'»06 

CO     lt)4    J=l,25 

lXC=IINnoF( J ) 

TOT AL=LABELB( IXC)*DISTIC( J,I ) 

IF( TOTAL. GE.LABELB( 1X8 ) )G0     TO    1 5A 

IBNOCBl  IXB)  =  lANQDBl  IXC) 

LABFLB (  IXB)=TOTAL 

CONT  INUE 

RFWIND    R 

V.P,ITE(  6,301)  IHOME 

FORMAT (  IHI,  lOX,  I2HH0ME     NODE     IS,  15///) 

WKITE( 9,52)  IHOME,NNOnES 

WRI TE(  9,302)  {  IANODB(  I  )  ,IBNODB( I ),LABELB(l )  , I  =  1,NN0DES) 

FORMAT (  10( 2I4,F5.0)  ) 

WRI  TE(  6,303  ) 

FORMAK  IHO,  20X,22HTREE     TABLE    ON     TAPE     262////) 

1F(  IHOME. EG. 727)G0    TO     300 

GO    TO     1000 

REWIND    8 

REW  no    9 

RFWIND     10 

STOP 

END 
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B3.   Trades 

The  nnechanics  of  Trades  will  be  outlined  briefly. 

After  all  the  data  had  been  read  in,  all  node  numbers 
were  changed  to  the  subscript  of  their  location  in  the 
lANODB  column  of  the  expanded  tree  table.   This  included  both 
the  lANODB  and  IBNODB  columns  of  the  tree  table. 

The  cumulative  totals  of  ITDF's  calculated  were  kept 
for  the  2,510  centroid  and  highway  links  making  up  the 
Indiana  network.   Each  link  was  described  by  the  subscripts 
representing  the  two  end  nodes. 

The  path  between  the  source  centroid  and  another  cen- 
troid was  traced,  the  path  being  stored  as  a  series  of  num- 
bers from  1  to  2,510,  each  representing  one  of  the  Indiana 
centroid  or  highway  links.   The  travel  desire  factor  was 
calculated  and  assigned  to  each  of  the  links  on  the  path. 
When  all  factors  had  been  calculated,  the  population  of  the 
source  centroid  was  changed  to  zero  and  the  next  source  cen- 
troid with  its  associated  expanded  tree  table  read. 

The  calculation  of  the  travel  desire  factors  were  sub- 
ject to  the  following  restrictions: 

1.  Only  cities  of  greater  than  5,000  had  interactions 
with  Interstate  cities. 

2.  Cities  of  less  than  1,000  had  interactions  only 
with  cities  within  150  miles  of  their  location. 

3.  Cities  between  1,000  and  5,000  had  interactions 
with  other  cities  of  less  than  5,000  only  if  within  300 
miles  of  their  location. 
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4 .   The  population  of  the  larger  city  of  a  pair  of 
cities  was  not  allowed  to  be  greater  than  10  times  the  popu- 
lation of  the  smaller  city. 

These  restrictions  were  imposed  in  an  attempt  to  over- 
come some  objections  to  a  gravity  model  approach  to  travel 
synthesis  as  outlined  in  an  earlier  chapter,  "Study 
Procedure,"  the  principal  objection  being  that  use  of  this 
approach  allows  a  virtually  unlimited  trip  generation  capa- 
bility of  each  centroid.   The  only  basis  for  the  numerical 
value  of  the  imposed  limitations  was  that  they  appeared 
reasonable . 
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3. 17*P0PULA)  )POPHOM=3.17*POPULA 
POPHOM 
OISTER 
DISTER 
1 


HAO( JJ)+SIGMAA 
HAT( JJ)+SIGMAB 
HAF( JJ  )  +  SIGMAC 


)**2.LT.500C.  )G0  TO  2'59 


196 


1200 


222 
223 


22'^ 


230 


2'tl 


262 


26  1 


260 

2A0 


282 
283 


290 


CO  220 
IF(  I13N 
CO  222 
IF( MUD 
IF(  PUP 
GO    TO 

co'vjt  in 

I^JCRMT 
CO  7  2'^ 
KruNT= 
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IF( IXC 
IX=IX+ 
NUKPEK 
CO  10 
JX= JX+ 

ir(  ixR 

GO  TO 
CON  TIN 
CO.MT  IN 
11=1 
[SAHE  = 
CO  280 
IF(  lEN 
IF( IX. 
CO  2  32 
IF( MOO 
KK=K 
GO  TO 
CONTIN 
CO  23^ 
IF( NOD 
POPHCM 
POPULA 
POPFX( 
CISTER 
SIGMAA 
SIGPAB 
SIGVAC 
IKNKNT 
CO  290 
JJ  =  NUM 
ALPHAO 
ALPHAT 
ALPHAF 


1  =  2 
OQB( 

N=l 
EX(N 
FX(N 
223 
UE 
=  1 

K  =  I 
K 

(K)  = 
00P( 
=  IBN 
6,13 
3C0 


127,2212 

1  )  .EO.OIGO  TO  220 

,  181 

).NE.I)GO  TO  222 

) .EO.O. )G0  TO  220 


,  100 

INCRMT 

INCRMT).EQ.O)GO  TO  230 
ODB(  INCRMT ) 
't  )KOONT 


K=2f KOUNT 
AKNU( K-1 ) 
AKND(  K  ) 

•  LT.IXOGO     TO    241 
AKND(K) 
AKNU( K-1 ) 

.GT.2067 )G0     TO    2A0 
NOEXR 

J  =  JEM, J  IM 
.NE.NODA( J ) )G0     TO    260 

.NE.NODB ( JX) )GC    TO    261 

1 

(  IX  )=JX 

240 

1 

•FQ.NODAI JX ) )G0    TO    262 

240 

OF 

OF 

IBNODB( I ) 

N=I  1,2212 
tDR(N).ME. ISAME)GO     TO    2  80 
FQ.0(G0     TO     280 

K=l,  181 
FX(K ).NE.N)GO     TO    282 


283 
UE 

K  = 
EX( 
=  P0 
=  P0 
K)  = 
=  LA 
=  P0 
=  SI 
=  SI 
=  IK 

M  = 
BER 
(  JJ 
(  JJ 
(  JJ 


KK, 
K  )  . 
PIN 
PEX 
0.0 
BEL 
PUL 
GMA 
GMA 
NKN 
1.  I 
(M) 
)=A 
)=A 
)  =  A 


18  1 

NE.NJGO 
( IHOME) 
(K  ) 


TO    280 


B(N)+01STEX(K) 

A«POPHOM 

A/CISTER 

B/DISTER 

T+1 

X 

LPHAO( JJ)+SIGMAA 
LPHAT( JJ )+SIGMAB 
LPHAF( JJ )+SlGMAC 
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28^. 

?80 

?2n 

299 

390 


1100 
AOO 

410 

A20 

1050 

1060 

300 

99  3 

$IBFT 


1^.2 

!<»'> 
1A5 
146 

147 

148 

14') 
15  1 
152 


CONTIN 
CON  TIN 
CHNTIN 
WR I T  E ( 
FORKAT 
P(1P  n( 
IF( IHG 
GO    TO 
WRITE( 
FORMAT 
WRIT6( 
FORMAT 
WKITE( 
FORMAT 
DO     105 
WRI TE( 
FORMAT 
REWIND 
RF  W  I  NO 
WRITE( 
FORMAT 
STOP 
END 
C    SUBPl 
SUI3R0U 
COMMCN 
J1M=25 
1F(  IXB 
IF(  UB 
JEM=1 
JTM=10 
GO    TO 
1F(  IXB 
JEM=10 
GO    TU 
IF( IXB 
JEN'  =  20 
GO    TO 
IF(  IXB 
JFM=jO 
GO    TO 
IF( IXB 
JE^'=40 
GO     TO 
IF(  U8 
JEM  =  'jO 
GO    TO 
JEM=60 
GO     TO 
IF( IXB 
IF(  TXB 
JFM=70 
GO     TO 
IF( IXB 

jrM=8o 

GO     TO 
1F(  IXB 
JEM=90 

GO    TO 

ir(  IXB 
jrf=io 


UF 

UE 

UE 

6,39 

(  IHO 

IHOM 

ME.F 

ICOO 

6,40 

(  IHl 

9,41 

(6(  1 

6,42 

{  IHO 

0     1  = 

6,10 

(  HO 

8 

9 
6,99 
(  IHl 


0 ) IHOME 

, 10X,6HIH0ME=, 15) 

e)=o.o 

0.727)GO    TO     llOO 

0) 

,5X, 31HINTERCITY     TRAVEL    DESIRE    FACTORS//) 

0)  ( ALPHAOI  n  ,ALPHAT( 1  ), ALPHAFI  11,  1  =  1, 2510) 

PE20.9)  ) 

0) 

,20X,18HFACT0RS    ON     TAPE     61) 

1,2510 

60)1  ,  ALPHAOI  I  ),  ALPHAK  I  )  ,ALPHAF(  I  ) 

,  3(  1PE20.9) ) 


3  )  IKNKNT 

,  10X,7HIKNKNT  =  ,  I  10) 


TINF  INOEXR 

IXB, J  EM, JIM,  I  HOME, NODES 
10 

.GT.1106)G0  TO  148 
.GT.  721  )GG  TO  142 

0 

150 

.GT.  789)G0  TO  143 

1 

150 

•GT.  865)G0  TO  144 

1 

150 

.GT.  927)G0  TO  145 

2 

150 

•GT.  986)G0  TO  146 

3 

150 

.GT.1049)GO  to  147 

2 

150 

2 

150 

.GT.1396)G0  to  156 

.GT.1148)G0  to  149 

3 

150 

.GT.I188)G0  to  151 

2 

150 

.GT.1222)G0  TO  152 

1 

150 

.GT.1262)G0  TO  153 

00 


198 


15-!. 

155 
156 

157 

158 

159 

163 

164 

165 

166 

167 

16B 


169 

150 


CD     TL)     150 
IK  IXR.GT 
J[:f:=llOO 
Cn    TO     150 
IF(  IXB.GT 
JFM=il99 
GO    TD     150 
JCM=12'5fl 
GO     TO     150 
IF{  UC.GT 
JF»'=U04 
cn    TO     150 
IF( IXe-GT 
JFM  =  150  } 
cn    TO     150 
1F(  IXB.GT 
Jt:M  =  l605 
GO    TO     150 
IF( IXB.GT 
JHM=1709 
cn    TO     150 
IF( IXB.GT 
JEf'=180a 
G(l     TO     150 
IF(  IXB.GT 
jrf'  =  1908 
cn    TO     150 
IF(  UR.GT 
JFM=2008 
GO    TO     150 
IF( IXB.GT 
JEM  =  21  n 
GO     TO     150 
IF( IXB.GT 
JFr'  =  2210 
GO     TO     150 
IF(  IXB.GT 
JFM=2311 
GO    TO    150 
JEM  =  24  11 

kftukn 

FND 


.1311  )G0    TO     15A 


ISS'.JGO    TO    155 


1AA2)G0    TO    157 


.1486)G0    TO    158 


15A0)GO     TO    159 


1583)GO    TO    163 


1622 )G0    TO    16A 


I663)G0    TO    165 


.1732)GO    TO    166 


1822)G0    TO    167 


1891  )G0    TO     168 


1984)00    TO    169 


tr.ATA 
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B4.   Border  Trees 

The  mechanics  of  Border  Trees  will  be  outlined  briefly. 

After  the  distance  arrays  and  node  sets  had  been  read 
in,  the  lANODB  array  of  the  tree  table  was  generated  and 
all  node  numbers  changed  to  the  subscript  of  their  location 
in  the  lANODB  array. 

The  distance  arrays  were  called  in  the  sequence  listed 
in  Table  Bl  depending  upon  the  border  state  in  which  the 
source  centroid  was  located.   Table  Bl  shows  the  operation 
performed  with  each  distance  array.   For  instance,  in  the 
sequence  for  a  source  centroid  in  Michigan,  the  first  time 
DISTML  is  used,  every  node  in  B  is  given  the  opportunity  to 
be  the  next  node  on  the  minimum  path  node  sequence  to  the 
source  centroid  from  every  node  in  C.   The  second  tim.e  DISTML 
is  used,  every  node  in  C  is  given  the  opportunity  to  function 
similarly  for  every  node  in  B. 

The  nature  of  the  operation  is  given  by  the  entries  in 
the  "From  Node"  and  "To  Node"  columns  of  Table  Bl  for  each 
distance  array  listed.   The  operations  performed  on  each 
distance  array,  except  for  the  first  one,  are  repeated  until 
no  label  changes  for  any  node  are  made  in  an  entire  sequence. 

At  this  point,  the  tree  table  and  the  minimum  path  node 
sequences  from  each  centroid  to  the  source  centroid  were 
written  and  the  next  source  centroid  was  read. 
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Table  Bl, 

Distance 

Arr 

•a^ 

'  Sequence  in  Bor 

der  Trees-By  State 

Michigan 

Ohio 

Distance 

Node 

Sets 

Distance 
Arrays 

Node 

Sets 

Arrays 

From 

'o 

From 

To 

Node 

Node 

Node 

Node 

DISTMC 

B,H 

M'* 

DISTLC 

F,C, 

,H 

L* 

DISTML 

C 

B 

DMBLLL 

H 

F,C,H 

DMBLLL 

F,C, 

,H 

H 

DKBLLL 

F 

F,C,H 

DISTML 

B 

C 

DMBLLL 

F,C, 

,H 

H 

DMBLLL 

H 

F, 

c, 

H 

DKBLLL 

F,C, 

,H 

F 

DKBLLL 

F 

F, 

c, 

H 

DISTML 

B 

C 

DISTMI 

E 

B 

DISTLI 

E 

C 

DISTLI 

E 

C 

DISTLK 

D 

C 

DISTMK 

D 

B 

DISTML 

C 

B 

DISTLK 

D 

C 

DISTMK 

D 

B 

DKBIII 

G 

E, 

G 

DISTMI 

E 

B 

DISTKI 

D 

E 

DLBKKK 

G,D, 

rF 

F 

DLBKKK 

G,D, 

,F 

F 

DLBKKK 

F 

G,D,F 

DKBIII 

E,G 

G 

DISTLK 

C 

D 

DISTKI 

E 

D 

DISTMK 

B 

D 

DIBKKK 

G 

G, 

,D, 

,F 

DISTKI 

E 

D 

DIBKKK 

G,D, 

,F 

G 

DIBKKK 

G 

G,D,F 

DLBKKK 

F 

G, 

,D, 

.F 

DISTMI 

B 

E 

DKBLLL 

F,C, 

,H 

F 

DISTLI 

C 

E 

DISTLK 

C 

D 

DISTKI 

D 

E 

DISTLI 

C 

E 

DKBIII 

G 

E,G 

DISTMI 

B 

E 

DKBIII 

E,G 

G 

DISTMK 

B 

D 

DIBKKK 

G,D, 

,F 

G 

DISTLC 

L* 

F, 

■  c, 

,H 

DISTKC 

K* 

G,D,F 

DISTKC 

K* 

G, 

,D, 

,F 

DISTIC 

I* 

E,G 

DISTIC 

I* 

E, 

,G 

DISTMN 

N 

B,H 

DISTMN 

N 

B, 

,H 

DISTLN 

N 

F,C,H 

DISTLN 

N 

F, 

,C, 

,H 

DISTKN 

N 

G,D,F 

DISTKN 

N 

G, 

■  D, 

,F 

DISTIN 

N 

E,G 

DISTIN 

N 

E, 

,G 

DISTNC 

N* 

N 

DISTNC 

N* 

N 

DLBJLJ 

C 

C 

DMBJMJ 

B 

B 

DMBJMJ 

B 

B 

DLBJLJ 

C 

C 

DKBJKJ 

D 

D 

DKBJKJ 

D 

D 

DIBJIJ 

E 

E 

DIBJIJ 

E 

E 

RETURN  TO  DISTML 


RETURN  TO  DMBLLL 


Centroids 
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Table  Bl. 

Continued. 

Kentucky 

Illinois 

Distance 

Node 

Sets 

Distance 
Arrays 

Node 

Sets 

Arrays 

From 

Tc 

) 

From 

To 

Node 

Node 

Node 

Node 

DISTKC 

G,D, 

,F 

K'' 

DISTIC 

E,G 

I* 

DLBKKK 

F 

G, 

D, 

F 

DKBIII 

G 

E,G 

DIBKKK 

G 

G, 

D, 

F 

DKBIII 

E,G 

G 

DLBKKK 

G,D, 

,F 

F 

DISTKI 

D 

E 

DIBKKK 

G,D, 

,F 

G 

DISTLI 

C 

E 

DISTLK 

C 

D 

DISTMI 

B 

E 

DISTMK 

B 

D 

DIBKKK 

G,D, 

-F 

G 

DISTKI 

E 

D 

DIBKKK 

G 

G,D,F 

DKBLLL 

F,C, 

rH 

F 

DISTKI 

E 

D 

DISTLK 

D 

C 

DISTMK 

B 

D 

DISTML 

B 

c 

DISTLK 

C 

D 

DISTLI 

E 

c 

DLBKKK 

F 

G,D,F 

DKBLLL 

F 

F, 

^c, 

H 

DISTMI 

E 

B 

DMBLLL 

H 

F, 

•  c, 

H 

DISTMK 

D 

B 

DKBIII 

E,G 

G 

DISTML 

C 

B 

DKBIII 

G 

E, 

,G 

DISTLI 

E 

C 

DISTKI 

D 

E 

DISTML 

B 

C 

DISTLI 

C 

E 

DISTLK 

D 

C 

DISTMI 

B 

E 

DMBLLL 

H 

F,C,H 

DMBLLL 

F,C, 

rH 

H 

DKBLLL 

F 

F,C,H 

DISTMK 

D 

B 

DLBKKK 

G,D, 

,F 

F 

DISTMI 

E 

B 

DMBLLL 

F,C, 

,H 

H 

DISTML 

C 

B 

DKBLLL 

F,C. 

,H 

F 

DISTIC 

I* 

E, 

,G 

DISTMN 

N 

B,H 

DISTMN 

N 

B, 

,H 

DISTLN 

N 

F,C,H 

DISTLN 

N 

F, 

■  c, 

,H 

DISTKN 

N 

G,D,F 

DISTKN 

N 

G, 

■  D, 

,F 

DISTIN 

N 

E,G 

DISTIN 

N 

E, 

,G 

DISTNC 

N* 

N 

DISTNC 

N* 

N 

DIBJIJ 

E 

E 

DKBJKJ 

D 

D 

DKBJKJ 

D 

D 

DLBJLJ 

C 

C 

DMBJMJ 

B 

B 

DIBJIJ 

E 

E 

DLBJLJ 

C 

C 

DMBJMJ 

B 

B 

RETURN  TO  DLBKKK 


RETURN  TO  DKBIII 


Centroids 


202 


enKDLR  TREFS 

NICHIGAN  AS  BORDER  STATE 

$FXECUTE        IBJ08 

tIBJOO 

SIRFTC    EXTCFN       SMaP 

Clf'EiMS  ICN  I  ANODBI  360),  I  BNODB  (  360  )  ,  L  ABENB  (  3  60 )  ,MCENTR(  39)  ,MMNODE(  IB 
1),DISTMC(18,39),MMNOBC(18),LLNODF(30),LLNOBE(3  0) ,LCFNTR{65) .DISTLC 
1(30,65)  ,KKMnUE(  13)  ,KK;^J0BF(  13),K,CENTR(34),DISTK:(13,3A)  ,I  I  node  (25)  , 
in\inBF(25),lCENTR(98),DlSTlC(25,98)  ,  MB  J  NOD  (  11)  ,MBJNOB(  11)  ,LRJNDD(2 
12).  LBJiMOBI  2?)  ,KBJNQD(  8  )  ,KBJNOB(  8),  I  BJNOD(  22)  ,  IBJNaB(  22), DIST  ML  (lit 
1  22  ). DISTMK(  11,  8), DISTMK  11, 22  ),DISTLK(22, 8), DISTLM  22, 22), DISTKK  8 
1,2?)  ,MBLNOD( 6),MBLN0B(  6),KBLN0D(2)  ,KRLNnB( 2) ,LBKNnD(2)  ,LBKN0B(2) , I 
1HKMOO( 3),  IRKNOb( 3) , KB  I  NOD ( 3) , KB  I  NOR (3 ) , DMBLLL ( 6 , 30 ),DKBLLL(2,33),D 
1LBKKK(2,13),DIBKKK(3,13),DKBIII(3,25) 

CIMF.\S  ICN  MNIATND(5  ),M^MTNB(  5),0ISTMN(5,  18)  ,LNATND(4)  ,LNATNB(^),DIS 
ITLM(4,30),KNATND(2),KMATNB(2),DISTKN(2,13),INATND(12),INATNB(12) 

ClriF^jSInN  DrSTIN{12,25),NNATND(280)  ,NNATNB(280)  ,N  AC  END  (280)  .NACENB 
1(280), CISTNC(2 80), LENGTH (360) 

CI  MENS  ICN    DMBJMJt  18,  18 ),DLBJLJ( 22, 22)  ,DKBJKJ(8,8),DIBJIJ(22,22) 

COMMON     IHOME.NODES, lANODB, I BNODB , LABENB , NI X,NOX,KIX 

REAL    LABENB 
COMMENT       LABENB    EQUIVALENT    LABELB 

INCFX=0 

NnC.ES=360 

MX  =  198 

NCIX  =  NI  X+1 

KIX=NIX+P6 

IXN IX    =    KIX    -     1 

IF(NIX .FO.l  )G0     TO     10 

READ('?,bl)(IANJ00B(I),I=2,NlX) 
10    RFAP(5,51)(IA^40^B((),  I=NOX,NODES) 

51  FnKMAT(20I^) 
3EAD(5,51)(MCENTR(I),I=1,39) 
CU     60     1=1,18 
RFAn{5,52)MMN0DE(  I  ),MN1X 

60    READ(5,56)(DISTMC(  I,J),J=1,MNIX) 

52  F0RMAT(2I5) 

56  FORMAT (  16F5.0) 

RFAD(^ ,51)  {LCENTR(  I  ),  1  =  1,65) 

CO  70  1=1,30 

RFA0(5, 52 )LLNGDE(  I  )  ,L\IX 
70  REAn(5,56)(DISTLC(I,J),J=l,LNIX) 

RFAD(5,5l)  (KCENTR(  I  ),  I =1, 34) 

Cn  80  1=1,13 

RFAO( 5,52 )KKN0DE( I  )  ,KNI  X 
8  0  REAn(5,56)(DlSTKC( I , J) , J=1,KN!X) 

RFAn(5,51)(lCENTR(  I), 1=1,98) 

CO  90  1=1,25 

REAn{5, 52 ) I INODE( I ) , INIX 
90  RCAn(5,56)(DISTIC( I,J),J=1,INIX) 

CO  100  1=1,11 

RFAOl 5,53)MBJNOD(  I  ) 

53  FnR^•AT(  15) 
WEAn(5,56)(DlSTML{I,J),J=l,22) 
RF A  0(5, 56)  (DISTMKI  I ,J)  ,J=1, 8) 

100  RLAD(5,56)(D1STMI(  I,J),J=1,22) 
cn  110  1=1,22 
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110 


120 


130 
54 


131 
132 
133 
134 

601 

602 

603 

604 

57 

605 

606 

607 

608 
59 

125 


141 

140 


RFAni5,53ILBJN00{  I  ) 
RFAD(5,56) (OISTLKt I , J) , J=l, 8) 
l<EAn(5,56)(  OISTL  I(  If  J).  J=lt22) 
LO     120     1  =  1, R 
Rt  An(5,53)KBJMOn(  I  ) 
KFAD(5,56)(DISrKI(  I,J),J=1, 22) 
KFA[)(5,51)(  IBJNnD(I),l  =  l,22) 
Cn     130     1=1,6 
l<FAD(5  ,53)MBLNOn(  I  ) 
I-;EAC(5,54)  (OMBLLK  I  ,  J)  ,  J=L,  30) 
Fi;Rr'AT(  ?0F'..0) 
LO     IJ 1     1=1,2 
l<PAri5,53)KRL'^OD(  I  ) 
KrAC(5,54) ( DKBLLL( I, J) , J=l, 30) 
rn     li2     1=1,2 
K[;Ar)(5,53)LBKM0n(  I  ) 
rttAD(5,5^)  (DLBKKK(  I , J)  ,J=1,  13) 
cn     133     1=1,3 
RFAn(5, 53 ) IBKNOOC I ) 
Kr\D15,54)(niBKKK(  I,J),J  =  1,13) 
cn     Ij't     1=1,3 
l<FAC15,53)KnilM0C(  I  ) 
READ(5,54)(DKHIII(I,J),J=1,25) 
CO     601     1=1,5 
RFAi::(5,53)MNATNn(  I  ) 
;<FAD15,54)(DISTMN(  I, J)  ,J=l,  18) 

1  =  1,4 

53)LMATNO(  I  ) 

ur Ar(s ,54) (nisTLN( I , J) ,j=i,3o  ) 

CO    603     1=1,2 

RFAIi(5,  53)KNATND(  I  ) 

HFAr)(5,54)(DlSTKN(I,J),J  =  l,13) 

CC    60^,     [  =  1,12 

KFAD(5,53) INATNCt I ) 

RCAn{5,54)(DISTlN(I,J),J=l,25) 

RCAn(5,57)  (NNA  T>4D  (  T  ),.NACEND(  I  )  ,  D  I  STNC  (  I)  ,  I  =  1  ,2  80  ) 

FniMAT ( 21 5, F 5.0) 

C(J    605      1=1,18 

RCAD(5,  59)(  Df'.HJMJI  I,J),J=1,18) 

C(J     606     1  =  1,22 

RrAn(5,59)(nLBJLJ(I,J),J=l,22) 

CO    607     1=1,8 

WFAn(5,59)(DKBjKJ( I,J),J=1,8) 

cn     608     1=1,22 

RC A r(^, 59 )(OI0JIJ(  I, J  ),J=1, 22) 

FORMAT (22F3.0) 

W«ITE(6,125)KBINnn(3),DKBIII(3,25) 

FHRN'AT  (  IHO,  I  5,  FIO.O) 

IF{ KalNODf 3).NF.4303)GO    TO     300 

IF(  Cl^D  11  l(  3,25  ).NE.  386.  )G0     TO     300 

rn     1^0     1  =  1,  18 

cn     141     J=iNUX,'>IODES 

IF( MMNCnEi I ) .NE. IAN0CB( J) )GQ     TO     141 

^MN0BF ( I)=J 

GO     TO     140 

CON  TINUG 

CnNTINUE 

CO    142     1=1,30 

CO     143    J=NOX, NODES 

IF(  LLNODEC  n.NE.  I  ANODBt  J)  )G0     TO     143 

nfiF(  I  )  =  j 


cn    M2 
HF^OI':), 


204 


v.^     .u     1^2 

1^3 

CONT  INUF 

1-^2 

CONTINUE 

CO     IA4     1=1,13 

CO     145    J=Nnx, NODES 

ir( KKNODE(  1  ).NE.  I&NODB( J)  )G0 

TO 

145 

KKNOOE ( I )=J 

GO    TO     1A4 

K5 

CP^JTI^4UE 

l^t* 

C(1NTINUE 

CO     146     1=1,25 

CO     147     J=NOX, NODES 

IF(  I  IN ODE (  I  ).NE.IANODB( J)  )G0 
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SIPFTC    SDCPl 
SUf'RUb 


736 


ao'f 


T*? 


800 


270 


290 


300 


NOBE( 

LABEN 

AL.GE 

(  IXB) 

(  IXB) 

TCH*1 

ITCHY 

UF 

HY.NE 

0 

1=1, 
JNOB< 

J  =  l, 
JNOB( 
LABEN 
AL.GF 
(  IXB) 
(IXB) 
TCh  +  1 
ITCHY 
UE 

HY.NF 
0 

1=1, 
JNOBI 

J  =  l, 
JNflBI 
L  AHEN 
AL.GE 
(  IXB) 
(  IXB) 
ICH+1 
ITCHY 

ur 

FY.NF 
0 

t  =  l, 
JNOBI 

J=l, 
JNOBI 

laben 

AL.GF 
(  IXB) 
(IXB) 
TCH+l 
I  TCHY 
UF 

FY.NF 
H.NE. 

K  =  l 
(K)  =  L 
6,270 
(  IHI, 
6,290 
(7(2X 
ATHb 
ICOO 


728 


730 


J) 

B(  IXC)+DMBJMJ(  I, J) 

.LABFNB(  IXB)  ) GO     TO 

=  IANODB(  IXC  ) 

=  TOTAL 

♦  I 

.0)G0     TO    809 

22 

I  ) 

22 

J) 

B(  I XC)  +  DLBJLJ(  I,  J ) 

.LABENB( IXB ) )G0    TO 

=  IA:>(0DB(  IXC  ) 

= TOTAL 

+  1 

.0)GO    TO    RC2 

a 

I  ) 

0 

J) 

B(  IXC)+DKRJKJ(  I , J) 
.LAPCNBI IXB ) )G0     TO 
=IANODB( I XC ) 
=TOTAL 

+  1 

.0)GO    TO    805 

22 

I) 

22 

J  ) 

n( IXC )+DIBJ  IJ(  I  ,  J) 

.LABFNBI I XB ) )G0     TO 

=  lANOCBI  IXC  ) 

=TOIAL 

♦  1 


.0)Gn     TO     ROA 

n)GO    T-)     11 CC 

NODES 

ABENBIK  ) 

)  IHOME,  ITER 

22H    THE     STARTING    NODE     IS    i I  5 ,5X  ,  II H I TF R AT  I ONS= # I  3 ) 

)  ( I ANOOBI  I  )  ,  IBNODBI  I  ),LABEN8 ( I)  ,  1  =  1, NODES) 

,21 3,F6.0) ) 


736 


7^2 


S»-'AP 
IINE    PATHS 
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CI  MEMS  ION     I ANODBI 360) , IBNODBt  36  0),LABFNB(360),  IBAKND(IOO) 

Cn^'MCN     I  HOME,  NODE  Sf  lANODB,  I  BNODB  ,  LABENB,  Nl  X.NOX.KIX 

REAL  LABENB,  IBAKDT 

WRI  TE( 6,460)  IHOME 
460  FOKMAT ( 1H0,6X,  16HTHE  HOME  NODE  I  S , I  5, 3X , 2 IHTHE  MINIMUM  PATHS  ARE/) 

CO  500  1=2, KIX 

IBAKDT=LABENB(  I) 

INCRMT  =  I 

CO       b20    K=1,100 

KnuMT=K 

IDAK^^D(K)  =  IANODB(  INCRMT) 

IF(  lENODBC INCRMT).EO.IHOME)GO    TO    550 

CO    5J0    L=NOX, NODES 

IF(  lONODBI  INCRMT ).NE. IANODB(L ) )G0    TO    530 

INCRMT=L 

GO     TU    520 
530    CONTINUE 
520    CONTINUE 
550    K0UNT=K0UNT+1 

inAKNC(KOUNT)= IHOME 

WRI TE( 6,560) IBAKCT, (  IB AKNO ( M  )  , M= 1 , KOUNT ) 
560    F0RMAT(1H    , 3X, F5. 0 , 3X , 201 5 ) 
500    CONTINUE 

RETURN 

END 
tDATA 
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VITA 
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